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THE THERMAL NEUTRON CAPTURE CROSS SECTION 
AND THE RESONANCE ABSORPTION INTEGRAL OF 
Th*? (LONIUM)! 


M. J. CABELL? 


ABSTRACT 


The thermal neutron capture cross section and the resonance capture integral 
of Th?*° have been determined by an activation method. Assumptions are made 
that the thermal capture cross section and the resonance capture integral of 
Co** are 36.5 barns and 48.6 barns respectively, that the half-life of Th?*° is 8.00 X 
10 years, and that Th®° behaves as a 1/v absorber in the subcadmium region. 
Based on these assumptions values of o2200 = 21.4+0.3 barns and of 937+32 
barns for the resonance absorption integral (from 0.5 ev to ~) have been ob- 
tained for Th, 

The half-life of Th?*' has been determined as 25.52+0.01 hours. 


INTRODUCTION 


Hughes and Harvey (1955) quote 27+2 barns for the thermal neutron 
capture cross section of Th?*° as determined by the pile oscillator method and 
35+10 barns for its capture cross section for pile neutrons as measured by an 
activation method, but no experimental details are available for these measure- 
ments. Jaffey and Hyde (1952) give two widely different values (61 barns and 
28 barns) for the thermal neutron capture cross section but the experimental 
uncertainties for each result are so large that neither can be preferred to the 
other. No measurements have been reported in the literature for either the 
resonance absorption integral of Th?” or the dependence of its absorption cross 
section on neutron energy. This paper describes activation measurements 
which have been made to determine the thermal neutron capture cross section 
more precisely; the resonance absorption integral has also been determined. 

The following nuclear reactions are involved: 


1Manuscript received March 3, 1958. 
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The experiments involved the irradiation in the NRX reactor of sub- 
microgram amounts of Th” (a) unshielded and (6) shielded with 0.025 in. 
cadmium. Cobalt monitors were attached to measure both the thermal and 
the epicadmium fluxes. After irradiation the thorium was purified by ion 
exchange and mounted on thin films for counting. Determinations of the 
absolute 6-disintegration and absolute a-disintegration rates of each source 
gave measurements of its Th! and Th*° content respectively. From these 
data and the known fluxes the effective capture cross section of Th**® for 
pile neutrons, its thermal capture cross section, and its resonance absorption 
integral were calculated. 

The method used is similar to that reported by Cabell, Eastwood, and 
Campion (1958) and possesses the advantage that chemical yields from the 
purification process, which vary with each experiment, need not be known. 


EXPERIMENTAL* METHODS 

(a) Starting Material 

Approximately 25 ug of Th**° in 20 ml of solution was available as a con- 
centrate from a Canadian pitchblende residue. The solution also contained 
macroquantities of rare earths and about 65 wg of Th*”. The thorium was 
purified in good yield and virtually free from all impurities by thenoyItrifluoro- 
acetone extraction (Barnes and Potratz 1954) and ion exchange (see section 
on chemical purification) and was finally obtained in 5 ml of 5 M hydro- 
chloric acid. This solution was used as a stock for all irradiations. Aliquots 
taken from it and examined by a-pulse analysis showed that less than 0.1% 
of the a-particles emitted had energies differing from the a-particle energy 
of Th™. 


(b) Itradiations 

Fifty-microliter aliquots (approximately 0.2 ug Th”°) of the stock solution 
were pipetted into small silica ampoules. Some of these were sealed directly but 
those to be irradiated under cadmium, and three others, were taken to dry- 
ness and were then sealed except for a capillary vent. Cobalt monitors, approxi- 
mately 10 mg in weight, were made from 0.005 in. cobalt wire and weighed to 
within +0.02 mg. These were wrapped helically around the sample ampoules 
which were then packed into super-pure iron irradiation capsules some of 
which were tightly lined with 0.025 in. cadmium. Reactor grade powdered 
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graphite was packed tightly into any remaining space and the capsules were 
then sealed. 

Except for one group in which three irradiations were carried out, the samples 
were irradiated in pairs one after the other (one cadmium shielded, the other 
unshielded) for 5400+2 seconds each in the same position in an empty fuel 
rod position of the NRX reactor during periods of constant pile power (con- 
stant to better than 40.5%). 


(c) Chemical Purification 

Cation and anion exchange chromatography were used to purify the 
irradiated thorium. 

Each cation column consisted of 1.00 g (dry weight) of 200-400 mesh 
Dowex 50 (hydrogen form) contained in a tube 10 cm long and 6 mm in dia- 
meter sealed onto a reservoir 10 cm long and 3 cm in diameter. Both column 
and reservoir were jacketed at 60°. Before loading, the columns were prepared 
by washing with water. 

Each anion column contained 0.25 g of 200-400 mesh Dowex 1 (chloride 
form) in a tube 6 cm long and 6-mm in diameter sealed onto a reservoir 5 cm 
long and 1.5 cm in diameter. The anion columns, which were operated at 
room temperature, were washed with 5 ml of 10 7 hydrochloric acid before 
loading. 

The irradiated thorium sample was dissolved in 5 ml of 4 M nitric acid and 
the solution was warmed to 80° before it was fed to a cation exchange column 
without application of pressure. When the feed had passed, 25 ml of hot 6 
nitric acid was poured into the reservoir and allowed to elute the column over 
a period of 1.5 hours. The effluent, which contained rare-earth impurities, 
was rejected. The resin was transferred to a conical flask by slurrying with 
20 ml of hot 0.2 7 ammonium oxalate solution and the thorium was extracted 
from it by boiling for 5 minutes. The resin was then filtered off and rejected. 

The filtered oxalate solution was transferred to a centrifuge tube and 6 
drops of a ferric nitrate solution (10 mg iron per ml) were added. The iron, 
which carries the thorium, was precipitated by adding excess ammonia and 
heating by immersion of the tube for 5 minutes in boiling water. The precipi- 
tate was collected by centrifuging, dissolved in concentrated nitric acid, and 
the solution was made alkaline once more. This precipitation—solution cycle 
was then repeated twice. The precipitate was washed with very dilute am- 
monia and then with water to remove any remaining ammonium nitrate. 

The precipitate was next dissolved in 2 ml of 10 M hydrochloric acid and 
fed to an anion exchange column at the rate of 0.2 ml/min to remove the iron 
carrier. After the feed had passed, the column was eluted with 1 ml of 10 AZ 
hydrochloric acid and the total eluate was collected in a small beaker and 
evaporated to dryness. Two milliliters of concentrated nitric acid was added 
and the solution was again evaporated to dryness. Finally the residue, which 
was invisible to the naked eye, was dissolved in 1 ml of 3 M nitric acid and the 
resulting solution was used for making sources for counting. The total time 
taken for each purification was about three hours. 
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(d) Absolute Disintegration-Rate Determinations 

Preliminary experiments were carried out to determine whether there were 
losses in counting Th?*° @-particles or Th®*! B-particles in a 4% counter owing 
to absorption of particles in the source mount. 

In the first experiment a series of gold-coated VYNS films of known thick- 
ness were prepared and an aliquot of the stock Th*° solution was evaporated 
to dryness on a film 44 wg/cm? thick. This was a-counted in a 4m counter and 
then another film was made to adhere to the back of the source before counting 
again. Additional films were then added, with counting after each addition, 
until a final total thickness of 209 wg/cm? had been built up. Over this range, 
44 to 209 ug/cm?, no change in a-disintegration rate was observed which was 
outside the statistical error of each count, i.e., +0.3%. 

In a second and similar experiment in which the Th*! 8-particles from a 
neutron-irradiated sample of Th?° were counted, no change in §-disintegration 
rate outside the statistical error of each count (+0.2%) was detected over 
a source mount thickness which was varied from 80 to 264 ug/cm?. Since for 
cross section measurements film thicknesses about 50 wg/cm? were used it was 
concluded that loss of 6-particles or a-particles due to absorption in these films 
was so small that it could be neglected. 

When a batch of irradiated Th?*® had been purified, three or four aliquots 
(sufficient to give between 300 and 1000 6-particle counts/second each but 
whose volume did not need to be known) were mounted on thin films and 
each was counted for 8-particles in a 4x counter. The progress of 8-decay was 
then followed for at least seventeen days by which time the Th?*!' content had 
decayed to a negligibly small amount and the (apparent) residual 8-particle 
count was only a small fraction of the background count. The (apparent) 
residual 8-particle count, which was due to the a-activity of the Th in the 
sample, was then subtracted from the B-activities determined previously to 
follow the decay of the Th**! component of the source. These corrected results 
showed no departure from exponential decay when plotted graphically and 
confirmed the purity of the product. From the corrected initial 8-particle 
count, the absolute 8-disintegration rate which the Th”! in the source would 
have had at the time of the end of the irradiation was calculated using 25.52 
hours (this investigation) for the half-life of Th**'. 

After B-particle counting had been completed each source was a-particle 
counted in the same counter at a lower applied potential which was in the 
middle of the a-particle counting plateau. The results for each film were then 
expressed as the ratio Rg. where 

absolute 6-disintegration rate of the Th?! in the source 
at the end of the irradiation 


(1) Resa = ——— en = —— 
absolute a-disintegration rate of the Th?*® in the source 





Standard deviations for the three or four separate determinations of Rg,« 
for any one irradiation were always less than 2%. The results for unshielded 
samples are distinguished by the subscript 1, i.e. (Rs,a)1, and those for shielded 
samples by the subscript 2. 
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(e) Measurements with Cobalt Monitors 

The determination of (mv), and (nv)... from the observed activity of 
cobalt monitors has already been described by Cabell, Eastwood, and Campion 
(1958). (mvo), is the “‘thermal’’ flux where n is the density of all neutrons 
below 0.5 ev (the cutoff energy for 0.025 in. cadmium), vo is the standard 
velocity, namely 2200 m/sec, and the subscript ¢ denotes ‘‘thermal” with this 
spectrum definition. The quantity (mv)... is the epicadmium flux per unit 
interval of In E. Throughout this paper a value of 5.28 years has been used 
for the half-life of Co® and the thermal capture cross section and resonance 
capture integral of Co have been taken as 36.5+1.0 barns (an average of 
the absorption and activation results) and 48.6 barns respectively (Hughes 
and Harvey 1955). 

Although (nv), and (nv)... were constant during each irradiation, (vo); 
varied between 6.83 and 6.99 X10" and (nv) ,/(nv)-.- between 29.9 and 38.3 
between separate irradiations. 


(f) Half-life of Th! 

Results from the decay for four half-lives of 18 separate sources of Th?*! 
prepared from five different irradiations and purifications were subjected to 
analysis by the method of least squares. The half-lives obtained are listed 
in Table I. The results gave an unweighted average of 25.52+0.01 hours where 
the error quoted is the probable error of the mean. A similar treatment of five 
sources for decay for six half-lives gave 11,2 = 25.59+0.05 hours. These 
results agree well with those obtained recently by Jaffey, Jerome, and Warshaw 


TABLE I 
DETERMINATION OF THE HALF-LIFE OF Th?! 





Irradiation T1/2 (hours) 
1 25.57 

25.44 

25.51 


25.43 
25.51 
25.43 
25.52 


to 


3 25.56 
25.47 
25.58 
25.51 


4 25.59 
25.42 
25.58 


5 25.66 
25.45 
25.56 
25.51 


Average 25.52+0.01 hours, 
where the error is the 
probable error of the mean 
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(1951) and by Knight and Macklin (1949). The former obtained the figure 
25.64+0.1 hours as an average of two determinations and the latter 25.514 
0.23 hours as a weighted average of five results. 


RESULTS 
We define an effective capture cross section for pile neutrons, oe, by the 
expression: 
‘ . Ati 
(2) ett = (Ra.a)iX1/ (nv) (1—e"), 


where (Rg.2)1 has been defined by equation (1) and ¢; hours is the time of the 
irradiation. \ and ), are the decay constants of Th?*! (0.02716 hours!) and 
Th*° respectively. Unfortunately A; is not known with a precision com- 
parable with that of the other factors in equation (2). To avoid introducing a 
large uncertainty at an early stage the results were expressed as the ratio 
Ger/A1 and are listed as such in column 4 of Table II. Since the groups of 
experiments were conducted over a period of time during which several fuel 
rod changes had been made it is not to be expected that oe/A, will be con- 











stant. 
TABLE II 
THE EFFECTIVE AND THERMAL CAPTURE CROSS SECTIONS OF Th?9 
Physical state 
of sample (ett,ca/A1) 
Group Experiment irradiated = (oets/A1) X10"? X 10 (@2200/A1) X 10" 

A 1 Solid 1.079+0.025 

2 Solid 1.871+0.014 7.93+0.28 
B 3 Solid 1.086+0.006 

4 Liquid 1.878+0.011 7.92+0.12 

5 Solid 1.836+0.012 7.50+0.13 
i 6 Liquid 0.991+0.017 

7 Solid 1.758+0.012 7.67+40.25 
D 8 Solid 0.985+0.007 

9 Solid 1.772+0.009 7.8740.12 


Weighted average 7.78+0.09 


Values of oete,ca/A1 are also listed in Table II]. The quantity oeg.ca, “the 
effective cross section under cadmium,”’ is defined by the relationship: 


(3) Sorr.ce = (Rae) XO Ga) (1). 


By the definition used for calculating (mv), from the monitor activities it 
will be seen that a2 m/sec for an absorber that, like cobalt, follows the 
1/v law, is simply ceg—err.ca- The values of (o¢—oete.ca)/A1 are therefore 
listed in column 6 of Table II as (e229) m/sec) /A1. 

The errors in oe/Ay and Gett.ca/A1 Which are given in Table II arise from 
two sources: an uncertainty in (mv), which was always 0.4% and an un- 
certainty in Rg.« for which the standard deviation of the three or four deter- 
minations has been taken (this varied between 0.3% and 2.0%). When an 
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average was obtained for o2290/A; each result was weighted proportionally to 
the reciprocal of the square of its error. 

For calculating an absolute value of og a value for the half-life of Th?*° 
must be assumed. The most precise value available, 8.23+0.25X10' years, 
was obtained by Curie (1930) but a more recent determination by Hyde 
(1949) has given 8.0+0.3 X 104 years. For want of better precision a value of 
8.0010! years has been assumed and a2 is expressed relative to this. 


Thus 
G20 = (7.78+0.09 X 10—"') K (2.747 K 10-) K 1074 = 21.4+0.3 barns. 


The resonance capture integral is related to the cadmium ratio (Req) by 
the expression: 


e say (MVo)e 2200 __ 
(4) Joe si eal (nv) e-c (Rea—1) ° 


Table III lists the values of Req and 
( o.dE F) /o2200 
0.5 


TABLE III 
THE RESONANCE CAPTURE INTEGRAL OF Th?3° 


obtained. 


( | i B) fom 
e 9.5 lees 





Group Experiments Rea 
A 1 and 2 1.735+0.041 42.942.4 
B 3 and 4 1.723+0.009 41.0+0.7 
3and 5 1.692+0.009 43.2+0.7 
Cc 6 and7 1.774+0.034 49.5+2.3 
D 8and9 1.799+0.013 47.6+0.9 


Since Rea = (Rs.a)1/(Rs.a)2 for equal irradiation times, errors in Req arise 
from both (Rg.a)1 and (Rg,2)2. Errors in (Rea— 1) were even larger and varied 
between 1.3% and 5.5%. The value of (nv)... was known with a precision of 
+1%. Using the weighted average, 


J o,dE/E = (43.8+1.4) X(21.440.3) = 937432 barns. 
0.5 
CONCLUSIONS 

The thermal capture cross section obtained in this work, 21.4 barns, is 
considerably lower than any figure reported previously. It is clear that the 
results of Jaffey and Hyde (1952), i.e., 61 barns and 28 barns, are much too 
high since these authors did not take into account any contributions due to 
epithermal neutrons and, with such a large resonance integral, their effect 
cannot be ignored. It is not possible to say why the new figure differs from the 
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values quoted by Hughes and Harvey (1955) since no experimental data are 
available for their results. 

Both the thermal capture cross section and the resonance absorption in- 
tegral obtained in this work are subject to correction when a more precise 
value has been obtained for the half-life of Th?°. With such a large resonance 
absorption integral it is also possible that a large resonance absorption peak 
is situated near the thermal region. If this is so, the tail of this peak may 
extend below 0.5 ev, making the assumption that Th**° behaves as a 1/v 
absorber in the subcadmium region invalid; a correction, probably not more 
than a few per cent, would then have to be applied. 
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NUCLEAR ORIENTATION EXPERIMENTS WITH Pr'”® AND 
Yb'> NUCLEI 


J. M. Danrets, J. L. G. LAMARCHE,? AND M. A. R. LE BLANC 


ABSTRACT 


Pr'? has been oriented as an impurity in cerium magnesium nitrate, and the 
anisotropy of the 1.57 Mev y-ray was measured as a function of temperature in 
the range 0.003° K to 1.0° K, both in zero external field and in a field of 320 gauss 
parallel to the trigonal crystal axis. The results are not consistent with what has 
been published about the spin Hamiltonian of Pr*** in this lattice. The magnetic 
moment of Pr! is very small, and it was not possible to see any saturation of 
the anisotropy at low temperatures. Values for the magnetic moment of Pr! 
were found for two different assumed decay schemes, these are 0.11 nuclear 


2 e 
magnetons for the spin assignments 2 “,2 +40, and 0.15 nuclear magnetons for 


the spin assignments 2 4,2 +50. Similar experiments were carried out on Yb". 
Measurements were made in a variety of external magnetic fields up to 700 
oersted, and at temperatures as low as 0.003° K. No anisotropy was observed 
for the 396 kev y-ray, nor for the 282 kev y-ray. The most likely explanation 
for this result is that the lifetime of the y-emitting state is about 107! seconds. 
The relation between the results of this experiment and others on the same 
isotope (Phil. Mag. 2, 1079 (1952) ) is discussed. 


1. INTRODUCTION 

The magnetic hyperfine structure method of orienting nuclei has been estab- 
lished by Bleaney ef al. (1954) and has already proved useful in establishing 
details of nuclear decay schemes, and in evaluating nuclear magnetic moments 
(e.g. Grace and Halban 1952). As part of a program to exploit the potentialities 
of such experiments, we have oriented Pr! by this method, and have observed 
and measured an anisotropy in the emission of the 1.57 Mev y-ray which 
occurs in the decay of this nucleus. A value of the magnetic moment of the 
Pr'® nucleus has been deduced from these measurements. 

We have also performed experiments to orient Yb'® as an impurity in 
cerium magnesium nitrate. When these investigations were started, it was 
discovered that a group at Oxford were preparing to align this isotope. At 
that time, very little was known about the effect of interactions on the angular 
distribution of gamma rays from an assembly of aligned nuclei, except that 
such effects could be very large (Grace et al. 1954). In addition, only one 
isotope, Co®, had been aligned in more than one crystalline environment 
(Bleaney et al. 1954; Ambler et al. 1953). Since the Oxford group planned to 
orient Yb! in the ethylsulphate lattice, and we planned to align it in the 
double nitrate lattice, it was thought that these investigations might comple- 
ment each other and yield information about the effects of interactions. The 
results of the investigations at Oxford have already been published (Grace 
et al. 1957), and are indeed very different from ours. 


1Manuscript received April 18, 1958. 4c ; 5 
Contribution from the Department of Physics, University of British Columbia, Van- 


couver, B.C. 
2Now at the University of Ottawa. 
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PART I—PRASEODYMIUM 142 
2. EXPERIMENTAL 

The experimental technique which we used is essentially that described in 
detail by Bleaney et al. (1954) (q¢.v.). A solution of a few milligrams of praseo- 
dymium nitrate in heavy water was irradiated in the B.E.P.O. pile at Harwell 
to an activity of some 8 millicuries. It was then brought to Vancouver by 
C.P.A.’s transpolar flight, and the active solution was poured into a saturated 
solution of about 20 g of cerium magnesium nitrate. Single crystals were then 
grown from this solution in a vacuum desiccator. Despite the short half-life of 
Pr'” (19 hours), it was possible within 2 days to obtain some 5 or 6 g of 
large, well-formed single crystals with a counting rate of about 200/second in 
the counting apparatus described below. These crystals, which are plates of 
trigonal symmetry, were mounted with parallel crystallographic orientation 
in a demagnetization cryostat, and cooled by adiabatic demagnetization. 
Temperatures of 3X10-*°K can be, and were, attained by the adiabatic 
demagnetization of this salt. After demagnetization, the magnet was removed, 
and two scintillation counters were placed in position round the specimen. 
During the period when the crystals warmed up from the low temperature 
attained on demagnetization to the helium bath temperature of 1.25° K, 
from 10 to 20 minutes according to the conditions prevailing during an experi- 
ment, the counting rate was measured in the two counters, and the temperature 
of the crystals was measured using a mutual inductance and a ballistic galvano- 
meter, as described by Bleaney eft al. (1954). The counters were thallium- 
activated sodium iodide crystals 1} in. diam.X1 in. long, pre-packaged by 
the Harshaw Chemical Co. In some of the experiments the cell containing the 
crystal was mounted directly onto the face of the photomultiplier tube (R.C.A. 
6342), in others it was coupled to this tube by a lucite light pipe 6 in. X 1} in. 
No significant loss of resolution was noticed with this latter arrangement. 
One scintillation crystal was placed on the trigonal axis of the radioactive 
crystals, about 10 cm from the source, and the other was placed at a similar 
distance in a direction at right angles to this axis. The pulses from each photo- 
multiplier were first amplified, then fed into a discriminator which passed 
only those pulses of energy greater than 1 Mev, and finally into a scaler. There 
were two such identical channels, and the pulses from each photomultiplier 
were counted on a separate scaler. After the active crystals had warmed to 
1.25° K, the counting rate was observed for a while to test the stability of the 
counters, and for normalization purposes. In some of the experiments, an 
external magnetic field of a few hundred gauss was applied to the radioactive 
crystals. This was produced by two water-cooled air-core coils in a quasi- 
Helmholtz arrangement. When these coils were used, the scintillation crystals 
were coupled to the photomultipliers by the light pipes, and the photo- 
multipliers were shielded with one layer of mu metal and two layers of annealed 
soft iron. A bucking coil was wound round the soft iron shield of one of the 


photomultipliers. 
We can define an anisotropy parameter e by the relation e = (J, —J,)/J,, 
where J, and J, are the normalized counting rates in directions parallel and 
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perpendicular respectively to the trigonal axis of the crystals. Values of « 
for the same temperature, but from different demagnetizations, were averaged 
to improve statistics. 

The absolute temperature of the radioactive crystals was obtained from 
magnetic susceptibility measurements, using the relationship given by Daniels 
and Robinson (1953). The Curie-Weiss A for our specimen was determined 
by comparing the measured value of the magnetic temperature 7* after 
demagnetization with the value calculated from the initial magnetic field 
and the initial temperature. Its value was 2.6 millidegrees, which is a reason- 
able value according to Daniels and Robinson (1953). 

The results of the measurements are shown in Figs. 1 and 2, in which ¢ 
is plotted as a function of 1/7; Fig. 1 is for the case where no external magnetic 


L ! 
€ | 








- 0-02 ; 


Fic. 1. Anisotropy of the 1.57 Mev y-ray of Pr. No external magnetic field. 


field is applied, Fig. 2 is for the case where an external magnetic field of 320 
oersted is applied parallel to the trigonal axis of the crystals. In each figure 
the solid curve is the curve of the form ¢€ = a/7? which best fits the data. For 
zero external magnetic field a = 4.95X10~7 °K~-*, for an external field of 320 
oersted a = 18.7X 1077 °K. 


3. DISCUSSION 


The double nitrates M2Mg3(NO3)12.24H2O, where M is one of La, Ce, or 
Pr, are isomorphous and form mixed crystals in all proportions. The para- 
magnetic resonance spectrum of the stable isotope, Pr'', as an impurity in 
lanthanum magnesium nitrate, has been observed by Cooke and Duffus 
(1955), who assigned to this ion the spin Hamiltonian: 


(1) KH = ¢\BH,S,+ AS,1,+42S:+A4,S,, 
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Fic. 2. Anisotropy of the 1.57 Mev y-ray of Pr‘. An external magnetic field of 320 gauss 
applied along the trigonal axis. 


where 
gy, = 1.55, 
A = 0077 cm", 
A = Average VaZ+A,? ~ 0.04 cm=, 
S = 4, 
\ 5 2 


and the z-axis is the trigonal axis of the crystals. The Pr*+*+* ion has an even 
number of electrons, and hence the energy levels need not be doublets (Kramers’ 
theorem). However the symmetry of the double nitrate crystal is such that the 
ground state is an accidental doublet. Strains and other distortions can cause 
a splitting of this doublet, and it has been suggested that this is the origin 
of the A-term. There is no unique value of A, a different value being appropriate 
for each ion, but an ‘‘average’’ value of A can be assigned. A is probably 
structure-sensitive. The most general spin Hamiltonian with trigonal sym- 
metry has terms g,6(/7,S,+H,S,) and B(S,P,+5S,l/,) in addition to those 
in equation (1). To the accuracy of published measurements, both g, and 
B are zero. In the case where S = 3, g, = B = 0, and A = 0, magnetic 
interactions between the ions do not affect the angular distribution of y-rays 
(Daniels 1957). In this case also, an external magnetic field applied along the 
z-axis has no effect on the angular distribution (Daniels 1957). For ‘‘high”’ 
temperatures (in this experiment, 7 > 0.005° K) the y-ray anisotropy 
e « 1/7*. The angular distribution of y-rays can be represented by an operator 
(1/4r+I’), and the angular distribution at any temperature is given by 
1 Spur (1’#*) 


1 
9 — —---+-- aoa. Pa ss 
(2) I(6, $) 4ar 2k°T” = Spur (1) 


(see Daniels (1957), in which some of the symmetry properties of I’ are also 
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given). Substitution into (2) of the Hamiltonian (1) shows that the A-term 
does not affect the 1/7? term in the angular distribution. It is thus difficult to 
understand the reason for the difference between the values found for ¢ with 
and without an external magnetic field. Many misgivings have been voiced 
that the spin Hamiltonian obtained from resonance experiments in a large 
magnetic field does not adequately express the state of an ion in zero 
magnetic field. This may be particularly so for praseodymium, which has no 
Kramers degeneracy. Under the circumstances we choose to use the values of 
€ in an external magnetic field to determine the magnetic moment of Pr'*. 
This is because: 

(a) The spin Hamiltonian (1) was determined in a large external magnetic 
field, and is probably more applicable to the case where there is an external 
field than to the case where there is none. 

(6) The measurements in this case were made at higher temperatures. 
The temperatures are more reliable, and the series expansion in powers of 
1/T (2) is more likely to be applicable. 

(c) The warming up times were longer and hence temperature inhomo- 
geneities are not so likely. 

(d) If the assumed spin Hamiltonian is inadequate, and for some reason 
magnetic interactions play a part in determining the anisotropy, an external 
field of 320 oersted should be sufficient to nullify these interaction effects, 
since the internal field in this salt is about 30 oersted. 

The possibility was considered that bremsstrahlung from the 2.166 Mev 
B-ray could be counted as y-rays and reduce the observed anisotropy. Calcu- 
lations based on formulae in Siegbahn (1955) showed that this effect is neg- 
ligible. The threshold of the discriminator was set at 1 Mev to get a good 
counting rate. In this case one must consider the reduction of anisotropy due 
to scattered y-rays. Such a calculation was made by Bleaney ef al. (1954) 
for the case of Co® in an apparatus which is essentially similar to ours. It 
was found there that this correction is quite negligible for values of « as low 
as 0.04, and we have therefore ignored it. It is very unlikely that the lifetime 
of the y-emitting state is greater than about 10-™ seconds, hence loss of 
anisotropy due to precession in the crystalline fields can be discounted. 

The decay scheme of 19.2h Pr! has been studied by various workers, in- 
cluding Gideon et a/. (1949), Jensen et al. (1950), Bartholomew and Kinsey 
(1953), Polm et al. (1954), and Sterk et al. (1955). The decay scheme suggested 
by Polm et al. (1954) is shown in Fig. 3. The arguments for this decay scheme 
are, briefly, as follows. Nd! is an even-even nucleus; its ground state there- 
fore is 0+. The 2.166 Mev 8-ray has log ft = 7.8, and the first special shape 
appropriate to a beta transition with AJ = 2, yes; the spin and parity of 
Pr'” are thus 2—. The other beta transition has log ft = 7.1 and allowed 
shape; it is therefore suggested that it is first forbidden with AJ = 0 or +1, 
yes. Since the first excited states of nearly all even—even nuclei so far measured 
have spin and parity 2+, the excited state of Nd!” in this scheme is assigned 
spin and parity 2+, and the y-ray is assumed to be £2. 

Since the spin Hamiltonian (2) indicates that at 0° K all the nuclei are 
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oriented parallel to the trigonal crystal axis, it is possible in principle to 
distinguish between certain possible spin assignments in the decay scheme. 


142 
Pr 
Songer e- 
\ 4% 
'; - 
\ 
96% 
E2 
57 MV 
Nd'*? O+ 


Fic. 3. Decay scheme of Pr!®. 


For example, the limiting value of ¢€ at low temperatures is +1 if the 0.586 Mev 
8-v pair carries off no angular momentum; +1/3 if the 8-» pair carries off one 
unit of angular momentum, the other assignments being as in Fig. 3. However, 
we could not attain temperatures low enough to see any saturation of the 
y-ray anisotropy. The observed positive sign of ¢ indicates that the y-ray is 
quadrupole. 

Substitution of (1) into (2) yields the formulae, ¢ = 34?/16k?7? if the 
B-» carries off no angular momentum, and e = 3.4?/32k?T? if the 6-v carries 
off one unit of angular momentum. Since the field at the praseodymium 
nucleus is due to the 4f electron shell, independent of the isotope, this field 
Is H = Aggy ly41/mig, = Agel 142/p142. These values are known for the stable 
isotope Pr'', and it is a simple matter to substitute into these formulae to 
find wigs. The result is: wigg = 0.11 nuclear magnetons if the B-v carries off no 
angular momentum, or wi42 = 0.15 nuclear magnetons if the 6-y carries off 
one unit of angular momentum. 


PART II—YTTERBIUM 175 
4. EXPERIMENTAL 

The starting material for these investigations was Yb2O3, separated by the 
ion exchange method, and supplied to us by Dr. F. H. Spedding of the Ames 
Institute One milligram was irradiated in the NRX pile at Chalk River, 
to an activity of about 5 me. On receipt of the irradiated sample, the quartz 
capsule was crushed under about 20 cc of water contained in a platinum 
crucible. A special jig was used for this operation. About three drops of 
concentrated sulphuric acid were added, and the whole was evaporated to 
dryness, first on a sand bath and afterwards under a heat lamp. The residue 
was then dissolved in about 20 or 30 cc of water and added to about 10 g 
of cerium magnesium nitrate. Crystals were grown from this solution in a 
vacuum desiccator. 
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These double nitrates of the rare earths exist only for the rare earths from 
La to Gd, and ytterbium magnesium nitrate does not exist as such. However, 
some ytterbium does go into the crystal lattice. After a crystal was removed 
from the solution, it was washed with distilled water to remove any solution 
from its surface, and dried with filter paper, to preclude the possibility that 
the activity is due to a layer of dried off solution on the surface. The crystals 
were generally clear enough for the observer to see that they contained no 
large inclusions of solution. On one occasion, a sample of Yb2O; became con- 
taminated in the irradiation, and the solution from which the crystals were 
grown emitted a gamma ray of about 700 kev energy, yet this gamma ray 
was completely absent from the crystals grown from this solution. In this 
way we are sure that the ytterbium does go into the crystal lattice; it is most 
reasonable to assume that it goes into a rare-earth position, although there 
is no direct evidence for this. The specific activity of the crystals was about 
1/30 that of the solution. 

During the irradiation of natural ytterbium, three isotopes are formed, 
4.2day Yb", 32.4day Yb", and 6.7 day Lu'”. However, about 10 times as 
much Yb" as Yb'® or Lu'” is produced (Cork et al. 1950), and the 396 kev 
y-ray of Yb’ has a greater energy than any other from any of the three iso- 
topes. Yb!®* has a y-ray of 300 kev, but this is less intense than the 282 kev y-ray 
of Yb">. The technique used to orient the ytterbium was exactly the same as 
that used to orient praseodymium, described in Part I of this paper. Measure- 
ments were made on both the 396 kev y-ray and also on the 282 kev y-ray; 
for the former a bottom cut discriminator was used, set to pass the photo- 
peak of the 396 kev y-ray; for the latter a single channel kicksorter was used, 
set to pass only the photo-peak of this y-ray. These measurements were made 
with zero external magnetic field, with fields of 260, 380, 510, and 700 oersted 
parallel to the trigonal crystalline axis, and with fields of 58, 100, and 160 
oersted perpendicular to this axis. In no case was any significant anisotropy 
of y-radiation observed at any temperature. The results are shown in Fig. 4. 


5. DISCUSSION 

The principal features of the decay scheme of Yb" are shown in Fig. 5. 
This decay scheme has been studied by a variety of investigators (de Waard 
(1955), Akerlind et al. (1955), Marty (1955), Mize et al. (1955a, 6), Cork 
et al. (1956)) and has been analyzed theoretically by Chase and Wilets (1956). 
We shall therefore assume that the decay scheme is established correctly. 
From measurements of the anisotropy of y-rays from Yb" oriented in ytter- 
bium ethylsulphate, Grace et al. (1957) deduced a value for the magnetic 
moment of Yb™; it is 0.15+0.04 nuclear magnetons; the sign is not known. 
The 396 kev y-ray and the 282 kev y-ray are both M2 with a small admixture 
of El. The experiments so far carried out do not give consistent values for the 


mixing ratios of these y-rays (see Grace et al. 1957). 

No one has observed paramagnetic resonance in Yb*+** in the double nitrate 
lattice, hence the spin Hamiltonian is not known. We have made an estimate 
of the constants in the spin Hamiltonian by the following method. We assume 
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_ Fic. 4a. Anisotropy of the 396 kev y-ray of Yb"®. Upper graph—no external magnetic 
field. Middle graph—external magnetic fields applied along the trigonal axis. Lower graph— 
external magnetic fields applied perpendicular to the trigonal axis. 
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_ Fic. 4b. Anisotropy of the 282 kev y-ray of Yb". Upper graph—no external magnetic 
field. Middle graph—external magnetic fields applied along the trigonal axis. Lower graph— 
external magnetic fields applied perpendicular to the trigonal axis. 
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that the crystalline field acting on the ytterbium ion is the same as that 
acting on the cerium ion; the coefficients for the expansion of this field in a 


Yb’? 
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EIl+M2 
396 KV 


i 
Jot 





ia 


Fic. 5. Decay scheme of Yb!®. 


series of spherical harmonics are given by Judd (1955a). We assume that the 
formula given by Elliott and Stevens (1953a) for the radius of the 4f shell, 


r® « (Z—55)-"4, 


is applicable. Using these values of crystal field parameters and ionic radii, 
we can work out an energy matrix for the ytterbium ion using the method and 
formulae given by Stevens (1952), Elliott and Stevens (1952, 1953a, 6), and 
Judd (1955a, b). The results of these calculations are as follows: whereas the 
free Yb*** ion has J = 7/2 and an eight-fold degenerate ground state, in the 
double nitrate environment, this degenerate level is split into 4 doublets 
whose energies are —2.26, +39.4, +0.4, and —37.5 cm—!. The state with 
energy —2.26cm™ isthe doublet spanned by |J, = 3/2) and |J, = —3/2), 
which lies lowest in the ethylsulphate. The lowest energy level is the one at 
— 37.5 cm~!. This is spanned by the wave function 


0.794 |7/2) —0.424 |1/2) —0.435 | — 5/2) 
and 
—0.794|—7/2)—0.424 | —1/2)+0.435 |5/2). 


Assuming a spin Hamiltonian of the form 
g \GH.S.+2,8(H2S:+H,S,)+ AS 1,+B(Slr+S yl y), 


we find g, = 4.2, g, = 2.9. The extent of the hyperfine structure is given by AJ 
and BI. Assuming that the magnetic moment of Yb! is 0.15 nuclear magnetons, 
AI = 0.024° K and BI = 0.017° K. It is true that it is rather optimistic to 
expect such a calculation to give a result which is quantitatively correct, but 
unfortunately no better estimate of the properties of the ytterbium ion in this 
lattice exists. On the basis of these estimates, we should expect to observe 
nuclear orientation under the conditions of our experiments. 
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In view of the established fact that sizeable anisotropies of both the 396 kev 
y-ray and the 282 kev y-ray have been observed from oriented Yb" nuclei 
(Grace et al. 1957), the absence of an anisotropy in zero external fieid in our 
experiment can be explained by one or more of the following hypotheses: 

(i) The ytterbium ion does not enter the lattice in a rare-earth position. 

(ii) The spin Hamiltonian is isotropic, i.e., g; = g, and A = B. 

(iii) Magnetic interactions in the double nitrate crystal destroy the align- 
ment of the ytterbium nuclei. 

(iv) The y-emitting state is long-lived and precession of the excited Lu'® 
nucleus in the internal atomic fields destroys the alignment before the emission 
of the y-ray. 

As for hypothesis (ii), the application of an external magnetic field in any 
direction should cause an anisotropy of y-emission to appear (Gorter 1948; 
Rose 1949; Ambler et al. 1953). Since this does not occur, the possibility 
of the other hypotheses must also be considered. 

Hypothesis (iii) is certainly applicable in this case. The theory of the effect 
of interactions has been given. by Daniels (1957), and an order of magnitude 
calculation based on the formulae given there shows that the magnetic inter- 
ference can be much greater than the forces which tend to orient the nuclei. 
Naively, the internal field at a rare-earth site is about 30 oersted; this will 
produce a splitting of the ground state of about 0.01° K, whereas the forces 
which tend to produce spatial orientation are of the order of magnitude of 
(1—B)I, or about 0.007° K. However, the application of an external field 
larger than the internal field of 30 oersted should cause an anisotropy to appear, 
and this is not so. 

Hypothesis (iv) is quite sufficient to explain our results, and also to explain 
the difference between our results and those obtained at Oxford. In the spin 
Hamiltonian for the Ybt+++ ion in the ethylsulphate, g, = 0 = B. Under 
these circumstances, the angular distribution of y-rays is not affected by 
precession of the y-emitting nucleus. This is easily seen, for the angular dis- 
tribution of y-rays is given by an expression of the form 


u aW,\C lyr), 


where g, is a statistical distribution (e.g., Boltzmann) factor, ¥, is a wave- 
function of the ion containing the y-emitting nucleus, and [ is an operator 
which represents the angular distribution of y-rays (see e.g., Daniels 1957). 
Because of the axial symmetry of the ion and its environments in a representa- 
tion with S, diagonal only the diagonal elements of T appear in the result. 
If there is precession for a time ¢t before the emission of the y-ray, |y,) must 
be replaced by exp|—i.AHt/h]\¥,). This is equivalent to replacing I by 
exp[i #2/h]| T exp[—i #1/h]. Now the spin Hamiltonian ” for ytterbium ethyl- 
sulphate is diagonal in this representation and hence the diagonal elements 
of exp[t Ht/h] T exp[—1 Wt/h] are independent of ¢. Thus the angular distri- 
bution of y-radiation in ytterbium ethylsulphate is not affected by precession. 
This is not so for the more general spin Hamiltonian which presumably applies 
to the double nitrate. It is difficult to estimate the lifetime necessary to wipe 
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out all the anisotropy which should be there. If we take as a rough criterion 
that the lifetime should be about as long as it would take for the excited 
Lu'® nucleus to precess in the atomic magnetic field, the result is about 
10-!° seconds. Chase and Wilets (1956) have explained that for both y-rays, 
the El transition is strongly inhibited, and the lifetime should be about that 
appropriate to an M2 transition. Weisskopf’s formula for the lifetime of a 
400 kev M2 emitting state gives 10-* seconds; since the transition in Yb!” 
is a collective transition, not a single nucleon transition, a lifetime of 10—!° 
seconds is more likely. A question arises whether the spin Hamiltonian of the 
ytterbium ion should be used during the precession, since after the y-emission 
the electronic configuration should be that of lutecium, which is not magnetic. 
This depends principally on the lifetimes of the 4f states. These are probably 
about 10~* to 10-* seconds, typical lifetimes for optical transitions, and in the 
absence of further information the use of the ytterbium spin Hamiltonian is 
justifiable. Thus, it appears that hypothesis (iv) is adequate to explain our 
results. 
6. CONCLUSIONS 

(i) Praseodymium 142 

An anisotropy of y-emission from oriented Pr'* has been observed and 
measured. Although the anisotropy is not completely understandable in terms 
of the published spin Hamiltonian Pr*++*, it seemed plausible to assign a value 
to the magnetic moment of this nucleus. The value assigned is 0.11 nuclear 
magnetons if the cascade 6-ray carries away no angular momentum, and 0.15 
nuclear magnetons if the 6-ray carries away one unit of angular momentum. 


(ii) Vtterbium 175 

No anisotropy was observed for either y-ray of ytterbium 175 in cerium 
magnesium nitrate. Although the possibility must always be kept in mind 
that ytterbium does not replace cerium in cerium magnesium nitrate, it seems 
that the most likely reason for not observing any y-ray anisotropy is that the 
y-emitting state has a lifetime of about 10-!° seconds. The observation of 
anisotropies from this same nucleus in ytterbium ethylsulphate is quite con- 
sistent with this hypothesis, since the internal fields in this salt are such that 
precession would not affect the angular distribution of y-rays. The estimate of 
10-'° seconds for the lifetime of the y-emitting state is consistent with all the 
other evidence about this nucleus. 
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THE ELECTRICAL CONDUCTIVITY OF 
MANGANESE ARSENIDE AND ANTIMONIDE'! 


GASTON FISCHER AND W. B. PEARSON 


ABSTRACT 


New measurements of the electrical resistivity of MnAs and MnSb are reported 
as part of a program of relating the electronic and crystal structure of com- 
pounds to their electrical properties. Modern experimental methods have been 
used in an attempt to obtain better samples and to avoid the formation of cracks 
which greatly influence the measured resistivity. The electrical resistivity of 
MnAs in the paramagnetic range where it is fairly independent of temperature 
is about 5X1074 ohm cm. In the ferromagnetic range it is about 210-4 ohm 
cm at 0° C with a temperature coefficient of 5X 1073 per ° C. 


INTRODUCTION 


Manganese arsenide, antimonide, and bismuthide each have the hexagonal 
nickel arsenide type of structure. The substances are interesting for the 
transitions which they undergo (Table 1). As part of a program designed to 
gain more understanding of the relationship between electronic and crystal 
structure on the one hand and electrical properties on the other hand, we 


TABLE I* 
TRANSITIONS IN MNAs, MNSs, MNBr 


MnAs_ 25°-50° C With hysteresis. Ferromagnetic <> paramagnetic. First order transi- 
tion 
~130° C A second order transition appears to occur in this region 
MnSb_ 313°C Ferromagnetic < paramagnetic. Normal second order transition 
MnBi = 340°-360° C Ferromagneiic < paramagnetic. First order transition resulting from 


some Mn atoms changing their sites in the crystal structure 


*Cf. Guillaud 1951; Bates 1934; Willis and Rooksby 1954; Bacon and Street 1955; Roberts 1956. 


have made careful measurements of the electrical conductivity of MnAs and 
MnSb. Although this property has been measured previously (Guillaud 
1951; Bates 1934) the details of the measurements of MnAs, particularly, 
are open to some doubt owing to the micro- and macro-cracks which permeate 
samples after they have been cooled through the ferromagnetic transition at 
25°-50° C, and furthermore, an absolute value of the resistivity has not so 
far been obtained. It was thought that some improvement over the earlier 
measurements could be obtained (1) by never allowing the samples of MnAs 
to cool through the transition at ~40° C before the measurements were made 
(this should avoid the formation of cracks in the specimen which were thought 
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Contribution from the Division of Pure Physics, National Research Council, Ottawa, 
Canada. 


Issued as N.R.C. No. 4830. 


Can. J. Phys. Vol. 36 (1958) 


1010 





FISCHER AND PEARSON: ELECTRICAL CONDUCTIVITY 1011 


to arise at the discontinuous anisotropic change of unit cell dimensions occur- 
ring at the transition); (2) by using, in the preparation of the sample, some of 
the methods now employed in the preparation and handling of semiconducting 
compounds. Although zone purification does not appear feasible, slow cooling 
and solidification of the melt from one end should lead to the formation of 
perfect stoichiometric crystals at least at the start of the process. Further- 
more, by a suitable distribution of temperature gradients along the sealed 
silica tubes in which samples are prepared, it should be possible to prevent 
distillation of the volatile components from the surface of the sample. 

We were also particularly interested in making measurements of samples 
of MnBi because the work of Roberts (1956) using neutron diffraction has 
established the nature of the change at 350° C which, as in MnAs, is a first 
order transition. Unfortunately, however, we have not succeeded in preparing 
samples of MnBi of sufficient quality for conductivity measurements. 


EXPERIMENTAL METHODS 

The principles outlined above were followed during the preparation of the 
samples im vertical or horizontal gradient furnaces. Pure hydrogen-treated 
electrolytic manganese was mixed with the required proportion of 99.9% As 
or 99.995% Sb and sealed in an evacuated silica tube. The mixture was slowly 
heated up to the melting point and then held some 25° C above this tem- 
perature for 1 hour. After a vigorous shaking it was solidified slowly and 
finally cooled either to 50° C or to room temperature over a period of 2 days. 

Microscopical and X-ray examination showed the MnAs samples to be 
homogeneous but MnSb samples contained about 4 volume per cent of other 
phases, even after they were annealed at 800° C for several days. The samples 
were polycrystalline and strongly ferromagnetic at room temperature. 

The electrical resistivity of the samples was measured by a potential method 
(Dauphinee and Mooser 1955). During measurement the specimens were 
kept under one atmosphere pressure of helium to prevent loss of As or Sb by 


evaporation. 


ELECTRICAL RESISTIVITY 

(a) MnAs 

The temperature dependence of the electrical resistivity of two samples 
(A and B) of MnAs which had never been cooled through the transition at 
~40° C is shown in Fig. 1 ((a) and (6)). However, even when the polycrystal- 
line samples were heated and cooled a number of times without passing 
through the transition temperature, it appeared that microcracks were 
formed since the specimens became softer and there were small changes in 
resistivity. This behavior probably resulted from the strongly anisotropic 
thermal expansion of the material (cf. Willis and Rooksby 1954) and its 
inability to be plastically deformed. On the other hand one sample (C) pre- 
pared in the vertical furnace remained very strong and showed no tendency 
to crumble after cooling to room temperature. This sample also was the only 
one to show a marked change of slope in the resistivity temperature curve 
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Fic. 1. Electrical resistivity of MnAs as a function of temperature. (a) Measurements 
obtained during the first heating cycle of a sample (A) which was never cooled below 40° C. 
The sample was prepared horizontally and the last portion to crystallize was measured. (6) 
Measurements obtained after several heating cycles: the resistivity starts to rise near the 
transition temperature. The sample (B) was taken from the same ingot used in (a) but the 
first portion to crystallize was measured in this case. 


(Fig. 1 (c)) in the region of 130° C where the presence of a further ‘‘transition’ 
had been suggested by lattice spacing variation and other evidence (cf. 
Guillaud 1951; Willis and Rooksby 1954). 

The electrical resistivity of samples A and B (~6.6 X 10-4 ohm cm at 130° C) 
was not much larger than that of sample C (~4.8X10-! ohm cm at 130° C) 
at temperatures where the samples could be compared; however, when samples 
A and B were eventually cooled through the transition at ~40° C their resis- 
tivity increased greatly, being then about a hundred times larger than that 
of specimen C. Such behavior is entirely consistent with the nature of the 
transition which is described in the following paper (Basinski and Pearson), 
since the effect on the resistivity of the fragmentation which occurs will 
depend on the relative grain sizes and orientations in the samples. 

In the paramagnetic range above 50° C the reasonable agreement of the 
resistivity of samples which had not been cooled through the transition with 
that of the sample which did not crumble on passing through the transition 
establishes the electrical resistivity of MnAs as about 5X10-* ohm cm and 
as being relatively independent of temperature. In the ferromagnetic range 
we believe that sample C, which remained mechanically strong and did not 
crumble, gives the true electrical resistivity of polycrystalline MnAs. Between 
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Fic. l(c). Electrical resistivity of MnAs as a function of temperature. This sample (C) 
was prepared vertically and it remained strong and did not crumble on cooling through the 
transition at about 40° C. 


—120° C and 20°C its resistivity varied linearly with temperature with a 
coefficient of +5X10-* per °C, the resistivity at 0° C being 1.710-4 ohm 
cm. 

Finally we must add that at the time when these measurements were made, 
we were not aware of the mechanism of the transformation at 40° C. It was 
rather that the peculiar variability of the electrical resistivity which we en- 
countered led us to make the detailed study of the transformation mechanism 
reported in the following paper (Basinski and Pearson). 

(6b) MnSb 

The variation of the electrical resistivity of two samples of MnSb during 
various heating cycles is shown in Fig. 2, and similar behavior has been ob- 
served in three other samples examined. In particular the behavior in the 
region of the Curie point is very variable. Nevertheless the measurements in 
the ferromagnetic range establish the temperature coefficient of resistivity 
(assumed to be constant) as about 3X 10-* per °C and indicate a resistivity of 
the order of 2X 10-* ohm cm for MnSb at room temperature. 


DISCUSSION 


The value of the electrical resistivity of MnAs in the paramagnetic range 
suggests that the substance is metallic and we therefore think that there is 
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Fic. 2 Electrical resistivity of MnSb as a function of temperature. (a) Measurements 
made on a sample during the first heating and cooling cycle. The heating curve can be con- 
sidered typical of MnSb. (b) Measurements made on another sample of MnSb during successive 
heating and cooling cycles, showing typical variations of resistivity encountered. 


direct interaction between the d electrons of the manganese atoms since, in 
the absence of such interactions, we would have expected MnAs to show 
semiconducting properties. However, the discussion of Pearson (1957) shows 
that the value of the axial ratio of the unit cell of MnAs, c/a = 1.53, is close 
to a critical region of values where a change from metallic to semiconducting 
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properties occurs in substances with the nickel arsenide structure. The inter- 
actions between the manganese atoms would not therefore be expected to be 
well developed, and this is supported by the electrical resistivity of MnAs 
of about 5X10-! ohm cm at 50° C, which is high for a good metal. The un- 
paired electrons on the manganese atoms should thus be regarded as mainly 
localized and a large contribution to the electrical resistivity from spin-dis- 
order scattering may occur in the paramagnetic range (cf. Mott and Stevens 
1957; Coles 1958). Reduced pr) curves for MnAs and MnSb are compared 
in Fig. 3 with those of other substances in which spin-disorder scattering 
is held to make a large contribution to the electrical resistivity. The dis- 
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Fic. 3. Reduced resistivity versus reduced temperature for several ferromagnetic and 
antiferromagnetic substances, @, is the transformation temperature as determined by electrical 
measurements; it may differ appreciably from the results of magnetic or structural observa- 
tions as in the cases of Dy and Mn. The data pertaining to the graph are: 





6-in°K pein Q“ em Reference 

Fe 1033 — Mott and Stevens (1957) 
Mn 65 140X10-* White and Woods (1957) 
Gd 289 -— Gerritsen (1956) 

Dy 170 120X10-8 Gerritsen (1956) 

AuMn 500 75 xX 10-8 Gerritsen (1956) 

Au;Mn 150 4510-8 = Gerritsen (1956) 

AsMn 321 505 X 107% Present work 

MnSb 580 295 X 1075 Present work 


continuous increase in resistivity at the transition in MnAs which accom- 
panies a similar increase in axial ratio indicates, as expected, a decrease in the 
charge carrier concentration. Because of the “‘critical’’ value of the axial ratio 
of this compound the carrier concentration should be sensitive to small changes 
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in the unit cell dimensions. In these circumstances the anomalous temperature 
dependence of the electrical resistivity above the transition at 40° C might 
also arise in part from an abnormal variation of charge carrier concentration 
and mobility with temperature. 
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THE NON-MARTENSITIC DIFFUSIONLESS TRANSITION 
IN MANGANESE ARSENIDE AT ABOUT 40° C'! 


Z. S. BASINSKI AND W. B. PEARSON 


ABSTRACT 


We have investigated the nature of the transition at about 40°C in man- 
ganese arsenide. It is already well known that the transition involves a dis- 
continuous change of the length of the a axis of the hexagonal unit cell, and that 
there is no change of crystal structure. The transition is accompanied by a loss 
of ferromagnetism of MnAs on proceeding to higher temperature. We find 
that the transformation takes place by a novel diffusionless non-martensitic 
mechanism. It appears that plastic deformation of MnAs is impossible and so, 
instead of the change following the martensitic pattern, the strain energy of 
the transformation is minimized by a fragmentation of the crystal grains into 
pencil-like subgrains whose direction lies parallel to the c axis of the crystals. 


INTRODUCTION 

Manganese arsenide has an interesting transition which occurs with hyster- 
esis just above room temperature (between 25° and 50°C). Both high and 
low temperature phases have the nickel arsenide type of structure; however, 
the discontinuous changes of lattice spacing along the a axis (Guillaud 1951; 
Willis and Rooksby 1954) and of the electrical resistivity (Guillaud 1951; 
Fischer and Pearson 1958) indicate a first order transition. The low tempera- 
ture phase is ferromagnetic, but the ferromagnetism disappears at the transi- 
tion point (Bates 1934; Guillaud 1951). The positive coefficient of susceptibility 
which is found between the transformation point and a temperature of about 
130° C has been interpreted by Guillaud as evidence of antiferromagnetism in 
this temperature interval, particularly as a ferromagnetic—antiferromagnetic 
change would be of the first order. However, a recent neutron diffraction study 
by Bacon and Street (1955) revealed no evidence of antiferromagnetism in 
the high temperature phase. 

Let us note that MnBi also undergoes a first order transition with rather 
similar characteristics at about 350° C. Roberts (1956) in a neutron diffraction 
study found that the Mn atoms were definitely in a paramagnetic state above 
the transition temperature. At the transition a movement of some 10% of the 
Mn atoms takes place randomly from the 0, 0, 0; 3, 3, 3 sites in the nickel 
arsenide structure stable at low temperatures to the trigonal bipyramidal 
holes at 3, 3, 3; 3, 4, 3. Himmel and Jack (1956) found that the high-temperature 
MnBi phase could be retained by quenching and that it was also ferromagnetic 
at room temperature, the Curie point being at about 200° C. Although there 
is an apparent similarity of the transitions in MnAs and MnBi it appears 
improbable that the mechanism of the change would be the same in each 
case, since we would not expect a normal diffusion-controlled reaction to 
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occur in MnAs at temperatures as low as 40° C. Indeed, in MnBi we know 
that a diffusion-controlled reaction can occur at an appreciable rate only at 
temperatures above 170° C, because not until this temperature is reached on 
heating up from room temperature, does the quenched metastable MnBi 
phase revert to the stable low-temperature phase. It thus appeared likely 
to us that the transition of MnAs would be a diffusionless type. All known 
diffusionless transformations involving an appreciable change of the unit cell 
dimensions are of the martensitic type, i.e., the strain energy is minimized 
during the change by fitting the parent and daughter phases coherently along 
a “habit”? plane which either remains macroscopically unrotated and un- 
disturbed (Wechsler, Lieberman, and Read 1953), or undergoes only a small 
dilatation (Bowles and MacKenzie 1954) during the transformation. To 
accomplish this, the strain responsible for the change of the lattice is generally 
homogeneous only over very small regions; macroscopically the adjacent seg- 
ments are either separated by fine slip, or undergo different deformations 
with the result that the daughter phase consists of fine lamellae of twins. As, 
however, MnAs is extremely brittle, showing no traces of plastic deformation, 
and as there are no known mechanical twins occurring in the material, it 
might not be expected to be able to undergo a normal martensitic type of 
diffusionJess transition. 

We have therefore decided to study the transformation in more detail to 
try to find out the mechanism of the change. 


EXPERIMENTAL INVESTIGATION 

In the first place we have attempted to retain the high-temperature para- 
magnetic phase by rapid quenching, but the ferromagnetic phase was always 
obtained when small samples were quenched into iced water, acetone-cooled 
with liquid nitrogen, or liquid nitrogen. In order to guard against the possibility 
that the metastable phase itself becomes ferromagnetic at low temperatures, 
the samples were also transferred into iced water after the initial quenching 
and then tested immediately, but the same results were found. 

A polished sample of MnAs was examined microscopically during the 
transformation on heating. No martensitic platelets could be seen, but in- 
stead irregular cracks appeared randomly over the surface. The time taken 
for a crack to form was too short to be observed visually. During the reverse 
transformation the cracks were observed to close up but they still remained 
visible. No new cracks were seen to form in further cycling through the trans- 
formation temperature. Fig. 1 shows a portion of the surface of the specimen 
after several transformation cycles. 

The change of magnetic properties during the transformation was used 
to study the nature of the change within the areas bounded by the surface 
cracks. The sample was covered with a colloidal solution of magnetite and 
microscopical observations were made under dark field illumination. During 
the transformation on heating, the small areas bounded by the cracks in- 
dividually became suddenly darker as the particles of colloid lifted from the 
surface when the grains became paramagnetic. These observations indicated 
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Fic. 1. Polished surface of MnAs after several transformations. Areas A, B, and C are 
the regions investigated by the back-reflection Laue technique. 





Fic. 2. Successive frames of the cine film of the transformation. Interval between the 
pictures 1/16 sec: X100 approximately. 
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that the transformation of a grain and the opening of the cracks surrounding 
it were simultaneous. 

As already noted, the changes also occurred sufficiently rapidly to appear 
instantaneous on visual observation. This was confirmed by cinematography 
at 48 frames per second which showed that areas had transformed completely 
between successive frames. The contrast on prints of these films was un- 
fortunately not sufficiently high to allow reproduction. It was found, however, 
that a slight tilting of the grains occurred during the transformation and that 
the change could be shown under critical oblique illumination conditions (cf. 
Lomer and Pratt 1951-52). Fig. 2 shows the transformation proceeding on 
the surface of a specimen during 8 successive frames taken at 16 frames per 
second. 

It was thus apparent that the rate of transformation of individual grains 
was very rapid, so making any diffusion-controlled mechanism very im- 
probable. To check whether the possibility of a diffusion-controlled mecha- 
nism could be completely ruled out by the rate of change of single subgrains, 
we made further use of the change of magnetic properties to fix an upper 
limit on the time taken for a grain to transform. A specimen of MnAs was 
placed inside a coil in a magnetic field and the signals obtained during the 
transformation were, after amplification, observed on an oscilloscope screen. 
Each individual change took place in a time of about 50 microseconds, which 
was of the same order as the over-all resolving time of the apparatus. We can 
therefore place only a lower limit on the rate of transformation. It could 
conceivably be very much faster. However, this rate is sufficiently high to 
rule out completely a diffusion mechanism for the transformation. 

Over most of the area of polished surfaces of samples of MnAs which had 
been transformed several times, the cracks formed an irregular network (e.g., 
see region A, Fig. 1). In some areas, however, they lay predominantly in one 
direction. An example of such an area is shown in region B, Fig. 1. The orienta- 
tions of the areas marked A, B, C, in this figure, shown in the stereographic 
projection (Fig. 3), were determined by microbeam back-reflection Laue 





Fic. 3. Stereographic plot of the orientations of the three crystals marked A, B, and C 
in Fig. 1. O—A, O—B, O—C. 
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photographs using a beam diameter of about 100 uw. In the areas B and C the 
c axis lay within 5° of the surface of the sample. In both cases also it was found 
that the direction of the ¢ axis coincided with the directions of the cracks. 
The c axis of subgrains in the region A, however, was inclined at about 35° 
to the plane of the microsection. It would therefore appear that the cracks 
run parallel to the ¢ axis in the original grains of MnAs, and the crystals 
fragment during the transformation into pencil-like subgrains. 

This conclusion was confirmed by the microscopical observation of the 
transformation in regions where the subgrains were aligned with their c axis 
parallel to the specimen surface. In these areas the regions which transformed 
individually were long and narrow rather than having the irregular shape 
shown in Fig. 2. 

The back-reflection Laue photographs generally showed asterism. This 
might correspond to a spread of orientations of the individual subgrains in 
the transformed material, since the beam diameter of about 100 uw was larger 
than the size of individual subgrains. For geometrical reasons such misorienta- 
tions can only correspond to small rotations of the singly crystalline subgrains 
around the c axis. In order to test these conclusions the following observations 
were made from back-reflection Laue photographs. 

(a) A specimen which was never cooled below 50° C showed perfectly sharp 
spots with no trace of asterism. 

(6) After the specimen was cooled to room temperature each spot was 
found to have degenerated into several closely spaced spots. 

(c) After several cycles through the transition the spots in both the high 
and the low temperature phases merged into continuous short (asterism) lines. 
In all cases the direction and amount of elongation corresponded to about 2° 
of rotation about the c axis. A Laue photograph taken looking down the c 
axis of subgrains in the area A of Fig. 1 showed almost no asterism, thus con- 
firming this conclusion. 

It is not surprising that such misorientations exist in the repeatedly trans- 
formed samples, since one could hardly expect that the pencil-like subgrains 
would strictly retain their original orientations under the conditions with which 
the transformation takes place. 


DISCUSSION 

The transition at ~40°C in MnAs which appears to be discontinuous, 
diffusionless, and non-martensitic is, as far as we are aware, unique. In so 
far as the expansion of the a axis of the unit cell* does not occur by the rather 
subtle mechanism required of a martensitic change, it might be called a 
“brute force transition.’””’ The change may be contrasted with those found in 
ferroelectrics which appear always to occur by a martensitic process. 

In a martensitic transformation the strain energy is minimized by fitting 
the parent and daughter phases along an unrotated and undistorted habit 
plane, but the mechanism of the transformation also involves some plastic 


*The c axis does not undergo any noticeable change during the transition. The a axis 
expands by about 1% on going to lower temperatures. 





BASINSKI AND PEARSON: TRANSITION IN MANGANESE ARSENIDE 1021 


deformation of the product phase. We believe that MnAs does not transform 
by the martensitic mechanism because such plastic deformation is not possible. 
Instead, the strain energy of the transformation in MnAs is minimized by a 
fragmentation of the crystal grains into pencil-like subgrains. The expansion 
of the @ axis in the transformation of MnAs produces dilatational stresses 
which are maximum across the planes parallel to the ¢ direction. It is therefore 
not surprising that all of the cracks lie parallel to this direction. It is at present 
difficult to visualize the mechanism of transformation inside each ‘‘pencil.”’ 
However, we presume that the strain produced during the transformation is 
absorbed elastically within each subgrain and that the boundary dividing 
the two phases travels along the ‘‘pencil’’ with a high velocity, which perhaps 
approaches the speed of sound. 

Finally we may note that the discovery of a diffusionless non-martensitic 
transformation reduces martensitic transitions to a subclass among diffusion- 
less transitions and the words should no longer be used synonymously. 
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FIELD-INDUCED ABSORPTION IN HYDROGEN! 


M. F. CRAWFORD AND R. E. MacDoNALp? 


ABSTRACT 


This paper reports an investigation of the fundamental rotation—vibration 
absorption induced in hydrogen by an electric field. The integrated absorption 
coefficients of the first four components of the Q branch have been measured 
for a range of field strengths at densities of 83.9 and 44.6 Amagat. The in- 
tegrated absorption coefficient of the S(1) line has been determined at the higher 
density for one high field strength. The relative intensities of the Q components 
are in good agreement with the predicted ratios. Accurate experimental values 
of the matrix element of the mean value of the polarizability, (a)o1, and the 
matrix element of its anisotropy, (y)o1 have been obtained. They are: (a@)o. = 
0.968 X 10-% cm3, and (y)o1 = 0.72X10-% cm® at 83.9 Amagat, and are only 
slightly dependent on density. The frequencies of the lines have been measured. 
Only Q(1) and S(1) show a measurable shift with density. The optical collision 
diameter has been determined and is very small, 0.26A. 


Il. INTRODUCTION 


The strength of an absorption line is proportional to the square of the change 
in the dipole moment associated with the corresponding transition. A homo- 
nuclear molecule such as hydrogen has no dipole moment and absorption 
corresponding to the fundamental rotation—vibration band does not occur. It 
is possible, however, to distort the charge distribution and the rotation- 
vibration lines can then be observed in absorption. This induced absorption 
is produced by collisions, and by the application of an external electric field 
as predicted by Condon (1932). The collision effect was first observed by 
Crawford, Welsh, and Locke (1949) and the field effect by Crawford and 
Dagg (1953). 

The subject of this paper is the field effect and it occurs when an external 
electric field is applied to the molecules in the absorbing path. The electric 
field induces a dipole moment in the molecule, the magnitude of which is 
determined by the product of the polarizability and the field strength. Since 
the polarizability is a function of the internuclear separation and the orienta- 
tion of the molecule the induced dipole will be modulated at the vibrational 
and rotational frequencies. Thus the rotation and rotation—vibration transitions 
can be observed in absorption. The amplitude of the modulated component 
of the dipole moment for a vibrational transition is proportional to the product 
of the rate of change of the polarizability with respect to the normal co- 
ordinate and to the field strength. Therefore, the selection rules are deter- 
mined by the matrix elements of the polarizability, and consequently are 
the same as those governing the Raman effect. 

The field effect is a new technique for determining molecular constants and 
for studying intermolecular interactions. Accurate values for the rates of 
~ IManuscript received April 2, 1958. 
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change of the polarizability with vibrational motion and the anisotropy of the 
polarizability can be obtained directly from the experimentally determined 
line strengths. Inherent in this experiment is the possibility of readily deter- 
mining these over a wide range of conditions. These so-called molecular 
constants are of considerable theoretical significance. For hydrogen many 
calculations of the polarizability and its rates of change are available. Inter- 
comparison of the theoretical and experimental values is a means of assessing 
the adequacy of the wave functions used. Further, accurate values of the 
rates of change are desired for other problems, in particular for the theoretical 
calculation of the collision-induced absorption (van Kranendonk and Bird 
1952). An optical collision diameter can also be evaluated from the integrated 
absorption coefficient and the integrated absorption, and is a function of the 
forces between molecules. 


Il. THEORY 
The integrated absorption coefficient, B’, for light with its electric vector 
in the direction ¢ relative to the direction of the applied field can be written 








| 9 
Sa'vN1| (x : 
(1) pa Yn Pe tn dr} , 
he |e 
where 
P, =a.E, = the component of the induced dipole moment in the 
direction of the light vector; 
E, = the component of the applied electric field in the direction of the 
light vector; 
a, = the polarizability of the molecule in the direction of the light 
vector; 


n,n’ = the totality of quantum numbers required to define the initial 
and the final state; 


l = path length; 

h = Planck’s constant; 

c = velocity of light; 

v = frequency of the transition (where » can be expressed in cm7! 
or sec™!, and B’ has the corresponding unit); 

N  =number of molecules per unit volume in the initial state. 


Equation (1) was formulated by Condon (1932), but the matrix elements 
of the induced dipole moment of the transition for hydrogen were not ex- 
plicitly evaluated. 

The dipole moment, P = aE, is a tensor relation. For a linear molecule the 
polarizability, a, can be expressed in terms of the polarizability parallel to the 
internuclear axis, a), and the polarizability perpendicular to the internuclear 
axis, a,. Thus the dipole moment can be written, 


P, = a Ej, 


2 
( P, a, E,. 
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Fig. 1 shows a molecule in the orientation 6, ¢, in the co-ordinate system in 
which the x axis is the direction of the applied field and the z axis the direc- 


y 





Fic. 1. Molecular co-ordinates. 


tion of propagation of the radiation. The components of the dipole moment 
in the x and y directions are, 
P, = P, cosé+P, sin 6 = E(a, cos*@+a, sin’6), 


P, = P, sin @cos ¢— P, cos @ cos ¢@ = E(a);—a,) sin 8 cos 8 cos ¢. 


(3) 


From equations (3) and the wave function of the vibrating rotator the 
matrix elements in equation (1) for the (v’ = 1, v = 0) vibrational transition 
can be written, 


J vor Re Wis'm’ dr _ E cos | Jor f Few cos’0 Voy’) sin 0 dé do 


(4) + (a, Jou | Po sin’@ Y(s'm’) sin 6d0 as| 
+E sin nor | —_ cos 6 cos @ sin Y¢y-m') sin 6 dé dq, 
where 


(ay )o1 = Sor ay Wi dry; 
(a,)o1 = Sor a, ~idry; 
(y)o1 = (a Jor— (a, Jor} 
VY. sm) = the normalized spherical harmonic eigenfunctions. 


Values of the matrix elements of the Q branch (AJ = 0) and the S branch 
(AJ = +2) for plane polarized radiation with its electric vector parallel to 
the applied field, « = 0, are listed in Table I. Table I shows that (y)oi, the 
matrix element of the anisotropy of the polarizability, is entirely responsible 
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for the absorption in the S branch and contributes to the absorption in all 
components of the Q branch except Q(0). 


TABLE I 
SQUARES OF THE MATRIX ELEMENTS OF THE INDUCED MOMENTS 


Transition (v’ = 1, v = 0) 


The sum over m and m’ of the 
J'<J squares of the induced moments 


Q components 


0—O0 E*Xa)oi? 

1 —]1 E?{3(a@)o12 +(8/75)(y)o1?] 

2<-2 E?|5(a@)o1? +(8/63 )(y)or?] 

3<-3 E?{7(a@)o1? +(8/45)(¥ )o1?] 
S components 

2-0 E*(¥)o1? X 4/45 

3-1 E*(y)o1? X 80/525 





The integrated absorption coefficient B’, equation (1), is proportional to 
the number of absorbers in the initial level and to the square of the matrix 
element of the induced moment. Under the condition of thermal equilibrium 
the relative number of absorbers in the initial level is given by the product of 
the Boltzmann distribution factor e~?7‘%+)/*? and the weighting factors 
introduced by the nuclear statistics and the degeneracy of the rotational 
levels. For hydrogen the odd rotational levels are weighted 3:1 relative to the 
even levels. Since the experimental results show that (y)o:? is about one half 
of (a)o:? its contributions to the squares of the induced moments of the Q 
branch are negligible. This means that the relative intensity of a Q com- 
ponent is proportional to the Boltzmann factor for its initial level times the 
(2J+1) degeneracy of the rotational state times the nuclear weighting factor. 
These calculated intensities are listed in the last column of Table III. 


Ill. EXPERIMENTAL 


In the field effect the observable quantity is an absorption contour. The 
apparatus required for the measurement of this quantity consists of a source 
of infrared radiation, an absorption cell, a means of applying and accurately 
measuring the electric fields, and a high resolution spectrometer. 


(a) Optical System 

Radiation from a 750-watt projection bulb, operated from a Sorenson line 
regulator, traverses a monochromating system and a small band pass in the 
region of 2.4 w is defined by the entrance aperture of the absorption cell. The 
absorption cell contains a parallel plate condenser which also serves as a 
light guide. The hydrogen obtained from commercial cylinders is admitted 
to the pressure cell through a needle valve and brass block. The pressure is 
read from a Bourden type gauge calibrated against a dead weight tester. 

The dipole moment induced in a hydrogen molecule has its major com- 
ponent in the direction of the applied field. An accurate evaluation of the 
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absorption coefficient therefore requires a determination of the degree of 
polarization of the radiation. This is effected by introducing a Wollaston 
double image prism between the exit aperture of the pressure cell and the 
entrance slit of the spectrometer. This prism gives an angular separation of 
about 0.04 rad between the two images of the exit aperture polarized in per- 
pendicular directions. One image or the other is focussed on the entrance slit 
of the spectrometer by a condensing lens. 

The Littrow type spectrometer has a plane grating as the dispersing element. 
The grating, obtained from Bausch and Lomb, has 300 lines per mm, a blaze 
angle of 26° 45’, and a ruled area 204 mm wide and 130 mm high. The auto- 
collimating mirror is a 13° off-axis paraboloid, 8 cm diameter, 100 cm focal 
length. The limiting aperture is set by the collimator. 

For several reasons the signal intensity available at the detector is relatively 
low. The small half-widths of the observed transitions, of the order of 0.13 cm=!, 
impose a restriction on the spectral slit widths that can be used if a reasonable 
peak absorption is to be recorded, and the signal at the detector is proportional 
to the product of the entrance and exit slit widths of the spectrometer. The 
available signal is further limited by the presence of the collision-induced 
absorption which is proportional to the square of the density, by the low 
efficiency of the light guide, and by the fact that the radiation is restricted to 
one plane of polarization. Accordingly, stringent requirements are placed on 
the electronic amplifying system. 

(b) Detector and Amplifying System 

A lead sulphide photoconductive cell is the detector. This cell has a sensitive 
area } in. X # in. and a resistance of 0.6 meg at room temperature. When the cell 
is cooled with a mixture of dry ice and acetone the resistance reaches a ter- 
minal value of 20 meg in one half hour. 

The block diagram, Fig. 2, portrays the main components of the amplifying 
system. Radiation incident on the detector is periodically interrupted by a 
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Fic. 2. Block diagram of the amplifying system. 
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mechanical chopper placed at the entrance slit of the spectrometer and the 
chopping frequency is 150 cycles/sec. A signal generated at the detector is 
applied in succession to a preamplifier, a twin T amplifier, and a phase detector. 
A reference signal from a photocell circuit, chopped at the same frequency 
as the wanted signal, is applied to the suppressor grids of the push-pull stage 
in the phase detector and ensures that out-of-phase noise voltages are rejected. 
The unidirectional output from the phase detector is filtered and the d-c. 
signal displayed on a 50-millivolt Brown “‘Elecktronik’”’ recorder. The ampli- 
fying system with a gain of 50,000 gives a full scale deflection for the Johnson 
noise of the input resistance. The twin T band pass filter, tuned to 150 cycles 
sec, is adjusted to a band width of approximately five cycles/sec to prevent 
overloading the phase detector with noise. The phase detector gives an effective 
band width variable from 3 to 3 cycle/sec and this is independent of varia- 
tions in the chopping frequency. The electronically regulated power supplies 
are fed from a Sorenson line regulator which is powered by an isolating 
line transformer. The isolating line transformer is required because in this 
laboratory the 60 cycle a-c. neutral line is a very poor ground. The PbS cell 
when cooled is operated with a polarizing voltage of 30 v supplied by a shielded 
B battery. 

(c) Measurement of the Field 

The absorption cell is that used previously in the investigations of the 
collision-induced absorption (Crawford, Welsh, and Locke 1949). To establish 
the electric field within the absorption cell a source of high voltage is applied 
through a spark plug to one side of a parallel plate condenser. The other plate 
of the condenser is contacted through a spring clip to the steel wall of the 
pressure cell which is grounded externally to a water pipe. 

The parallel plate condenser consists of two plane strips of glass 82.4 cm 
long, 1.2 cm wide, 0.3 cm thick, coated with silver and mounted in two Lucite 
spacers as shown in Fig. 3. The grooves in each spacer were milled simul- 
taneously with two cutters, and are slightly wider than the thickness of the 
glass plates. The plates are held in the grooves by set screws bearing against 


a 





Fic. 3. Cross section of the condenser. 
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strips of brass shim stock inserted between the outer surfaces of the plates 
and the outer walls of the grooves. The plates were aligned by adjusting the 
set screws and their parallelism checked by measuring their separation with 
a travelling microscope at many positions along the top and bottom of the 
assembled unit. The mean deviation of the separation did not exceed 0.01 mm. 
The condenser with its highly reflecting plate also functions as a light guide. 

Since the absorption coefficient is proportional to the square of the electric 
field strength, the field must be maintained constant and its strength measured 
accurately. The field was produced by a power supply obtained from Neutronic 
Associates. It is regulated to 0.1% and its output can be varied continuously 
in the range 5-30 kilovolts. The potential difference between the condenser 
plates was determined by connecting in parallel with the plates a voltage 
divider consisting of 14.99+0.02 meg resistor in series with a 1000-ohm 
standard and measuring the voltage across the standard with a potentiometer. 


(d) Measurement of Integrated Absorption Coefficients 
For the field effect we define a nominal absorption coefficient k(v) by the 


relation 
(5) I(v) = To(ve*", 
where 
Io(v) = the signal in the absence of the electric field; 
I(v) = the signal when the field is applied; 
l the path length; 
I/Iy = T(v) = the transmission. 


The absorption is defined as 


(6) re f ~ (1—T(s))ee. 


The distribution of Jo(v) with frequency is a function of the collision-induced 
absorption, the spectral distribution of the radiation source, and the response 
of the spectrometer. In Fig. 4 the dotted line represents the distribution of 
I(v) over part of the hydrogen fundamental region. The dip in the absorption 
contour at 4155 cm™! divides the Q branch of the collision-induced absorption 
into two parts and each has a half-width of ~250 cm-!. With the application 
of the electric field four lines appear which are not present in the zero field 
trace. These lines, labelled Q(0), Q(1), Q(2), Q(3), are components of the Q 
branch, AJ = 0, of the fundamental vibrational transition (v’ = 1<—v = 0); the 
initial value of J is designated by the number in parentheses. The field effect 
lines are exceedingly sharp and have a half-width which is mainly instrumental 
even at a pressure of 1500 p.s.i. By evaluating fR(v)ldv = f log -(Io/I)dv over 
the absorption contour of a line the nominal integrated absorption coefficient, 
B, is obtained. This nominal integrated absorption coefficient in general differs 
from the integrated absorption coefficient, B’, unless the resolving limit of the 
spectrometer is appreciably less than the half-width of the transition proba- 
bility. However, it will be shown that B’ can be obtained from B. 
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Fic. 4. The Q components of the fundamental vibrational band. 


B was determined from a point-by-point scan of the absorption contour. 
The grating was set at a fixed angle in the vicinity of the absorption; the field 
was turned on and off several times to ensure a good average value for 
log (Jo/I). The grating was then rotated through a small angle and the pro- 
cedure repeated. The total angle through which the grating turned in travers- 
ing a line was about 5 X 10~‘ rad. Consequently, the angular separation between 
adjacent settings of the grating had to be measured to within 10~* rad to 
obtain a precision of a few per cent in the integrated absorption coefficient. 
The driving mechanism of the grating did not permit a determination of the 
angular change with this precision. 

The required precision was obtained by an interferometric method. The 
interferometer consists of two flats and one face of each is coated with a 
dielectric film of high reflectivity in the green region of the spectrum. One of 
the flats is mounted so that it moves with the grating while the other remains 
stationary. With a small separation between the coated surfaces, narrow 
localized fringes are visible when the interferometer is illuminated with 
collimated monochromatic light. The angular separation between adjacent 
settings of the grating was determined from the number of fringes that crossed 
a fiducial mark on one of the plates during the rotation of grating from one 
setting to the next. The angular change is given by (8), which is derived from 


(7) NX = 2t cos (r+8), At = 4.46, 


where 


N = order of interference; 
A = wavelength; 
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t = separation between the inner faces of the flats at the fiducial mark; 

ry =angle of refraction measured at the air space between the two 
flats; 

6 = angle between the inner faces of the flats; 

Aé@ = change in 6 due to rotation of grating; 

d =distance from the axis of rotation of the grating to the fiducial 
mark. 


The angles r and @ were adjusted to approximately zero. For this condition a 
small change in @ due to rotation of the grating changed ¢, but had negligible 
effect on the cosine term. Therefore 


(8) Ad = X.AN/2d, 


where A.V is the number of fringes that crossed the fiducial mark. Angular 
rotations of 2X10-> rad were measured by this method to within 2107 rad. 


IV. RESULTS 


As shown by the discussions of Dennison (1928), Nielsen, Thornton, and 
Dale (1944), and by Kostkowski and Bass (1956) the nominal integrated 
absorption coefficient and the integrated coefficient are in general unequal 
but related. The relation involves the half-width of the transition probability 
under the experimental conditions, the resolution of the spectrometer, and 
the peak absorption. Under two conditions the nominal integrated absorption 
coefficient equals the integrated absorption coefficient. These are: (1) when 
the instrumental half-width is small compared to the half-width of the transi- 
tion probability, (2) when the absorption is small enough so that only the 
linear term in the expansion of the exponential is required. Since the experi- 
mental conditions were not compatible with (1), the integrated absorption 
coefficients were obtained under condition (2) by extrapolating the nominal 
coefficients to small absorption. 


(a) Nominal Integrated Absorption Coefficients—Q Components 

The nominal integrated absorption coefficients were determined under 
various experimental conditions using plane polarized light and « = 0. Values 
of B for four field strengths are recorded in Table II for the (1 — 1) component 


TABLE II 
NOMINAL INTEGRATED ABSORPTION COEFFICIENTS—(Q(1) COMPONENT 
Transition (v’ = 1, v = 0) 





B, cm." 
B, cm™ B/E*X10', cm71/ (calculated 
E*X10-4(e.s.u./em)? (measured) (e.s.u./em)? eq. (10)) 
6.35 271 42.7 .273 
4.77 .229 47.8 .227 
2.82 .151 53.6 .150 
1.81 .102 56.5 -101 





For a: absorption, A = 0.248 cm™!; minimum transmission, Tmin = 62%. 
Conditions: 1.98 mm spacer; density 83.9 Amagat; 27° C; 82.4cm path length; 0.46 cm slit 
width, 
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of the Q branch. Representative contours from which these nominal absorp- 
tion coefficients were evaluated are shown in Fig. 5. Each B value is an average 
of at least two traces and each point on a contour, Fig. 5, is itself an average 
of a minimum of three measurements of log,(Jo/J). An average deviation of 
about 1% for the smaller coefficients is a realistic measure of their precision. 


0.50 








= 4 2 
040 E*=6.35x 10*(esu/cm) 


7 £2477 xio%esu/cm)? 


€7=2.82x10*(esu/cm)* 
0.20 


£7=1.81x10*(esu/cm)* 


0.4 0s 12 1.6 
cm! 


84 AMAGAT Q(1) 


Fic. 5. Contours of the Q(1) components, 


Theoretically the integrated absorption coefficient, B’, is proportional to 
the square of the field strength. It is evident from column 3 of Table II that 
the nominal integrated absorption coefficient departs progressively with in- 
creasing field strength from the predicted constant relationship. In Fig. 6 a 
plot of B vs. E? illustrates this deviation from linearity. It is anticipated from 
condition (2) that when the resolving limit is larger than the transition half- 
width of the component the departure from linearity will be a function of 
the integrated absorption coefficient. Therefore B can be represented by a 
power series in E?, 

(9) B = ak’?+q,E'. 
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Fic. 6. Nominal absorption coefficient vs. field strength squared. 


Because of the small deviations associated with the experimental results 


graphical methods are adequate to determine the constants. Fig. 7 is a plot 
of B/E? vs. E*. This linear relation gives 


ay) = 61.5X 1077 cm—!/(e.s.u./cm)?, 


(10) 
a, = —2.85X10-"! cm-!/(e.s.u./cm)4. 
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Fic. 7. B/E* vs. E?. 


Substitution of these values of ap and a; in equation (9) yields column 4 of 
Table II. Thus for the experimental conditions noted in Table II, the depen- 





Vea 
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dence of the nominal integrated absorption coefficient on the field strength is 
well represented by the equation 


(11) B = 61.5X107 E?—2.85X10-"E*. 


The coefficient of the first term on the right hand side of equation (11) is the 
value of B’, the integrated absorption coefficient, per unit field strength 
squared. 

To appraise the accuracy of do, the slit width and the separation of the 
plates of the condenser were varied. The spectral slit width was increased 
from 0.46 to 0.86 cm~ and the plate separation decreased from 1.98 to 1.37 mm. 
It was somewhat surprising to observe that with the application of the electric 
field an absorption was detected over the whole of the spectral range trans- 
mitted by the monochromator (0.2 yu). Further tests showed that the effect 
persisted and with about the same magnitude (6% absorption for 15 kilo- 
volts), even beyond the region of the collision-induced absorption. This anoma- 
lous absorption was evident in nitrogen to about the same degree as in hydro- 
gen and was just detectable in the case of helium gas. Returning to hydrogen, 
the dependence of the anomalous absorption on the field strength was found 
to be linear. With the Wollaston prism set to select only the component of 
the radiation with the electric vector perpendicular to the applied field, the 
absorption recorded was only about one third that of the parallel component 
for the same experimental conditions. The effect was not detectable when 
the larger spacer, 1.98 mm, was used, even though an equal field was applied. 
It appears that the field slightly decreases the reflectivity of the silvered con- 
denser plates. A slight decrease in reflectivity has a greater effect on the signal 
when the narrower spacer and wider slit are used because a larger portion of 
the signal is then due to radiation multiply reflected in the light guide. 

By measuring the absorption well outside the region of the Q(1) line and 
applying this as a constant correction to the absorption recorded in the region 
of the Q(1) line, the result obtained was 


(12) ayo = 60.0X 1077 cm—!/ (e.s.u./em)?. 


The agreement between the values of a in (10) and (12) is very good. There 
appears to be little evidence for a systematic error in the results associated 
with the slit width or the condenser plate separation. 

Another value of ao was obtained from a series of measurements with the 
1.98 mm spacer and the grating in a different geometrical arrangement. The 
value of ay obtained was 61.0X 1077 cm~!/ (e.s.u./em)*. 

The nominal integrated absorption coefficient of the Q(0) component was 
measured at two field strengths and that of the Q(2) and of the Q(3) at only 
one high field strength. Since the absorption in these components is small 
it was not necessary to make extensive measurements at a number of field 
strengths as required for the Q(1) component. From the nominal integrated 
absorption coefficients of these three components their integrated absorption 
coefficients can be obtained by using the ratio of the coefficients of the powers 
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of E as determined in the analysis of the data for the Q(1) component. Equation 
(9) can be rewritten 
(13) B = aoE?+ (a1/a0*) (aoE?)?. 


Since the absorption is small the second term on the right hand side of equa- 
tion (13) is small and aE? can be approximated by B, giving 


(14) Bis aoE? + (a;/a,?) B?, 


where the ratio a;/a,° is obtained from the analysis of the Q(1) data (10). 

The adequacy of equation (14) was tested and proven by applying it to 
the nominal integrated absorption coefficients of the Q(0) component obtained 
at the two field strengths. It was also used to determine the vaiue of ao for 
the other two components. 

The average value of do with its average deviation for each of the four 
components of the Q branch is listed in Table III, column 2. The ratios of 
these experimental values are given in column 3 and the ratios of the theoretical 
values in column 4. The good agreement between columns 4 and 5 shows that 
the experimental measurements and their analysis give accurate values of the 
integrated absorption coefficients. 


TABLE III 
INTEGRATED ABSORPTION COEFFICIENTS—Q COMPONENTS 
Transition (v’ = 1, v = 0) 
(83.9 Amagat; 27° C; 82.4 cm path) 





Ratio 
J'’<J ae X10? cm=!/(e.s.u./em)? Experimental Theoretical 
0<—O0 12.0 +0.2 1.000 1.000 
1-1 61.0 +0.3 5.08 6. i1 
2<2 11.1 +0.3 0.925 0.915 
3-3 8.15+0.2 0.679 0.705 


(b) Dependence of the Absorption Coefficient on Density 
B of the Q(1) component was measured at a lower density, 44.6 Amagat, 
and its dependence on field strength was found to be 


(15) B = 31.3X107 E°—0.984X10-" Et. 


From equation (15) ao equals 31.3107 cm~!/(e.s.u./cm)*. The ratio of 
do to the density at a density of 44.6 Amagat is 0.702; at a density of 83.9 
Amagat (Table III) the ratio is 0.727. This 4% difference in the ratios indicates 
that the rate of change of the polarizability increases slightly with increasing 
density. In view of the accuracy of the experiments, it is highly probable that 
this increase with density is real. It is qualitatively consistent with the obser- 
vations of MacNaughton (1948) on the variation of Raman scattering per 
molecule with the density of the scattering medium, but the magnitude of 
the effect in hydrogen is too small relative to the experimental precision to 
warrant quantitative deductions. 
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A small correction should be applied to ao because the local fields at the 
absorbing molecule are slightly different than the static and radiation fields. 
This correction for hydrogen, neglecting fluctuations, reduces do by less than 
2% for the higher and 1% for the lower density (van Kranendonk 1957). 


(c) Anisotropic Absorption—S(1) Component 

The matrix element of the anisotropy of the polarizability, (y)o1, can be 
determined more accurately from the S branch than from the Q branch. 
The contribution of (y)o1 to the absorption is quite small in both branches, 
but it is the sole contributor in the S branch for any polarization. Most favor- 
able measurements can be made on the S(1) component because at room 
temperature it is the strongest. The absorption was first detected using a high 
field strength obtained with a small plate separation. Because the anomalous 
absorption observed with a small plate separation introduces a correction of 
the same order as the absorption due to (y)o1, better accuracy was obtained 
with a large plate separation even though the absorption was smaller. 

The absorption was measured with the light vector parallel to the applied 
field and six determinations of the integrated absorption coefficient were 
made at high field strengths. Since the peak absorption was only 2}% the 
nominal integrated absorption coefficient equals the absorption coefficient. 
For a density of 83.9 Amagat, a path length of 82.4 cm, and a plate separation 
of 1.98 mm, the value obtained is 


(16) dy = 1.89+0.11X 1077 cm~'!/(e.s.u./cm)?. 


(d) Frequencies of the Absorption Lines 

The measurements by Crawford and Dagg (1953) on the components of 
the Q branch indicated a frequency shift, the magnitude of which was depen- 
dent upon the J value involved in the transition and the density. More accur- 
ate measurements of the frequencies were made in this investigation to test 
the dependence indicated by their preliminary results. 

The components of the Q branch were measured by superposing an argon 
emission spectrum on the absorption spectrum and recording both on the 
same trace. Seven argon lines conveniently bracket the second order of the 
Q components and are suitable as secondary standards (Plyler, Blaine, and 
Tidwell 1955). The maximum separation between adjacent standards is 
150A in NV) and interpolation within this interval is precise. The value of the 
refractive index required for converting wavelength to frequency was obtained 
from the formula proposed by Edlén (1953). The frequency of each Q com- 
ponent was determined from six scans and had an average deviation of 0.05 
cm~! and a maximum deviation of 0.1 cm~!. Measurements were made at 
densities of 83.9 and 44.6 Amagat for the Q lines. The frequency of the S(1) 
line was measured relative to neon standards, but the small absorption per- 
mitted measurements at the higher density only. 

The results of our frequency determinations and the recent low density 
Raman effect measurements of Stoicheff (1957) are contained in Table IV. 

Both the present static field and Raman frequencies quoted in Table IV 
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TABLE IV 
COMPARISON OF FIELD EFFECT AND RAMAN EFFECT FREQUENCIES 
Transition (v’ = 1, v = 0) 
(units cm~?) 





Field effect 








Component — —_——_—_______—__—_——— Raman effect* 
J'<J 83.9 Amagat 44.6 Amagat 1.5 Amagat 
0-0 4161.20 4161.21 4161.13 
1-1 4155.02 4155.13 4155.20 
2<2 4143.43 4143.44 4143.39 
3-3 4125.93 4125.90 4125.83 

(4712.86 
3<1 4712.80 4713 .08° 


(4713.01¢ 








*Stoicheff (1957). 
’Cumming (1952). 
©Feldman (1955). 


are accurate to 0.05 cm~'. Intercomparison of the measured frequencies in- 
dicates a small, density-dependent shift of the Q(1) component to lower 
frequency. The dependence of the shift on density is linear within the limits 
of error. The frequencies of the other three components of the Q branch are 
all greater than the corresponding Raman frequencies, but the shifts are less 
than the experimental errors and are not density-dependent. Only the S(1) 
component of the S branch was measured in the static field experiment. A 
measured value of the frequency of the S(1) component is not contained in 
the paper by Stoicheff, but from the listed molecular constants a value can be 
calculated. The frequency of this line has been measured in the Raman 
effect of the gas at low density by Cumming (1952) and by Feldman (1955). 
Their values are 4713.08 and 4713.01 cm™, respectively. Comparison of the 
three Raman effect values and the field effect value listed in Table IV shows 
that the S(1) line at a density of 83.9 Amagat is shifted in the same direction 
as the Q(1) component and by about the same amount. Thus the static field 
results establish that only the Q(1) and S(1) components show a density 
dependent frequency shift and even for these components it is very small. 

Raman effect measurements in the gas at high density are not available for 
comparison, but Raman results for liquid and solid hydrogen (Allin, Feldman, 
and Welsh 1956) with an ortho/para ratio equal to the equilibrium ratio at 
room temperature show for the Q(1) line a shift of about 11 cm and in the 
same direction as that recorded for the gas in the field effect. But the shift is not 
linear with density up to the density of the liquid. 


V. INTERPRETATION 
(a) Matrix Elements of the Polarizability 
The anisotropic and isotropic parts of the matrix elements of the polariza- 
bility can now be determined from the experimental data. Inserting the 
squares of the matrix elements listed in Table I and the appropriate values 
of ao, equation (1) becomes for the S(1) component 
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2 _ 84 yNI80(y)n _ —7 1 2, 
B/E’ = cs 505 0 = 1.89 107° cm™"/(e.s.u./cm)*; 
and for the Q(1) component 
3 oxy 
B/E’ = SNF set Soni | = dy) = 61.0K10™ cm~?/(e.s.u./cem)’, 
where 
N = 49.7 X 10"/cm, 
l = $2.4 em, 


»(S(1)) = 4713 cm-, 
v(Q(1)) = 4155 cm=!, 


From these equations 
(y)o. = (0.72+.02) X 10-* cm, 
and 


(a)or = (0.968+.01) X 10-* cm?. 


These values are for the gas at a density of 84 Amagat. The density depen- 
dence discussed in Section IV(d) indicates that at zero density (a): is smaller, 
but by less than 2%. 

The results of at least eight theoretical calculations are available for com- 
parison with the experimental values of (@)o1 and (y)o1 and are summarized 
in a paper by Ishiguro, Arai, Kotani, and Mizushima (1952). Since their 
calculations are the most elaborate and take into consideration both the 
electrical and the mechanical anharmonicities our results are compared with 
their predictions in Table V. The experimental values of (@)o: and (y)o1 are 
appreciably lower (30%) than the theoretical predictions although the theo- 
retical and experimental ratios of (ay)o1/(a@,)o1 are in much better agreement. 
An accurate calculation of second order molecular constants such as the rates 
of change of the polarizability places stringent requirements on the theory, 
and the comparison in Table V shows that even for hydrogen the best calcula- 
tion to date is not very satisfactory. 


TABLE V 
EXPERIMENTAL AND THEORETICAL VALUES OF THE MATRIX ELEMENTS 
Transition (v’ = 1, v = 0) 











Field effect Theoretical* 
(a)o1 0.968 XK 107% cm? 1.39X10-% cm? 
(yo 0.719 0.901 
(a\ or 1.45 1.99 
(a, Jou 0.73 1 .09 
(ay or /(ay Jor 1 .99 1 83 





“Ishiguro, Arai, Kotani, Mizushima (1952). 


Our value of (@)9: can also be compared with the preliminary value reported 
by Crawford and Dagg (1953), and with the value obtained from the measure- 








1038 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


ments of the intensity ratio of Rayleigh to Raman scattering (Stansbury, 
Crawford, and Welsh 1953). The preliminary value is appreciably higher, 
1.2(5) X 10-* cm; the difference is due mainly to the inaccuracy of the method 
used in the earlier experiment to measure small angular rotations of the 
grating. The value of (a@)o: derived from the intensity ratio experiment is 
1.1 10-*° cm* compared to 0.968 X 10-** cm’. This agreement is satisfactory 
in view of the much larger error quoted for the scattering experiments. 


(b) Optical Collision Diameter 

The absorption lines of hydrogen even at 1435 p.s.i. are very narrow, sym- 
metrical, and not shifted appreciably. The interaction between molecules is 
small and one is justified in applying the classical theory of impact broadening 
to derive an optical collision diameter. The following treatment is based on 
the paper by Dennison (1928). 

With the transmission defined by 





a kyl 
(17) Lome ee 
Dennison shows that for a spectrometer of infinite resolution 
(18) bd os 1.353 B'No*/(Bm)"" ee a 
”  (7.60N°o'/42Bm)+(v—)? Bx?’ 

where 

a = the optical collision diameter; 

m = the mass of the molecule; 

B= Per: 

B’ = the integrated absorption coefficient; 

N = the number of molecules per cm’. 


For a spectrometer of finite resolution the area under the absorption curve 
is 

(19) A= { (1—T»)d», 

independent of slit width. Using equations (17), (18), equation (19) can be 


written 
(20) dm f Ot eee 


o is readily obtained by substituting the experimental values of B’ and A 
in an expanded form of (20). For the higher density, B’ = aoE? = 61.0X6.35 
<X10-? = 0.387 cm—!, A = 0.248 cm, Trin = 62% (Tables II, III), and 
o = 0.24X10-8 cm. For the lower density (Section IVd), B’ = 31.3X6.35 
X 10-3 = 0.199 cm-!, A = 0.145 cm, Tam = 76%, and o = 0.27X10-° cm. 
This agreement is satisfactory, and 0.26X10-° cm is taken as the optical 
collision diameter. 

The value of o can be checked by calculating minimum transmissions, 
Tm, and comparing these with the measured values. When the resolving 
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limit of the grating is appreciably smaller than the spectral slit width a tri- 
angular distribution function is a good representation of the instrumental 
intensity distribution. For this condition 


ae 2 " By 
(21) Tata = 2 J faecal te. 


where 


x = (v—v 9); 
vo = frequency at the center of the line; 
w = spectral slit width. 


Using ¢ = 0.26X10-8 cm, values of Tian were calculated for a sequence of 
values of w. The slit widths 0.43 cm and 0.40 cm! give transmissions in 
agreement with the measured values for the higher and lower densities, re- 
spectively. The spectral slit width determined from the measured widths of 
the entrance and exit slits was 0.46 cm~!. The difference is the order of the 
error in the direct measurement of the slit width. Thus the observed contours 
give a consistent value of the optical collision diameter. For the experimental 
conditions listed in Table II the transition half-width obtained from equation 
(18) by substituting ¢ = 0.26X10-' cm is 


The small value of the collision diameter, which is about one third the inter- 
nuclear separation, is qualitatively consistent with the interaction energy 
between two hydrogen molecules and shows that only S-type scattering is 
effective in producing broadening. 
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NUCLEAR LEVEL SPACINGS! 


A. G. W. CAMERON 


ABSTRACT 


T. D. Newton’s shell-dependent level spacing formula depends upon the 
densities of single-particle orbits near the Fermi level in the nucleon gas. It has 
been found that these densities can be computed by taking appropriate averages 
of the second differences between adjacent atomic masses, using the writer’s 
semiempirical mass formula from which electrostatic and pairing energy terms 
have been omitted. With this procedure, observed nuclear level spacings have 
been fitted with a roct mean square error factor of 1.74. This fit shows that the 
level spacings are proportional to (2J+1)7! to a good approximation. Since the 
average density of single-particle orbits depends on the nuclear excitation energy, 
and since its computation takes a long time even with an electronic computer, an 
approximation formula with five coefficients has been fitted to the computed 
orbit densities for each Z and N in the ranges 10 < Z < 109and 10 < N < 189. 
These coefficients are tabulated. 


1. INTRODUCTION 


For many years it has been customary to represent nuclear level spacings 
by a formula of the type 


(1) D = c, exp—C2E!”, 


where D is the spacing of levels of a given spin and parity at excitation energy 
E, and c, and cz: are constants which depend upon the nucleus. With empirically 
determined values of c; and cs, this formula could be expected to give an 
order-of-magnitude estimate of level spacings in nuclei not too near closed 
shells of nucleons. Recently Newton (1956) has shown that a more elaborate 
formula of this type, derived for a Fermi gas model of the nucleus, can fit 
nuclear level spacings with a root mean square error factor of 3, provided 
one replaces the density of nucleon orbits at the Fermi level by the density 
of nucleon orbits in the shell model (Klinkenberg 1952). 

Although Newton’s achievement represents an enormous increase in 
accuracy over previous level-spacing formulas which do not depend on shell 
effects, nevertheless the remaining error factor of three is too large for the 
rather precise calculations required in nuclear astrophysics (Cameron 1957a). 
Therefore, the writer undertook the present investigation in an attempt to 
improve Newton’s level-spacing formula following a hint from Newton that 
the density of nucleon orbits at the Fermi level might be better calculated 
using nuclear masses instead of Klinkenberg’s shell-model energy spacings. 


THE LEVEL-SPACING FORMULA 


The derivation of the level-spacing formula from the Fermi gas model has 
been given in detail by Newton (1956). According to this derivation one can 


1Manuscript received November 13, 1957. 
Contribution from Atomic Energy of Canada Limited, Chalk River, Ontario. 
Issued as A.E.C.L. No. 643. 
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write the density of levels with magnetic quantum number M at excitation 
energy U in a nucleus with Z protons and N neutrons as 


2 2 

: addi eee aa a ae 
(2) P(UZNM) = (2) “(fifefafs)” “exp 3 Gt 2 mer" 
where 

ti _ Gt, 

fe = Gyt, 

fs = mGt, 

z r x 

= —GHh— -P, 
has 2 
G = Gz+Gy, 


t is the nuclear temperature, Gz is the density of proton orbits at the Fermi 
level, Gy is the density of neutron orbits at the Fermi level, and 7? is a weighted 
average of the squares of the magnetic quantum numbers of the nucleons at 
the Fermi level. The density of levels of spin J is 


(3) D-(UZNJ) = P(UZN, M = J)—P(UZN, M = J+1). 


The nuclear temperature is 


_ 3[3,(9 i 
(4) t= S18 3 ; 


This fairly complicated expression has been used in these calculations rather 
than the simplified expression given by Newton (1956), in order that the 
level-spacing formula (2) will not approach zero at zero excitation energy, an 
unattractive feature of Newton’s formula in many calculations. Although the 
derivation of the level-spacing formula does not hold for small values of U, 
the formula does go smoothly to a finite value at U = 0 when the nuclear 
temperature is given by equation (4). 
Following Newton, the following two approximations are made: 


(5) m? = const. A*’, 
(6) eg ee a 
Wlaal "| «a 1” sa“G ~ 


The level-spacing formula can thus be written: 


2 1/2 2 
(7) D(UZNJ) = const. er+1y| Gea'a'(® -3)/'| exp— Gt. 


This formula has one adjustable constant, provided Gz and Gy can be 
accurately determined. Newton used a second adjustable constant in his 
procedure for determining values of Gz and Gy. 


3. THE DENSITIES OF NUCLEON ORBITS 


Consider a nuclear potential well containing a set of degenerate nucleon 
orbits. As one adds nucleons to these orbits, to a first approximation, the 
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nucleons all have the same binding energy. If we now destroy the degeneracy 
and introduce a spacing d between the energies of successive orbits, then to a 
first approximation the difference between the binding energies of successive 
nucleons will also be d. These differences between successive nucleon binding 
energies are equal to the second differences between the masses of the nuclei 
so formed. This simple model forms the basis of the work reported here. 

The writer has recently published a revised form of the semiempirical 
atomic mass formula (Cameron 19575), in which shell and pairing corrections 
have been determined by an empirical procedure. The mass formula contains 
volume and surface energies and Coulomb and Coulomb exchange energies. 
Not all these terms are desirable in the proposed task of determining nucleon 
orbit spacings from second differences of masses. In the process of forming 
excited nuclear states, nucleons are raised from occupied orbits to unoccupied 
ones, and the energy differences involved are not the same as the binding energy 
differences which result when additional particles are added to the nucleus. 
For this reason it was decided to drop the Coulomb and Coulomb exchange 
energy terms from the mass formula and also the pairing energy terms. When 
a proton is raised to an excited orbit the electrostatic energy of the nucleus is 
unchanged to a good approximation. The pairing energy is a collective pheno- 
menon which has no direct relation to the spacing of nucleon orbits. 

The modified mass excess formula can then be written as: 


(8) M'(Z, A) = 8.3674 —0.783Z+E,+E£,+S(Z, N), 
where 
Ate 
E,=- 17.0354 1—1.846 Lf ae | A, 
[ A—97\2 2 2 ; 
E, = 25,8357) 1— 1.712 424). |e , 


and S(Z, V) is a shell correction energy which was empirically determined. 
The above terms are expressed in Mev energy units. The shell and pairing 
energies were written as 

(9) S(Z, N) = S(Z)+S(N), 

(10) P(Z, N) = P(Z)+P(N). 


P(Z) and P(\) are zero for odd values of Z or .V, respectively. The quantities 
S(Z)+P(Z) and S(.V)+P(N) have been tabulated (Cameron 1957). The 
pairing corrections P(Z) and P(.\V) are given here in Table I. 

Newton (1956) has shown that a fictitious excitation energy U must be 
used in the level-spacing formula. This fictitious energy is based on the energy 
the ground state would have in an odd-odd nucleus where there is no pairing 
energy. We therefore write 


(11) U = E+P(Z, N), 


where E is the actual excitation energy and P(Z, N) is the sum of appropriate 
negative energy terms in Table I. 
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TABLE I 


PAIRING ENERGIES FOR NEUTRONS AND PROTONS 
(If either nucleon number is odd, the pairing energy is zero for that nucleon) 











ZorN  P(Z) P(N) ||\ZorN  P(Z) P(N) ||ZorN  P(Z) P(N) 
10 —2.25 —2.70 || 74 —1.23 1.20 |} 188 —0.80 
12 —=2.10 9:48. || 76 —0.85 —0.94 140 —0.68 
14 =9.15 Buh 78 —0.98 —1.00 142 —0.64 
16 —1.54 —1.74 || 80 —0.72 —1.05 || 144 —0.58 
18 — {42 195 82 —0.80 —0.54 146 —0.55 
20 —1.51 —1.72 || 84 =0.77 —0.60 148 —0.57 
22 1,73 —1.63 || 86 —0.89 —0.75 || 150 —0.57 
24 —],44 —1.41 || 88 —0.92 —0.75 || 152 —0.55 
26 —1.45 —1.29 || 90 —0.80 —0.85 || 154 —0.60 
28 =—1:37 —1.47 || 92 —0.81 —0.97 || 156 —0.58 
30 —-1.09 —1.32 || 94 —0.69  —1.02 || 158 —0.58 
32 — a — 546 96 —0.70 —1.05 160 —0.61 
34 —1.42 —1.44 || 98 —0.76 —1.06 162 —0.63 
36 = ea —1.46 || 100 —(079 1.07 164 —0.65 
38 —1,20 —1.52 102 —0.80 —1.06 || 166 —0.66 
40 —1.00 —1.51 104 —0.74 —1.05 168 —0.65 
42 —J.16 —| 47 106 —0.73 —1.02 170 —0.65 
44 = 98 —1.45 108 = 72 —0.97 172 —0.64 
46 —1.38 =], 28 | 110 —0.72 —0.91 || 174 —0.64 
48 $58 —i 23 112 —0.72 =f} 33 176 —0.63 
50 2.22 5 SF 114 —0.71 =), 74 178 —0.61 
52 = 04 —0.62 || 116 —0.69 —0.66 180 —0.59 
54 —1.11 —0.76 || 118 —0.68 —~0; 61 182 —0.55 
56 —1.13 —1 a 120 —0.66 —0.61 || 184 —0.39 
58 =i, 21 —1.23 || 123 —0.61 —0.90 186 —0.36 
60 —1.43 —1.40 124 —0.42 —(0.52 188 —0.40 
62 —1.15 —1.36 || 126 —0.36 —0.81 || 190 —0.40 
64 —0.99 —1.30 128 —0.41 —0.68 || 192 —0.40 
66 —0.91 —1.29 || 130 —0.49 —0.72 || 194 —0.40 
68 —0.92 —1.%4 132 =O. 27 196 —0.40 
70 —1.00 —1.28 || 134 —0.68 198 —0.40 
72 —1,11 —1.24 136 —0.67 || 200 —0.40 











We now define the orbit separation energies as 
(12) d(Z) = [M'(Z, A)-—2M'(Z—-1, A—1) + M'(Z—2, A—2)],. 
(13) d(N) = [M’(Z, A)—2M'(Z, A-—1)+M'(Z, A—2)],. 


Here the bracket with subscript v indicates that the quantities must be 
evaluated near the valley of beta stability where the shell corrections were 
determined (these expressions vary slowly as one crosses the valley of beta 
stability). 

When this procedure was first contemplated it was thought that the right 
hand sides of equations (12) and (13) would have to be multiplied by some 
constant factor of order unity in order to correct for the errors in orbit spacings 
introduced by a calculation in which one adds particles rather than excites 
existing particles. The necessity for such a multiplicative constant would 
have been demonstrated if the resulting level-spacing formula did not have 
the correct dependence on A. However, the excellent fit which has been ob- 
tained shows that no such constant is needed. 

Orbit separation energies d(Z) and d(N) were calculated close to the bottom 
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of the valley of beta stability for 1 < Z < 139 and 1 < N < 239. These 
calculations required some further extrapolations of the writer’s shell correc- 
tions (Cameron 19576). 

It is necessary to average the orbit separation energies in the vicinity of the 
Fermi level. The average must be weighted so that large weight is given to 
orbits close to the Fermi level and small weight to those far from the Fermi 
level. The Fermi level is taken to coincide with the highest filled proton and 
neutron orbits. Following Newton (1956), we choose weights according to 
the formula: 


(14) w(E) = In(lt+exp—E/t), 


where E is the energy difference between the orbit in question and the Fermi 
level. It was found necessary to give equal total weights to the bound and to 
the unbound orbits. The average spacings are thus given by the expressions: 


(15) d(Z) = L d(Z) wi(Z) Ze y w(Z0, 
(16) d,(Z) = Ye d(Z.) wi(Z,) y m(Z0, 
(17) - any wixy / > 2 wi(V»), 
(18) d,(N) = - d(N,) w,(N;) =: w,(N), 
where 
. she \ 

(19) wi(Z,) = In| 1+exp—> 2 UZ)4 ’ 
(20) w(Z) =In2=w,(Z+1), 


zi-1 
(21) w,(Z:) in| 1-+exp—24 1 » d(Z; it, 
and similar expressions for w,(.V;), w,(.V), w,(.V;), and w,(V+1). 
The average orbit densities are then 


Gz = 2/[d,(Z)+d,(Z)], 
Gy = 2/[d,(N)+d,(N)]. 


In equations (15) to (18) Z,, Z,, Ni, and N, are limits at which the sums 
are cut off. In this work it has been customary to cut off the sums for w < 0.01. 

It may be seen from equation (4) that the nuclear temperature depends 
upon the total orbit density G, which in turn depends upon the nuclear tem- 
perature through the weights w. It is therefore necessary to determine G and 
t simultaneously by an iterative procedure. One assumes a value of G, cal- 
culates ¢, calculates G, recalculates ¢t, and so on until successive values of ¢ 
agree sufficiently closely (it has here been required that two successive approxi- 
mations to ¢ agree to within 1%). 


ll 


ll 
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4. THE FIT TO OBSERVED LEVEL SPACINGS 


The most accurate experimental information about nuclear level spacings is 
that derived from low energy neutron interactions with nuclei. The number 
of partial waves is then very limited (usually only s-waves), and the levels 
which are formed have only two spin states. An attempt was made to fit 
observed level spacings in 82 nuclei. Most of the observed spacings were taken 
from a paper on radiation widths by Stolovy and Harvey (1957). The remain- 
ing spacings were for nuclei used by Newton (1956) which were not included 
in the paper by Stolovy and Harvey. The writer has re-estimated the spacings 
in these latter nuclei through an inspection of BNL-325 (Hughes and Harvey 
1955) and its supplement (Hughes and Schwartz 1957). Levels were counted 
only up to energies at which it appeared likely that nearly all levels had been 
resolved. In cases where levels above 50 kev were included in the count, correc- 
tions were made for p-wave interactions (these were important for light 
nuclei and also near Pb’). These estimates are given in Table II. 


TABLE II 
SPACINGS OF LEVELS FORMED BY S-WAVE NEUTRON 
INTERACTIONS WITH NUCLEI AS ESTIMATED BY 
THE WRITER FROM BNL-325 
(The Z and A refer to the compound nucleus 
formed by neutron bombardment) 


a A Dovs 
9 20 110000 ev 

il 24 150000 
13 28 60000 
15 32 50000 
17 36 10000 
21 46 2000 
23 52 2000 
25 56 2000 
27 60 3000 
43 100 14 
50 113 150 
50 117 150 
50 119 300 
50 120 150 
59 142 110 
64 158 2 
65 160 But 
67 166 7.5 
72 180 8 
81 204 1000 
82 207 25000 
82 209 400000 
83 210 40000 
92 235 17 
92 237 17 
95 242 0.43 





When attempts were made to fit equation (7) to these observed level spacings 
by the procedure which has been outlined, the results were for the most part 
very encouraging. However, in certain regions the fit was poor. These regions 
were the transition regions near heavy closed shells in which the addition of a 
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small number of nucleons is sufficient to cause spherical nuclei to become 
strongly deformed, or vice versa. In these regions the calculated values of 
d(Z) and d(N) were negative. It was thus evident that any close correspon- 
dence between the second differences of the masses and the spacings of nucleon 
orbits had been destroyed by collective effects in the transition regions. 
Accordingly, typical values were assigned to d(Z) and d(.NV) in the transition 
regions, usually about two-thirds of the values that d(Z) and d(.V) would 
have if the shell corrections had been zero. When these substitutions had 
been made, a good fit was also obtained in the transition regions. 

It was found also that in certain regions away from shell edges the d(Z) 
and d(NV) fluctuated randomly between moderate and very small values. It 
was known from the work done in fitting masses (Cameron 19575) that these 
fluctuations had no physical significance. If the fluctuations were allowed to 
remain, they would cause large fluctuations in level densities calculated for 
very low excitation energies. Therefore the d(Z) and d(N) were smoothed 
locally in such regions. This had very little effect on the fits at the excitation 
energies corresponding to neutron binding energies. A table of d(Z) and d(V) 
is given as Appendix I. 


10 
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Fic. 1. Ratio of observed to calculated level spacings for 82 nuclei, plotted as a function 
of mass number. 





Figure 1 shows the ratio of the observed to the calculated level spacings as 
a function of mass number. It may be seen that every region of mass number 
has been fitted about as well as any other. The normalization to the points in 
Fig. 1 gives the final level spacing formula: 


z 


"ee 
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4 ae 3.26 10° de $\.1" r 
(24) D(UZNJ) = “(2J+1) | GsGuG A = a 2)e] exp—"3 Gt. 


Factors of error for the points in Fig. 1 are defined as follows: 


log Dors. 


(25) Log (error factor) = = 
Daeate 








The average error factor is 1.53. The root mean square error factor is 1.74. 
This latter value should be compared with the root mean square error factor 
of 3 obtained by Newton (1956). It is evident that a further improvement in 
the fit to level densities has been achieved by the procedure outlined in this 
paper. A not insignificant part of the error factor of 1.74 may be due to un- 
certainties in observed level spacings. Many of the values used in the fit are 
uncertain by at least a factor two. 

Figure 2 shows the ratio of the observed to the calculated level spacings 
plotted against the spin of the target nucleus which is bombarded by neutrons. 
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Fic. 2. Ratio of observed to calculated level spacings for 82 nuclei, bombarded by neutrons, 
plotted as a function of the target spin. 


The spins of the compound nuclei will differ from that of the target nucleus 
by +3. It may be seen that the fit is about equally good for all target spins, 
thus showing that level spacings vary inversely as (2J+1) to a good approxi- 
mation. 

It should be noted that the level spacing formula has not been tested over 
a very wide range of excitation energy. 
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5. APPROXIMATIONS TO THE ORBIT DENSITIES 


Although the method outlined in this paper for determining Gz(t) and 
Gy(t) has proved quite successful, it is nevertheless a very clumsy and lengthy 
method. The calculation of one level density takes between 1 and 10 
minutes on the Chalk River Datatron, the longer time being characteristic 
of nuclei near Pb®*. Under these circumstances it is not possible to use this 
procedure as a tool in the solution of other problems which may require the 
calculation of many level densities. Therefore the regular procedure has 
been used to calculate Gz(t) and Gy(t) for a range of values of ¢t (generally 
corresponding to a range of excitation energy of zero to 100 Mev). These 
values have then been fitted by the general approximation formula: 


ae - ao(k) +ai(k)t+a2(k) yt" ‘baa(k)t- task) 
(26) G,(t) = Se 

1+¢+¢°+0°4+0° 
where k stands for Z or \V. The five coefficients ao(k),... , as(k) were fitted 
by a least squares procedure to 15 values of G,(t). The resulting approximation 
usually fitted these values to within 1%. 

These fits were made for 10 < Z < 109 and 10 < V 189. The resulting 
coefficients are given in Tables III and IV. It may be Bie that some of the 
coefficients are negative. However, it has been verified that no set of co- 
efficients gives a negative value of G,(t) within the useful range of ¢. 


TABLE III 
COEFFICIENTS IN THE APPROXIMATION FORMULA (26) FOR G2(t). 
(Gz(t) is in units of Mev~!) 





iz ao a a2 a3 a, 

10 2.54090 —6.88500 7.06180 —2.66610 0.58669 
11 2.05760 —5.25000 5.36770 — 1.91450 0.48746 
12 1.82440 —3.50200 3.93670 — 1.08140 0.49654 
13 1.75700 — 2.00900 2.97640 —0.49370 0.50516 
14 1.23250 —1.19700 2.27190 —0. 15630 0.54597 
15 1.07400 — 1.43700 2.35670 0.11140 0.55716 
16 1.07760 —0.66700 1.64230 0.59240 0.56247 
17 0.97450 0.22700 1.00450 1.00930 0.57339 
18 1.02940 0.77900 0.57730 1.37540 0.59761 
19 1.38120 0.12700 0.91440 1.46700 0.64208 
20 0.82010 0.53200 0.25000 1.72890 0.68103 
21 0.73380 1. 12900 —0.35370 1.99080 0.68045 
22 1.78550 —0.21700 1.38770 1.61690 0.77693 
23 1.30070 1.59400 0.17310 1.96610 0.79870 
24 0.94390 2.95310 —1.14950 2.41430 0.78967 
25 0.77270 3.52960 —1.62140 2.39850 0.84887 
26 —0.04100 6.96580 — 5.70350 3.84030 0.71714 
27 0.13470 5.62540 —4.67970 3.19150 0.84693 
28 —0.01630 4.66900 —4.14950 2.84330 0.93519 
29 0.57050 1.29030 —0. 13600 1.06950 1.21312 
30 0.97740 0.77100 1.03660 0.64720 1.28997 
31 1.33880 0.94890 1.38730 0.86970 1.29596 
32 1.34420 1.49790 1.43370 1.15010 1.29318 
33 1.79310 0.87790 2.57380 0.97090 1.36046 
34 2.30220 0.42990 3.05330 1.12470 1.38150 
35 2.24120 0.50170 3.28720 1.19380 1.42727 
36 1.83260 1.29290 2.40790 1.67410 1. 40286 
or 1.73610 1.22640 2.55910 1.65670 1.45010 
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TABLE III (continued) 
Z a a az a3 as 
38 1.80250 0.98900 2.80540 1.57840 1.50325 
39 1.72880 1.66180 1.85380 2.06680 1.46949 
40 1.76890 1.64160 1.91960 1.98760 1.51824 
41 1.80940 1.66110 2.03090 1.88230 1.57321 
42 1.89230 1.34410 2.70580 1.37500 1.70617 
43 1.89830 1.51060 2.58460 1.32130 1.75270 
44 1.88300 1.77150 2.34820 1.24600 1.80983 
45 1.93330 1.83000 2.31840 1.08200 1.88211 
46 2.04170 1.64150 2.38290 0.80770 1.97991 
47 2.26540 1.03400 2.50520 0.60540 2.05252 
48 2.60780 —(. 43680 3.28500 0.21940 2.13629 
49 2.41160 — 1.55210 3.97610 —(0.43440 2.29490 
50 0.73028 1.61310 —0.79680 1.11850 2.13641 
51 0.62134 1.76260 —1.21720 1.77750 1.93396 
52 2.20146 —0.70960 2.41220 0.51580 2.15497 
53 2.20045 0.97640 1.46610 1.24090 2.08108 
54 2.02370 1.96340 1.53190 1.48450 2.08343 
55 1.81430 2.68210 1.98800 1.44040 2.13392 
56 2.15520 2.44970 2.96980 1.29650 2.18503 
57 2.78490 1.95430 3.96360 1.30290 2.21492 
58 2.88500 3.45450 2.32200 2.53780 2.05018 
59 3.54270 2.91280 2.95300 2.68350 2.05714 
60 4.94570 —0.47280 6.80350 1.23850 2.32293 
61 5.18540 —1.17040 7.72260 1.00950 2.41191 
62 3.87360 2.00980 4.18400 2.62370 2.20457 
63 2.84670 3.89310 2.06030 3.55600 2.09311 
64 2.65694 3.18680 2.85390 3.09850 2.20282 
65 2.66221 2.17600 3.83120 2.57700 2.32124 
66 2.48543 2.14120 3.53000 2.72270 2.31589 
67 2.31180 2.36530 2.95080 2.98440 2.29285 
68 2.37716 2.09170 3.19930 2.77600 2.35510 
69 2.60136 1.56730 3.66300 2.48930 2.42417 
70 2.49254 1.76560 3.74210 2.26540 2.49443 
71 2.35504 2.15640 3.69930 2.04360 2.57392 
72 2.37008 2.51860 3.65590 1.81040 2.65373 
73 2.52609 3.16450 2.79100 2.20410 2.59819 
74 3.02744 3.07870 2.69460 2.17460 2.63103 
75 3.40231 3.45160 1.96100 2.41220 2.62283 
76 3.61399 4.59460 —0.78900 4.08380 2.29499 
77 4.28974 2.91840 0.67410 3.14380 2.48983 
78 5.08195 —0.52930 4.57210 0.66070 2.99123 
79 5.49290 — 4.31320 9.15820 —2.46220 3.63057 
80 4.76203 — 5.54530 10.60750 —4. 23370 4.06301 
81 1.56832 2.19680 —0.98880 1.08570 3.15736 
82 0.63248 0.91420 —0.02280 1.09140 3.20923 
83 0.85952 0.70290 1.71860 — 1.25420 3.58958 
84 3.88686 —6.40780 12.39100 — 6.02830 4.42853 
85 5.62571 —7.94070 14.63030 — 5.86640 4.36500 
86 5.60035 —2.93880 7.46750 —0.53400 3.30689 
87 5.37494 1.21960 2.01890 3.42210 2.55312 
88 5.15555 3.89580 —0.87940 5.60320 2.17708 
89 4.57913 7.29850 — 5.00680 8.41830 1.66300 
90 4.94414 6.02950 —2.23660 7.06140 1.99504 
91 5.84768 1.79290 4.48440 3.45800 2.76167 
92 4.96580 3.80280 2.68080 4.35710 2.68089 
93 3.99003 5.79620 0.75160 5.43370 2.52601 
94 3.52739 5.76460 1.31870 5.09940 2.62699 
95 3.36583 5.06740 2.51560 4.48820 2.77705 
96 3.46175 4.07180 3.65780 4.03910 2.87281 
97 3.45633 3.86930 3.66190 4.05240 2.92873 
98 3.50179 3.72480 3.44800 4.25230 2.91469 
99 3.55396 3.56250 3.37660 4.31570 2.95171 
100 3.62648 3.31610 3.52580 4.19890 3.02965 
101 3.72031 2.86080 3.95910 3.84590 3.15330 
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TABLE III (concluded) 






















































Z a a a2 a3 a4 
102 3.76769 2.40420 4.40690 3.50280 3.27771 
103 3.67740 1.88570 5.31650 2.84130 3.46659 
104 3.27821 2.14940 5.42540 2.51880 3.59350 
105 2.85502 2.93760 4.49230 3.06580 3.48965 
106 2.89336 2.54480 5.18630 2.70070 3.55147 
107 3.15908 2.23410 5.43010 2.73370 3.54030 
108 3.31585 2.58470 4.91310 3.18420 3.44506 
109 3.35546 3.46910 3.62430 4.20510 3.18402 

TABLE IV 
COEFFICIENTS IN THE APPROXIMATION FORMULA (26) FOR Gy(t) 
(Gy(t) is in units of Mev~) 

N do a a2 a3 as 
10 —2.01000 6.37300 —3.46800 1.99730 0.37526 
11 —2.27300 7.19400 —3.75600 2.11950 0.43772 
12 —1.72100 5.93700 —3.11900 1.96670 0.50522 
13 0.27400 2.81200 — 1.42300 1.55070 0.59158 
14 — 4.62700 12.49300 —7.28600 2.68610 0.58200 
15 —2.24700 7.28500 —4.71400 2.12100 0.66448 
16 0.75300 0.89600 —0.20200 0.96220 0.80110 
17 —0.01900 3.65100 — 1.76400 1.48780 0.80672 
18 0.88100 1.63700 —0.16700 1.23570 0.88170 
19 1.15900 0.55300 0.90100 1.04210 0.95641 
20 0.70400 0.65700 1.04100 1.13680 0.99477 
21 0.75900 1.42600 0.76300 1.25750 1.01410 
22 1.38900 2.97800 0.11500 1.77780 1.02907 
23 2.35200 1.91200 1.19100 1.66390 1. 12262 
24 2.55600 0.97700 2.20400 1.22460 1.25306 
25 1.13800 4.08400 —0.61800 1.92310 1.25547 
26 0.65300 4.85000 — 1.82800 2.17130 1.29349 
a 0.33800 5.13400 —3.16000 2.47540 1.32796 
28 —0.52290 5.88800 —4.81400 2.92910 1.33400 
29 —0.62330 5.77100 — 4.76000 3.07560 1.30398 
30 0.85800 3.07600 — 1.68500 2.08310 1.44368 
31 1.74660 1.68900 0.48200 1.61130 1.55552 
32 1.84010 1.65500 1.38800 1.53770 1.62671 
33 1.88110 1.68200 1.96400 1.61890 1.66548 
34 1.93890 1.74200 2.41200 1.73940 1.70445 
35 2.23470 1.42300 3.12600 1.74670 1.75216 
36 2.78700 0.31000 4.59800 1.42840 1.85296 
37 3.28660 —0.79500 6.15600 1.02330 1.96976 
38 3.11860 0.09400 5.49900 1.43550 1.98554 
39 2.89480 0.87300 5.15900 1.62700 2.02530 
40 2.54000 2.43100 3.73100 2.20970 2.02364 
41 2.55200 2.85700 3.37900 2.36250 2.06658 
42 2.69850 3.63100 2.42200 2.75370 2.07943 
43 3.64080 2.55300 2.99100 2.53680 2.17423 
44 4.10290 2.08000 3.13400 2.34950 2.27746 
45 3.57860 2.68800 2.55000 2.14260 2.40308 
46 2.89790 4.15700 0.39100 2.79400 2.38510 
47 2.88360 3.39100 0.49000 2.44650 2.51017 
48 2.97650 2.00700 0.69300 2.11510 2.62291 
49 2.20430 2.26500 —0.56300 2.25900 2.68696 
50 0.36630 4.99800 —4.91500 3.70550 2.57086 
51 0.59230 2.56900 — 1.84300 2.96190 2.56207 
52 1.78910 1.04800 0.45000 1.95430 2.74485 
53 1.71180 2.86200 —0.46400 2.46000 2.70429 
54 1.65080 4.64700 — 1.35500 3.01100 2.64281 
55 2.70740 3.86600 —0. 13200 3.09680 2.62880 
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TABLE IV (continued) 


1051 





N ao ay, a2 a3 a 
56 3.92370 2.35000 2.25100 2.70700 2.72411 
57 4.11380 2.75400 2.77700 3.00040 2.74027 
58 3.77330 3.94800 2.58000 3.42240 2.74288 
59 3.32810 6.17300 0.45000 4.84870 2.57430 
60 3.31030 6.45600 0.86400 5.00610 2.59945 
61 3.04530 7.99200 —0.52200 6.02320 2.49465 
62 3.23470 8.12500 —0.49300 6.25530 2.48990 
63 4.50300 5.44300 2.31900 5.51490 2.65379 
64 5.43750 4.07200 3.70100 5.34360 2.72528 
65 4.62960 7.05800 0.87300 6.69820 2.59748 
66 4.10510 8.61300 —0.25100 7.04340 2.63621 
67 5.06110 6.23400 2.19100 6.19140 2.80129 
68 5.94120 3.96300 4.56000 5.24150 3.01265 
69 4.70240 7.40300 0.98400 6.66220 2.87188 
70 3.26820 11.79600 —3.81100 8.66880 2.63130 
71 3.45980 11.59100 —4.14800 8.84410 2.64101 
72 4.82300 7.67400 —0.15700 7.03040 2.96799 
73 5.65400 5.97700 0.78300 6.54580 3.08789 
74 4.87440 8.97600 —3.34100 8.23430 2.88701 
75 4.31010 10.22700 — 4.85000 8.43420 2.94242 
76 5.55520 6.24490 — 1.69400 7.04660 3.17342 
77 8.20490 —2.35650 6.45500 3.33950 3.77604 
78 9.78170 —8.31620 11.34000 0.92570 4.18195 
79 8.89640 —8.09010 9.69600 1.08890 4.22318 
80 5.55150 —1.41700 0.63770 4.40940 3.78019 
81 1.76200 5.69660 —8.38780 7.98270 3.30502 
82 0.79360 2.24950 —3.12110 5.85450 3.67937 
83 0.90040 1.91210 —2.22470 4.90410 3.77064 
84 2.21280 3.35830 —3.81270 5.35250 3.69503 
85 4.48590 1.27410 —0.34820 4.15720 3.89296 
86 7.58190 —3.79360 6.37050 1.79000 4.27251 
87 10.53330 —8.99900 13. 45800 —0.67680 4.68776 
88 11.96770 —10.24560 15. 88800 —1.13360 4.79335 
89 12.31040 —8.40750 14.01100 0.61690 4.49846 
90 12.22800 —7.04840 13 .50300 1.21760 4.44878 
91 11.58070 —5.09040 12. 44700 1.89840 4.39150 
92 9.86250 —(. 43420 7.93500 4.19100 4.05640 
93 7.81890 4.43690 3.46800 6.20680 3.78750 
94 6.14920 7.43030 1.30400 7.02070 3.71880 
95 4.93470 9.37510 —0.23400 7.71410 3.63770 
96 4.29580 9.16020 1.07300 6.88460 3.82260 
97 4.14940 7.94990 3.20900 5.78790 4.03630 
98 4.27180 6.03760 6.30000 4.02360 4.41810 
99 4.28080 5.44040 7.13400 3.65850 4.50110 
100 4.38820 5.09990 7.36900 3.65780 4.52350 
101 4.70600 4.44620 7.80600 3.58260 4.54500 
102 5.11340 3.38080 9.29100 2.70850 4.75280 
103 5.32590 3.21600 9.23500 2.79070 4.75830 
104 5.38920 3.38840 9.20400 2.63050 4.84500 
105 5.47650 3.28700 9.55200 2.25120 4.96700 
106 5.51350 3.34390 9.76200 1.92910 5.07630 
107 5.27980 4.30570 8.80300 2.14540 5.09780 
108 5.28230 4.47900 8.82400 1.82290 5.21230 
109 5.49430 3.85470 9.93700 0.91550 5.44680 
110 5.74000 3.27800 10.68400 0.25920 5.61740 
111 5.85290 3.05900 11.13300 —0.43680 5.82450 
112 5.75270 3.44930 10.89600 —0.80600 5.99010 
113 5.72910 3.71210 10. 43500 —0.94970 6.08100 
114 5.50360 5.33220 6.89920 1.27590 5.57640 
115 5.40300 5.92930 5.61720 1.65350 5.56110 
116 5.42110 5.93650 5.03990 1.62540 5.62550 
117 5.49550 5.87270 4.07270 1.96120 5.59070 
118 5.59090 5.57640 3.27710 2.16750 5.59820 
119 5.62980 5.36060 2.01680 2.70800 5.52860 








CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 






TABLE IV (continued) 

































N ao a a2 a3 a4 
120 5.66620 4.73380 1.00260 3.27500 5.44026 
121 5.67060 3.50780 0.78880 3.31940 5.45726 
122 5.18893 3.31780 —0.98520 4.30800 5.28631 
123 4.44115 3.12600 —2.67410 5.34260 5.09107 
124 3.94805 0.46200 —0.53990 4.34570 5.27667 
125 1.66941 3.65770 —5.70440 7.07060 4.80341 
126 0.80912 0.58190 0.84720 4.24580 5.14117 
127 1.23839 0.82540 2.48570 —2.45000 6.60541 
128 3.23842 —1.91680 8.80140 — 4.66090 7.07604 
129 7.23749 —13.92430 28 .01920 —14.47400 9.02634 
130 11.37358 —26.88090 48. 45080 —25. 29430 11.25946 
131 13 .62251 —31.13200 55.11310 — 27 .82330 11.77850 
132 14.76010 —30.84020 54.71950 — 26. 28380 11.48710 
133 14.57480 — 25 .36240 46.77530 —20.65240 10.42810 
134 14. 15950 —19.75880 38. 44600 —14.77130 9.28260 
135 13 .23160 —13.01450 28 .77300 —8.30630 8.03760 
136 11.99790 —5. 85650 18. 35800 — 1.44620 6.67990 
137 11.18760 — 1.62000 12.58600 2.51450 5.91660 
138 10.87090 —0.72080 12. 16300 3.06420 5.86780 
139 10.52770 —0.89560 13.83700 2.18780 6.12870 
140 9.41760 1.44650 11.83300 3.35750 5.97910 
141 6.33490 12.06220 —2.49200 11.77350 4.36050 
142 5.37460 14.86770 —5.20500 13.17730 4.15420 
143 6.92580 5.94960 9.31000 4.42380 6.02650 

‘144 6.41020 6.33040 9.87500 3.99240 6.19510 
145 4.68490 13 .02210 0.54500 9.56500 5.11670 
146 5.05440 11.82440 1.90300 9.17060 5.19250 
147 7.93780 —0.00180 19.08100 —0.65990 7.21970 
148 11.11410 —13.60990 39 45300 —12.59010 9.73350 
149 13 .45580 —25.31500 57 .93300 —23 .82260 12. 16380 
150 11. 24360 —19.48420 5100000 —2). 10400 11.53720 
151 6.87400 —5.01680 31.93900 —9.37830 9.44750 
152 2.51800 9.04440 14.40000 0.05660 7.68050 
153 2.85130 7.70580 15.49100 —0.50290 7.74770 
154 7.65620 —7.34620 36. 24600 —12.40810 10. 25200 
155 9.76270 —10.98020 39.81700 —13.79610 10.53150 
156 9.86980 —6.55190 31.67600 —8.35830 9.41860 
157 9.72860 —2.26000 24. 28900 —3.50890 8.43560 
158 9.32240 2.32390 16.87700 1.15380 7.49220 
159 9.17610 4.98160 12.90700 3.47850 7.08120 
160 9.08990 6.68370 11.17000 4.09000 7.09430 
161 8.76630 9.46520 7.17000 6.28290 6.72320 
162 8.63630 11.34140 4.46300 7.66880 6.52450 
163 8.63930 12.75230 2.16600 8.90000 6.35900 f 
164 8.63570 14.39110 —0.83600 10.52040 6.08610 f 
165 8.73020 15.52130 —2.68200 11.43760 5.96600 
166 9.29130 14.77110 —2.35400 11.12320 6.11070 
167 10.30140 12.39950 —0.09200 9.76800 6.44140 
168 11.06960 11.16790 0.06400 9.79280 6.46290 
169 11.46850 11.61530 —2.64700 11.72790 6.06890 
170 12.20170 9.83190 —2.29900 12.08730 5.82800 
171 13 .02650 7.11960 0.00800 10.72030 6.13250 
172 13.27010 6.48600 —1.00100 11.33450 5.99250 
173 13 .26560 6.07640 —2.18000 12.02210 5.82560 
174 13.17270 5.42440 —2.90500 12.28090 5.80400 
175 13.21756 3.20470 —1.11600 10.83180 6.12770 
176 13 .33235 —0.38930 3.26500 7.53960 6.90970 
177 13 .47239 —5.17700 9.82000 2.59790 8.06580 
178 13.61737 —11.34520 19.27830 — 4.59480 9.78480 
179 13.78681 —19.06110 32 .00070 —14.35500 12.14410 
180 13 .32149 —24.81640 41.72400 —22.31430 14. 13590 
181 11.47239 —24.29050 40 .76040 — 23 .23280 14.49410 
182 6.65768 —8.54540 12.91000 —5.69670 10.56670 


1.96431 .31070 —11.70070 9.67710 .08568 
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TABLE IV (concluded) 











N ao ay az a3 a 
184 0.61563 2.40930 —6.54590 10. 29300 6.72776 
185 1.01590 0.87910 2.18890 —6.38400 11.07745 
186 2.28950 0.29780 5.39480 —6.51590 11.12806 
187 3.98508 —2.98150 13.31490 — 10.26400 12.07910 
188 5.98064 —8.22120 24.19440 15.98570 13 .45220 
189 9.85338 — 24. 06020 53 .69850 _ ~ 35 02800 17 .99470 





In using these approximations it is still necessary to determine G and ¢ 
simultaneously by an iterative procedure. It is recommended that an initial 
value of G(t) be chosen as 


(27) G(t) initial = a4(Z) +a,4(.V) 


Near closed shells it is often better to choose an initial value about one half 
of the amount given by equation (27). 


The writer is indebted to Dr. T. D. Newton for several stimulating dis- 
cussions. 
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APPENDIX I 
A table of the nucleon orbit spacings d(Z) and d(N) calculated from second differences of 
masses. [d(Z)], and [d(NV)]|, are reference spac ings obtained from the mass formula from which 
the shell correction terms were omitted in addition to the Coulomb, Coulomb exchange, and 
pairing enneny terms. {d(Z)], and [d(V)], are the spacings calculated as described in the text, 
and [d(Z)|a and [d(N)]Ja are the spac ings, smoothed and in ple aces adjusted as desc ribed in the 
text, which were used in the final calc ae of level apa ings. All apne ings are in Mev. 











“(dU N Ya 





Zor N Ze “(d(Z) Ne - (d(Z) De la Nk ta N Ne 
0 40.000 40.000 
1 30.000 30.000 
2 20.000 20.251 20.000 
3 8.935 8.935 8.961 8.961 
4 5.465 5.465 5.713 5.713 
5 4.247 4.247 4.035 4.035 
6 3.415 3.415 3.247 3.247 
7 2.875 2.875 2.742 2.742 
8 2.494 2.494 2.390 2.390 
9 2.212 2.212 2.126 2.126 
10 1.992 1.992 1.920 1.920 
11 1.865 10.145 10.145 1.754 0.854 0.854 
12 1.715 0.765 0.765 1.554 1.544 1.544 
13 1.585 1.555 1.555 1.447 1. 407 1.407 
14 1.478 0.778 0.778 1.356 0.946 0.946 
15 1.411 2.351 2.351 1.275 1.655 1.655 
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Zor N [d(Z)], (d(Z)}u [d(Z Ja (d(N)]- [d(N Ju [d(N)}a 
16 1.327 1.587 1.587 1.205 1.955 1.955 
17 1.253 1.433 1.433 1.110 1.150 1.150 
18 1.187 1.157 1.157 1.057 1.067 1.067 
19 1.146 0.996 0.996 1.010 1.030 1.030 
20 1.090 0.610 0.610 0.967 1.037 1.037 
21 1.041 2.031 2.031 0.928 1.648 1.648 
22 1.010 0.440 0.440 0.869 0.649 0.489 
23 0.968 0.728 0.728 0.840 0.010 0.245 
24 0.930 0.790 0.790 0.809 0.309 0.310 
25 0.904 0.794 0.794 0.782 0.612 0.536 
26 0.870 0.800 0.800 0.741 0.651 0.651 
27 0.849 0.809 0.809 0.718 0.658 0.658 
28 0.821 0.991 0.991 0.698 0.718 0.718 
29 0.794 1.564 1.564 0.679 1.749 1.749 
30 0.775 1.365 1.365 0.646 1.066 1.066 
3l 0.750 0.790 0.790 0.630 0.530 0.530 
32 0.735 0.765 0.765 0.615 0.575 0.575 
33 0.714 0.704 0.704 0.587 0.557 0.557 
34 0.698 0.518 0.518 0.573 0.553 0.553 
35 0.678 0.448 0.448 0.560 0.510 0.510 
36 0.664 0.554 0.554 0.550 0.450 0.450 
37 0.648 0.598 0.598 0.526 0.406 0.406 
38 0.636 0.606 0.606 0.516 0.356 0.356 
39 0.620 0.580 0.580 0.507 0.397 0.397 
40 0.607 0.577 0.577 0.487 0.397 0.397 
41 0.595 0.565 0.565 0.479 0.389 0.389 
42 0.584 0.554 0.554 0.469 0.399 0.399 
43 0.571 0.541 0.541 0.455 0.345 0.345 
44 0.561 0.541 0.541 0.446 0.166 0.166 
45 0.550 0.530 0.530 0.439 0.269 0.269 
46 0.539 0.519 0.519 0.425 0.325 0.325 
47 0.530 0.500 0.500 0.418 0.358 0.358 
48 0.519 0.459 0.459 0.411 0.321 0.321 
49 0.513 0.403 0.403 0.400 0.320 0.320 
50 0.503 0.233 0.233 0.393 0.323 0.323 
51 0.495 2.145 2.145 0.388 1.578 1.578 
52 0.488 0.498 0.498 0.376 0.666 0.666 
53 0.478 0.488 0.488 0.371 0.501 0.501 
54 0.472 0.482 0.482 0.366 0.596 0.576 
55 0.465 0.515 0.515 0.356 0.516 0.440 
56 0.456 0.536 0.536 0.352 0.132 0.259 
57 0.449 0.379 0.379 0.347 0.257 0.226 
58 0.444 0.354 0.354 0.339 0.259 0.259 
59 0.438 0.338 0.338 0.334 0.274 0.274 
60 0.430 0.320 0.252 0.329 0.269 0.269 
61 0.426 0.046 0.171 0.323 0.263 0.263 
62 0.420 0.270 0.230 0.317 0.287 0.287 
63 0.413 0.333 0.317 0.311 0.261 0.261 
64 0.409 0.379 0.379 0.308 0.148 0.148 
65 0.402 0.362 0.362 0.304 0.154 0.154 
66 0.397 0.397 0.397 0.297 0.247 0.247 
67 0.393 0.413 0.413 0.293 0.233 0.233 
68 0.387 0.447 0.447 0.290 0.140 0.140 
69 0.383 0.403 0.403 0.283 0.123 0.123 
70 0.378 0.388 0.388 0.282 0.252 0.252 
71 0.374 0.434 0.434 0.274 0.264 0.264 
72 0.369 0.419 0.419 0.272 0.232 0.217 
73 0.366 0.416 0.416 0.271 0.171 0.157 
74 0.362 0.362 0.362 0.264 0.054 0.145 
75 0.358 0.298 0.298 0.262 0.302 0.226 
76 0.353 0.293 0.293 0.256 0.246 0.238 
77 0.348 0.258 0.258 0.256 0.156 0.150 
78 0.345 0.225 0.225 0.252 0.042 0.082 


79 0.343 0.103 0.200 0.246 0.086 0.077 





a 
& 
he 





Zor N 


(d(Z)], 


0 
0 
0 
0 
0 


0. 
o. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


0 
0 
0 
0 
0 
0 
0 


.338 
335 
331 

327 
.324 
320 
318 
316 
313 
309 

304 
302 
300 
298 
.293 
. 292 
.291 

. 286 
. 285 
. 282 
.281 

.276 
.273 
.273 
.270 
. 268 
. 266 
. 263 
. 261 

.259 
. 256 
. 255 
254 
252 
247 
. 249 
244 
.244 
. 242 
241 

. 239 
238 
235 
233 
. 232 
. 230 
. 229 
. 226 
224 
. 223 
221 
. 222 
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(d(Z)}u 


0. 
0. 


038 
075 


.051 
.997 
394 
.120 
.042 
.116 
.013 
.059 
. 234 
.172 
.170 
. 238 
. 253 
. 292 
. 281 
. 286 
. 285 


282 


. 281 
.276 
.273 
.273 
. 290 
.338 
. 366 
.333 
dll 
. 299 
. 296 
. 265 
. 254 
. 242 
.227 
.219 
. 204 
.194 
.172 
151 
.139 
. 162 
.065 
.083 
.132 
. 180 
.029 
. 826 
. 164 
143 
161 
112 


(d(Z)]o 


0.200 
0.200 
1.051 
1.997 
0.394 
0.200 
0.200 
0.200 
0.200 
0.200 
0.234 
0.172 
0.170 
0.238 
0.253 
0.292 
0.281 
0.286 
0.285 
0.282 
0.281 
0.276 
0.273 
0.273 
0.290 
0.338 
0.366 
0.333 
0.311 
0.299 
0.296 
0.265 
0.254 
0.242 
0.227 
0.219 
0.204 
0.194 
0.172 
0.151 
0.150 
0.150 
0.150 
0.150 
0.132 
0.180 
1.029 
1.826 
0.164 
0.143 
0.161 
0.112 
0.150 
0.150 
0.150 
0.150 
0.150 
0.150 
0.150 
0.150 


(d(N)]- 


0. 
0. 
0. 
0. 
0. 
0. 


244 
241 
238 
234 
231 
231 


0.225 
0.224 
0.221 


0. 


219 


0.214 
0.214 


0. 
0. 
0. 


212 
208 
207 


0.203 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


202 
198 
197 
197 
193 
192 
189 
187 
186 
185 
181 
179 
179 
177 
177 
172 
171 
170 
169 
166 
164 
165 
162 
160 
158 
159 
157 
154 
155 
152 
151 
149 
150 
148 
147 
145 
143 
143 
140 
139 
140 
139 
136 
137 
135 
133 
132 
131 


(d(.N lu 


0.094 
—0.009 
0.968 
1.524 
0.761 
0.261 
0.185 
0.034 
0.001 
—0.061 
0.024 
0.044 
0.122 
0.148 
0.147 
0.183 
0.202 
. 228 
.217 
227 
223 
. 222 
.199 
.197 
186 
195 
181 
189 
.199 
187 
. 187 
172 
181 
.180 
.179 
186 
184 
185 
.152 
.020 
.442 
—0.331 
0.147 
0.214 
0.245 
0.252 
0.861 
1.589 
0.060 
0.098 
0.087 
0.025 
—0.097 
0.003 
0.030 
0.059 
0.090 
0.079 
0.056 
0.047 
0.025 
0.083 
0.212 
0.141 


eccocoscocescococcoso 


cocoeceo 
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[d(N)a 


0.066 


0 
0 
1 
0 
0 
0 


0. 
0. 
.090 
.090 
0. 
.086 
.126 
147 


0 
0 


0 
0 
0 


0. 
. 202 
. 228 
.217 
227 
. 223 
. 222 
.199 
.197 
. 186 
.195 
.181 
.189 
.199 
. 187 
. 187 
.172 
.181 
. 180 
.179 
. 186 
. 184 


261 
724 
.524 
-590 
. 296 
.192 


090 
090 


090 


183 


185 
180 


.180 
.180 
.180 
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_42orN \d(Z)|- [d(Z)]u (d(Z)|a (d(N)}, [d(N)}u [d(N Ja 


. 132 0.132 0.132 
.129 0.169 0.169 
. 130 0.230 0.230 
.129 0.149 0.149 
. 126 0.116 0.116 
.127 0.077 0.077 
.125 0.035 0.035 
.123 . 133 0.133 
.123 . 163 0.163 
.122 492 0.492 
.122 112 0.112 
.120 . 100 0.100 
119 . 109 0.109 
.118 . 108 0.108 
119 . 109 0.109 
116 116 0.116 
118 . 108 0.108 
114 .114 0.114 
115 115 0.115 
114 .114 0.114 
.112 112 0.112 
113 ats 0.113 
.112 112 0.112 
Lil 101 0.101 
ahi .O91 0.091 
. 109 089 0.089 
.110 .090 0.090 
107 077 0.077 
0.108 .078 0.078 
107 077 0.077 
. 106 066 0.080 
. 106 .036 0.080 
. 104 306 0.080 
. 105 . 205 0.080 
. 103 .057 0.080 
101 .O19 0.080 
103 .003 0.080 
.102 .032 0.080 
.099 .029 0.080 
101 051 0.080 
.099 . 859 859 
.099 .979 1.979 
.097 O17 0.017 
.099 .O19 0.019 
.099 .019 0.019 
.095 .005 .050 
.096 —0.024 0.050 
.094 —0.116 0.050 
.094 —0.186 0.050 
.095 —0.165 0.050 
.094 0.014 0.050 
.094 0.064 0.064 
.0O91 0.071 0.071 
.093 0.103 . 103 
.091 0.111 Ill 
.093 0.113 0.113 
.089 0.109 . 109 
.089 0.109 . 109 
.090 —0.070 .090 
.088 0.088 .088 
.089 0.089 089 
.088 0.088 0.088 
.086 0.086 0.086 
.088 0.088 0.088 





144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 


cooceocoececo 


eeocococoeccocoo SEB SB OeeeMeAeSSeSeeSseeeeeeese = 
eoococoocecoeoceceooecoececooececececece|ce|cece|cecececscoceoceco 








ZorN 


208 
209 
210 
211 
212 
213 
214 
215 
216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
232 
233 
234 
235 
236 
237 
238 
239 


(d(Z)}, 
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\d(Z)]u 


(d(Z la 





(d(N)I, 


0.086 
0.085 
0.086 
0.085 
0.084 
0.084 
0.084 
0.084 
0.081 
0.082 
0.081 
0.083 
0.082 
0.079 
0.080 
0.081 
0.081 
0.078 
0.078 
0.079 
0.078 
0.079 
0.076 
0.077 


[d(.N lu 


0.086 
0.085 
0.086 
0.085 
0.084 
0.084 
0.084 
0.084 
0.081 
0.082 
0.081 
0.083 
0.082 
0.079 
0.080 
0.081 
0.081 
0.078 
0.078 
0.079 
0.078 
0.079 
0.076 
0.077 
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0.086 
0.085 
0.086 
0.085 
0.084 
0.084 
0.084 
0.084 
0.081 
0.082 
0.081 
0.083 
0.082 
0.079 
0.080 
0.081 
0.081 
0.078 
0.078 
0.079 
0.078 
0.079 
0.076 
0.077 
0.075 
0 075 
0.075 
0.075 
0.075 
0.075 
0.075 
0.075 











SCATTERING FROM A SMALL ANISOTROPIC ELLIPSOID! 


R. A. Hurp 


ABSTRACT 


Scattering of an electromagnetic wave by a small ellipsoid having tensor 
permeability and dielectric properties is dealt with by expanding the fields as 
power series in A~'. Consideration has been given to the first three terms of the 
expansion. 


INTRODUCTION 

The problem of solving Maxwell’s equations in an anisotropic medium 
for other than cylindrical coordinate systems is difficult, and to the author's 
knowledge, has received little attention in the literature. Jones (1945) has 
considered the case of an anisotropic ellipsoid in a static field, and some work 
has been done for ferrite-type anisotropy in spheres by Berk and Lengyel (1955), 
Ament and Passerini (1956), and Artman (1956, 1957). 

In this paper we shall study the diffraction of a plane electromagnetic wave 
by an ellipsoid which has permeability and dielectric properties described 


by tensors 
w= (uiy), € = (Ee), 1,j = 1,2,3. 


It is assumed that the dimensions of the ellipsoid are small compared to the 
wavelength inside the ellipsoid. The solution is obtained by expanding the 
various field components as power series in \~! in the manner used by Steven- 
son (1953 a, 6) for scattering from isotropic ellipsoids. 

The first three terms of the series are considered. Explicit determination 
of the fields is not, however, carried out for the most general permeability and 
dielectric tensors, since this requires the inversion of matrices of rather large 
order. When somewhat less general tensors are considered explicit evaluation 
of the fields becomes feasible. 

The case where the tensors have diagonal form is worked out, and in the 
paper following this, the coefficients are determined for a ferrite in which the 
steady magnetic field is applied along one of the principal axes of the ellipsoid. 


GENERAL 
We consider an ellipsoid, having the permeability and dielectric properties 
described above, with principal semi-axis Va; (7 = 1, 2,3), and whose equation 
in rectangular coordinates x; is* 


(1) =: x, /a; a 


1Manuscript received March 3, 1958. : 
Contribution from the Radio and Electrical Engineering Division, National Research 


Council, Ottawa, Canada. 
Issued as N.R.C. No. 4793. A 
*Summations and indices range over 1, 2,3 unless otherwise indicated. 
Can, J. Phys. Vol. 36 (1958) 
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We introduce the elliptical coordinates £; defined by the equations 


(2) Zz x; /(a:+é;) = 1, 
where & > —@3, ~—Gs3 - £5 os —d2, —a2 > &; = —@i, and a\>a2>43. Thus 
the surface of the ellipsoid is described by &, = 0. 

On the ellipsoid it is assumed that there is incident a plane electromagnetic 
wave 6, # whose rectangular components are 


(3) H, = mexp(—ik z a irtint) 
J 

(4) é,= Zolvexp ( —i8 b> apertint) , 
Jj 


with Z) = V 0/€0 and k = 27/\. The time dependence e‘”‘ will henceforth 
be omitted. Assuming that the magnetic field has unit amplitude, and since 
the fields & and #& are orthogonal to each other and to the direction of 
propagation, it follows that 


(5) Ll = Dme = Ln? = 1, 
(6) Llm, = Loman, = Dn; = 0, 
plus nine relations of the form 

(7) 1, = mon3—NMs3. 


Following Stevenson, the incident field is expanded as a power series in k. 


(8) HK, = > (ik)" HO? ; 


n=0 


Comparison of equations (3) and (8) yields 


(9) Ho) = (-1y'm( 5 wx) / (n!). 


We denote the external scattered fields by E, H, D, B, the internal fields by 
e, h, d, b, and expand each as power series in k 


(10) H = >° (ik)"H®, etc. 


n=0 


Substitution in Maxwell's equations yields the set of equations 


(11) VxXH® =VXxXE® =V-B® =V-D® = 0, 
( Vv xXH™ = De», 

(12) +VXE® = —cB®™-», 
| V-BM =V-D”M =0, forn>0O, 


where ¢ is the velocity of light in vacuum. The same set of equations holds for 
the internal field components. The M.K.S. unit system is employed. 
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Continuity of field components at £; = 0 leads to 


[ H+ LH?” —h” Jan = 0, 

| [B+ 103 —b” Joorm = 0, 
| [E+ 2™—-—e™ an = 0, 

| [D™ + 606 —d nom = 0. 


(13) 


We seek step-by-step solutions of equations (11) to (13) as far as m = 2 
subject to the condition that the scattered fields should vanish at infinity.* 


THE STATIC SOLUTION, MAGNETIC FIELD 


When n = 0, the equations which have to be solved are 


(14) VXH® =0, 

(15) V-H® = 0, 

(16) Vx h® =0, 

(17) V-b® = 0, 

with 

(18) oY” = x wish”, 


and # = m,; subject to boundary conditions (13). 
From the analogy with the isotropic static case, we assume solutions in 
the form of gradients of potential functions, thus 


(19) h© = Vg, 
(20) H® = Vo, 
with ¢ given by 
(21) ¢” = pi CiX iy 
and 
(22) a =  Awili(h),t 
where c,, A, are independent of the variables x,;, and 
= du 
(23) ne = f ota 
(24) R(u) = [(ai+u) (a2+u) (as+u)]?. 


*That the fields vanish at infinity for the first three terms was shown by Stevenson. A 
similar analysis shows that this is also true for anisotropic scatterers. 
tSee, for example, Whittaker and Watson, Chap. XXIII. 
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For continuity of the tangential field components it is sufficient to ensure 
that the potentials are continuous across £; = 0. This gives 


(25) zZ CX, = ba AJIet>d M ;X i, 


where the x; are evaluated at ¢; = 0, and J; = 7,(0). Equation (25) is satisfied 
independently of and &; if 
26) CC, = A J ;+m,. 


To evaluate the normal components, we have, for the normal component 
of any vector V with rectangular components V;;, 


(27) Ve, = U, > V; Ox ,/ OE; 


where U, isa scale factor whose value is of no consequence. From consideration 
of equation (2) it follows that 


(28) Ox ;,/ OE, = 3x;/(E;+a;). 

Hence, using equations (18), (21), (27), and (28), and putting &; = 0, 
(29) bey = 3U, > X Mijls/Ai, 

and for the incident field “ 

(30) HM) = 1U,>> mx ;/Q;. 

Now taking the £; gradient of & and putting £; = 0 we obtain 

(31) Hi? = Ui > (41;—A)x:A1/ai 


where 4 = (a,a2a;)~?. 
The continuity of normal components then yields 


(32) Do Mises = wo(Li—24) A ct pom,. 
J 
Eliminating A; between equations (26) and (32) we find 
(33) Dd [His—wo(1—-24/T,)8ij)e) = 2uom A/T, 
J 


where 6;; is Kronecker’s symbol, 
Oi: = I, 
5;, = 0, when 7+). 


Calling the determinant of the system A and the minors M4, the solutions 
of equation (33) are given by 


(34) Cj = Quod * = mM ;;/I;. 


The A, are then readily found from equation (26). 
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The zero order electric fields are written in the form 
E®=vs” May oO 
with 
o” = dA wili(ts), o” = » CX 4. 


The coefficients A ;, €; are found from the A,, ¢; by making the substitutions 
m, — Zoli, wij — €i3, aNd wo — €o in equations (33) and (26). 

In future a bar over a quantity will indicate the corresponding electric 
component. 


THE FIRST ORDER MAGNETIC FIELD 


From equation (12) the equations to be solved are 


(35) VXH® = cD, 
(36) V-H® = 0, 
(37) VX h® = cd, 
(38) V-b® = 0, 
with #7) = —m, x xm, and 
where the components of D and d™ are given by 
e — 
(39) a? = oo) Ale) + z. A jx; atest ’ 
j Xi 
(40) dy = Do eax. 


J 
We seek a solution to equations (35) and (36) of the form H® = G+4FO 
where G™ is a particular solution of VXG® = cD, V-G@ = 0, and 
F is the gradient of a harmonic function. The particular function G™ can 
be taken from Stevenson’s paper. Thus 


(41) GP = DY Gases (Es), 
J 

where Gii = 0, Gi; = —G351; 

and 

(42) Giz = e0cAg, (P)*. 


To find the complementary function we put F® = V6, where ® is a 
potential function of degree 2. The most general exterior harmonic function 
of this degree is 


(43) & = BX,T,(§:)+B'Xs'Ty' (b1) +xivars » ByJi(&1)/x4 


*The symbol (P) following an equation will henceforth be used to indicate that the cyclic 
permutations of the equation hold. 
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where B, B’, and B;, are coefficients to be determined, and 

R*(v)0(v), 

z x / (ait) I, 


1 


X, 
8(v) 


‘ du 
To) = Jen RO 
Ji(G) = oo ’ 


and », v’ are roots of the quadratic equation 


(44) Bv) = dF (ait+v)~* =10. 
To satisfy equations (37) and (38), the interior fields are taken to be of 
the form 


(45) hy = Do gas; 
J 
and therefore 
(46) a = » High jX ke 
<3 


Using equations (37) and (45) we obtain 
(47) £2—Z12 = cd, (P). 
Substitution of equation (46) into equation (38) gives 


(48) » High si = O. 
Matching of the normal fields at ; = 0 leads to 
- 2X N 4M 4X xX 5/Ajs+ » Gi] jx ix;/ai 
(49) +2BR*(v)(T+4/») 2) xa (atv) '+2B'R'Y) (7, +4/") 
x X xa; (atv) 
+2 1X 2X3 » Bix, [J :(8(0)—a;')—2Aa,] = po” Do *euttis8n/ Oe 
Continuity of the tangential components yields 
- d mm sx ix ;/(as+é)+ » G iy] sx x 5/(aité) 
(50) +2BT. Rv) D0 xc'(atv) (ait) '+2B'T, RG) 
XD xe actyy (até) 
+x1X2%3 x B,J {8(¢)— (ait+é) |] = » £13 /(a, +8), 
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where & stands for either of tangential variables £ or £3. Equating coefficients 
of x,x, in equation (49) gives 


(51) —nym,+2BR*(v)(T,+4 /v)/(ai+v)+2B'R°(v’) (T,+4 /v’) /(ai+v’) 


=) 
= Ko Zz Miskiir 
J 


and 


(52) —nm,/a;—nym;/a;+Gi;I;/a;:+G;1;/a; 


+B [J,(1/a;+1/a;)—2a,A*] = w'( 5 MikSkj fai+ >, antsi/4,) : 
C k 


In equation (52) the indices 7, 7, / are unequal. The tangential equations (50) 
yield 
(53) —nym,+2BT,R?(v)/(aitv)+2B’T,’'R?(v’)/(aitv’) —gii = 2d/ai, 


(54) —nymjt+Gy1 t+BJk = 83: (i, 7, R unequal). 


In equation (53) a new unknown ) is introduced because the variables 
x ;?/(a,+€) are not all independent, being connected by 


55) > xa, (a,+¢)* = 0. 


Relation (55) is established by differentiating the equation of the ellipsoid (1) 
with respect to either of the tangential variables and using equation (28). 

It is apparent that equation (54) and the electric version of equation (32) 
imply equation (47). Also, by summing equation (51) with respect to 7, using 
equation (44) and the orthogonal property of the m; and mj, we find 
di jHigsi = 0 and hence equation (48) is satisfied. We are thus reduced to 
15 equations in the 15 unknowns B,;, B, B’, X, and gi. The gi; can be 
eliminated from equations (51) to (54) leaving a set of six equations. The 
6 X 6 matrix of the set has all but three of its positions non-zero, so inversion 
of it will be, in general, a tiresome task. 

Before proceeding to the third term, we need expressions for the second 
term electric field. These are obtained by writing —A,, —c; for A; and é;; 
and making the substitutions w;; @ €;), wo <@ € and m; — Zp 1, elsewhere. 


THE SECOND ORDER MAGNETIC FIELD 
The fields of order k? satisfy the equations 
(56) VX H® = cD, 
(57) V-H® = 0, 
(58) VX h® = cd, 
(59) V-b® =0, 


with an incident field component given by 


(60) HY = im ( n,) . 
J 
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As before, solutions in the form of H® = F®-+G are sought. The particular 
function G® may be taken from Stevenson’s paper: 


(61) G? = — BALD 2° (Er) — T(E) — eo Baler? Jo (Es) +%2*Si (Bs) 
+3x3°I 33 (1) — Is (1) ] + €ocBoxrx2Js(€1) + €0c€ roxsJ1 (£1) + 0p /dx1, (P), 


where 


¥ = =. A ii, 


1 


Yi = (1/10)xsTu(s) 2 xs'— (1/5)xul i), 


*- du 


I = oe 
(€1) 2, R(u) 
~ du 
I = | — ee ’ P ’ 
11(&1) Je, (a;+u) R(u) ( ) 
C, = (a2—a3)[(a.+v)B+(a,+v’)B’], (P). 
In equation (61) the quantities B;, B are the coefficients of the first order 
electric field. 
The complementary function takes the form F® = V@®> where #® is 
the third degree exterior harmonic function. 
(62) oe = z x [OK (&:)E,+0,'K ( (&:) Ey'| + Exivox3K (€1:) + + D xi i(é), 
with 
= 2) x7/(@;+¥)-1, 
J 


_ \— oo A 
~ Je, (aitu)(u—v;)R(u) ’ 


. ee did 
Bib) + Lae ’ 


and the »,, v;’ are defined as the two roots of the equation 


(63) >> (1+28,;)/(aj;+v;) = 0. 


J 


E,, E,’, E, and D; are 10 coefficients which have to be determined. 
For the internal field we assume an expansion of the form 


(64) he? = 430 dinx pte: 
jk 


with di jx = d ix}. 
Substitution in the curl equation (58) yields 


(65) = €:jxdkje = cz E18 jer 
ace j 














1066 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


where €;,, is the permutation symbol; €;, = 0 if any of i,j, & are equal; ei, = 1 
if i, j, k are an even permutation of 1, 2, 3, and otherwise ¢;, = —1. 
Since d;, is symmetric in j, k and €;,, is skew in these indices, it follows that 


>: €1jx1n jt = 0. 
t.j.k 
Hence 


(66) » €:8 31 = O, 


which is satisfied in virtue of the electric analogue of equation (48). 
Taking the divergence of b® and equating coefficients of x, to zero gives 


(67) d Mish jx = 0. 


The normal field components are now calculated. When the coefficients 
of x13, x1%2?, x1%37, x1 etc. are successively equated to zero, we obtain 


(68) molgiE:/(aytvi) +g Ed /(aytvs’) +X ij] — (1/4a%) = Mindy j; 
- 
— (1/2a;)(1—6;;) =. Med ij = d;,/a;, 
k 


(69) molg:Eitg,E,’— (31,—A)D;/a,—Ui]+ (1/2a;) =. Mires = Nis 
- 

where 

as = ay [4myny2— 1 (A 1, + 0¢B\J2)+ (1 '10)A,(31,-A)], 


Xy2 = ay [Fmyn2?— t(Ayl2+ eocB, Ji) + (1/10) 4,(3—A)] 
+ (1/2a2)menin2, | 





= sh (P) 
X13 = ay" [4myn3?— 4(A 1134+-8¢e0cBiT33) + 3e0cBiJot+ (1, 10)4,(434,-—A)] | 
+(1 mitted 
U;, = (1/4a;) (A+ eocBis) + (1/10) A1(1;4+-2A —T/ay), } 
a (1/2a,—1/;)K;—A/ay?’, 
with K, = K,(0), and where A; is introduced because the variables x; are 


constrained by equation (1). 

The E,, E,’ are eliminated from equations (58) by multiplying by 1+26,,, 
summing on j, and using equation (63). This gives \; in terms of the co- 
efficients d;,x, thus 


(70) (vit lA; = a e( 2X ict 7 K-E Hidejj—3 >, Mindy ii 
J kj k 


—4ai; >) (1—8:)) a jder)/ay, 
F.8 


where y; = 2+4:/a2+a;1/a3,_ (P). 
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An additional relation is obtained by setting the coefficient of x:x2x3 equal 
to zero: 


(71) E(KB(0)—A*)+mynon3/a,+-menyn/a2+msynn2/a3 
tee) C.J./a; = wo De (m1 id ses Med 31 tMs 1d 112). 

Matching the tangential fields leads to a more complicated set of equations. 
From the terms in x,3/(ai+&), x1%2?/(a2+&), x1%3?/(as+£&) we have 
(72) (3/2)[E1.Ki/(ait+v1) + E1'Ki'/(ait+e1’)|+ Yu-— fdiut™1 = 2Mi/ai, 
(73) E,K,/(a2+11)+£1'Ky'/ (a2+01') + Yi2—3derz = Mi/ae, 
(74) E,K,/(a3+01) + Ey’ Ky'/ (as+01') + Yis—4dsis = Mi/as. 
From the terms involving x,°/(ai+€), xixe/(ai+&), %1%3"/(ai+€), and 
x1/(ai1+£&) we find 
(75) (3/2)[E.K1/ (ai +n) +£y'Ky /(ai+n1’)]+ Yu-tdiu—M1 = 2A1/a, 
(76) 3$[E1.K1/(a2+01) + Ei’ Ky'/(ao+01')|4+Z12—Fdizx = Ai/de, 
(77) 43{[21Ki/(a3+1) + Ei’ Ky’ /(a3t+1')]4+2Z13— 4diss = Ai/as, 
(78) —}(£,Ki+2)'Ky')+Vit3Dii— je. = — Ar, 


where 


VYi2 = 4monyn2+(1/10)AWh, 
Y;; = Lm3nyn3t+ eocB,J24+ (1 10) Aili, 
Zi2 = Lm yno?—1(AyIo+ €0cB,J;)+ (1 20) Ait, / 


Z13 = Lmyn3?— 3(A113+3e0cBi33) + (1/20) Ahi, | 
V; 1(4,J+ecB,I3) — (1 /10O) Al. / 


Yu = Lmyn,?— 3(Aq11+ €0cB,J2) + (3/20) Ah, | 


There are two more sets of cyclic permutations of equations (72) to (78). The 
coefficients A; and V/, have been introduced because of the restraints contained 
in equations (1) and (55) respectively. Vj; is needed to ensure that the two 
sets of tangential equations are consistent. 

An additional relation is obtained from the terms involving x;X2x3/(ai+€): 


(79) EK+mynon3+ ecC1J1—dies = 0, (P). 


By using equation (51) (electric) and equation (65) it may be shown that the 
cyclic variations of equation (79) are redundant. 

To determine A,, 1,, N; we proceed as follows: equation (72) is multiplied 
by 2 and added to equations (73) and (74). Equation (63) is used, thus 
eliminating FE; and E,’ and yielding 


(80) —48 > djy+2M1 = (11 +2)Mi/a1. 
d 


Similarly from equations (75) to (77) 
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(81) 1m, = 4AA it 1A J- deocBiA /a3— diy iad MN, = yiAi/ai. 
j 

Furthermore, from equations (72) and (75) 

(82) Ay = M,—a,.\,. 


Equations (80) to (82) may now be solved for Ai, M,, Vi. It turns out that 
only A; is required for subsequent use, since equations (73) and (74) can be 
shown to be implied by equations (76) and (77). It is found that 


(83) 4(yi tl) Ai = —H > di, j+4(mt Ay, 
J 


where 
(84) 4(yit1) Ai’ = a(r—2)"42q €3)2j2— Z, ade] 
\ j j 
ss 
+yi[mi—A (12+ Is) — 2c Bid /as) ( . 


If we now eliminate E,, E,’ from equations (68), (76), and (77) we get 


(Quogi/K1) {[(A1/a2) — Z12+ 3122] (a2+01') — [( Ai /as) — Zis+ 40133] (as+1’)} 
(85) a nae paid i12— 5 (A2/a1) a pusdin | (101/02) 


hayden } a; Mid ias— 4 (a3/a1) 2, sidin [a +r, a3). 





A similar equation also holds except that »; and »;’ are interchanged throughout. 
Two more equations of cyclic form are obtained when the other E;, FE,’ are 
eliminated. 

If E is eliminated between equations (71) and (79) we have 






; =¥ 
(86) po z; (urd 23 t+ Meid i31t+Msid 2) — mM ngn3/a,— Monn; /d2—M3NN2/A3 
t 





—e¢ 2) CJi/a; = (x a-1—A"/K) (di23—mnon3—eocC J). 





The above seven equations involve only the unknowns d;,. Together with 
the eight equations (65), and the three equations (67), we have a set of 
18 equations in the 18 dj. However, equations (65) and (67) are simple in 
form and may be used to eliminate 11 of the d;,,, and we are left with a 7X7 
set of equations in which all the matrix positions are non-zero. 

Once the di; have been determined the external coefficients E;, E;’, E 
follow readily from the tangential equations (76), (77), and (79). Knowing 
the E’s, the D;, e; may now be calculated from equations (69) and (78). If the 
e, s are eliminated, the matrix which must be inverted is of order 3X3. 

To find the electrical components of order 2, we make the substitutions 
A,-»A,, B,+-—B,, B-+—B wherever they occur explicitly, and make the 
substitutions y;; © €;;, uo  €0, mM; -> Zoli, Zoli ~> —m;, elsewhere. 
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THE FAR FIELD 


Since the exterior fields for an anisotropic scatterer have the same form as 
in the isotropic case, the expression derived by Stevenson for the far field may 
be applied here, provided that the exterior coefficients are modified as above. 

If R, 6, @ are the polar coordinates of the field point in the wave zone, then 
Stevenson shows that 


( p\\ g-kR 

Ey = ZH, = (26, 2. “y —— 
(87) 

= _(.1_ a a2) - 

i = (, indae “°a/) R’ 
where 
(88) P = 3k°S\?+k‘[—4S1"— (1/18) $3" + (1/18) S3""] 
and 


ge =-- (4 3) = a;A iy 
i” ae pia; +a,ax3 x gi/ai, 
(89) (2, 3) 2 z Pde 


S28) = (2/5) = A al ps aja; — (1/5) > a,(1 +25,) | ’ 
5 j 


t 


where 
a, = x,/R, 
and 
pi = —(6/5)[(a2+v) (astv)B+ (arty’)(as+v’)B'], | 
qi: = (2/5)[eoc(a2—a3)A,—3B)], | (P). 


n= (2/75) (3 (2a, —d2—a3)A1—5e9cB,—50D,].* J 


The expression for P is obtained from equation (89) by substituting electric 
for magnetic coefficients. 
DIAGONAL DIELECTRIC AND PERMEABILITY TENSORS 


When the dielectric and permeability tensors have the form 


€) 0 Bi 0 


*The formulae differ from Stevenson's because we use a slightly different definition for the 
coefficients and a different time dependence. 
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explicit determination of the fields becomes quite easy. It is found that the 
coefficients for the first term are 


~ 2Quom,A 
we) ot ime) Li 2ed * 
(91) A; 3 (uo— Mi) m: 


(ui— Mo) Li +2y0A ; 


For the second term we find 


a3(1 — 2/0) (nym2—ecA 113) +ae(1 —H3/mo) (nxm3+ecA1I2) 





oi P), 
(92) ms Ji[a3(1—p2/po) +a2(1—pus/u0)]—2a1A ( ) 
‘ 2 SiT/— '—SoTy' 
(93) 2BR (v) = - ve lhe 7; ’ 
where 
uj(aj+v) 
_ T= /uo)+A/v, ui T, athe” - 
T, = + 1=1,2, 
aitv Mo Qi > u,/a; 
j 
_ anna, Mijn sm; 
S; = (l—yi/uo)nmit- At 1 = 1,2, 


a ipo a u;/a; 
j 


the primed quantities being obtained by replacing v by v’ everywhere. 
The internal coefficients g;; may be readily obtained from equations (53) 
and (54), using equation (93) and the value of \ obtained from 


2d My 1 2BR’ (v)T, Bj 
a BP ip sso i Ae i. 
om Ho 2X a; Ho 2 wns st Ho ue aj;t+v 


2BR'Y')TY sy 
* Ho 7 aytv"" 


We shall not write down the complete set of coefficients for the third term 
since only the coefficients D; are required to find the far field. The D; are given 
by 


(95) 10[(u1—mo) Ji: +2u0A D1 = mi(uo—m) DY am i 


2 2903 — €3232 
+ Seda ete = ‘ ase ha, 44 2u—m) I+) (3a1:+a2+43)[(uo— Mi) L1— QUA | 
2uiai(11— 2. A)\ 


+ 40a 1A — > 
Vite 


2 a) i ). 
( +2e0¢ 1) (uo M1) +3 fo) (, (P) 

The electric coefficients may be obtained from the magnetic terms by 
carrying out the substitutions referred to above; the far field may then be 
evaluated by using equations (87) to (89). The resulting expressions are 
rather lengthy, so will not be written down here. 





ane 


——eeQQQreaa 
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THE MAGNETIC FIELDS OF A FERRITE ELLIPSOID! 


R. A. Hurpb 


ABSTRACT 


Approximate expressions are found for the internal and the adjacent external 
magnetic fields of a small ferrite ellipsoid under plane wave excitation. Con- 
sideration is given to the variation of apparent susceptibility with the size of the 
ferrite. 


INTRODUCTION 
In the preceding paper a power series solution in 1/A (A = wavelength) for 
the scattering of a plane wave from a small anisotropic ellipsoid was outlined. 
In the present paper these results are used to compute the magnetic field in 
and near a ferrite ellipsoid. 
A plane wave is assumed to be incident on a ferrite particle whose permea- 
bility has the form 
be —1k 0 
(1) w= itk be 0 
0 0 Ms 
with a scalar dielectric constant. Terms as far as \~' are found for a general 
ellipsoid; the term in A~? is computed for a spheroid whose axis of symmetry 
is in the x; direction, i.e., the direction of the static magnetic field applied to 
the ferrite. 


THE FORM OF THE FIELDS 
A small ferrite ellipsoid with permeability properties described by equation 
(1) is located at the origin of the rectangular coordinate system x; (2=1, 2, 3). 
The equation of the ellipsoid is 


3 
> x//a, = 1. 


i=1 


Upon the ellipsoid there is incident a plane electromagnetic wave 6, #” 
whose rectangular components are 


KH, = m,exp( —ik y a#,+iet) ; 
j 


C, _ Za exp —it 2. en 
j / 


where k = 27/X and Zo = (uo /€o)?. 


1Manuscript received March 3, 1958. 
Contribution from the Radio and Electrical Engineering Division, National Research 
Council, Ottawa, Canada. 
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The internal fields of the body are denoted by e, h; and the external scattered 
fields near the ellipsoid by E, H. 

According to the previous paper the rectangular components of these fields 
are given by 


(2) hy=c,tik > gu —B] } - dine mrte, +0), 
J jk 


(3) H, = a0 /dx,+ik[ao™ /ax,+Go]—k*[ab /ax,+G0?]+0(k*), 


with d;,,=d,;, and where 


gp” = z A x7 ,(é1), 


1 


| 
| o&” = BX,T,(&:)+B'X,'T,' (fs) +xwws D Bixi'Ti(és), 
t 


Oo 
] 


> «OE Ki(ti)+ D> x0/E/K i (E:) +ExxersK (é1) 
+ = Dx ,1,(&), 


Gr = eoc[ A gx2I2(E1) —AaxsI3(€1)], 


t 


nt : i ‘ 
| G\? = “whe £2 1(8)— Te) | Sea Baber’ FG) +0280 
| 


+-3x3°Is3(21) = 13(€1) | +c Boxix2J3(€1) +eocxex3(a2—a 3) 
| X [ (arty) B+ (ai+v’) Bi (&:) +8/dn. As, 


and in which 
vi = (1/10)x,7, (1) a x —(1/5)x.0(é:). 
J 


In the above expressions, c is the velocity of light in vacuum, & is the 
radial-type ellipsoidal coordinate (¢;=0 is the surface of the ferrite), and 
the coefficients ¢;, gi), dix, €:, Ai, B, B’, Bi, Ei, E’, Di, A,, B,, B, B’ are 
independent of the variables x;, £:. In equation (4) the x; components of the 
vectors G) and G® are given; the other components are obtained by cyclic 
permutation of all subscripts excepting of course the subscripts on €o, #o, and 
£,. The remaining functions are given by 


* du : du 
I(é) = J R(u)’ Ti) = J (a,tu)R(u) ; 






Pen _du el fj {echude 
hel jee. ae FP 
- a du ‘ _ { du 
ia (u—v)R(u) ’ K(f) = | R*(u)’ 

: P du 
Rh = J (a;+u)(u—v,)R(u) ’ 
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2 2 
=R 0] z, Xi Hoctr)-1 | ’ 
i 


6, = c® xj /(as+v,)—-1, 
Jj 


X 







where R(u) = [(a1+u)(a2+u)(a3+u)]*. The range of the integrations is from 
&, to ©. The primed functions 7,’, K;,’, etc. are obtained by substituting v’ and 
v, for v and v_ respectively. The quantities v, v’ are the two roots of the 
quadratic equation 







Zz, (a,;+v)™ =0 





and »y,, v;’ are the roots of 






Dd (14+26,;)(a;+v,)* = 0. 


j 






The electric field components for the first two terms (which are needed 
to compute the higher order magnetic fields) are given by expressions identical 
in form with equations (2), (3), and (4) except that the coefficients ¢;, gi;... 
are replaced by é;, 9:;... etc.; and the vector G™ is replaced by G®, G™ 
being obtained from G“” by making the substitutions A ,— — A ,and €) — uo. 










THE ZERO ORDER TERM 


A, are found to be 


p Suef atl ins) 
cq, = Ao i, + I. ’ 





The coefficients c;, 






























ie a ee mf) 
©) eee a | 
3 = ___2poms = 
5 Quo (us—mo)Ls/A ’ 
(6) A, = (¢;—m),)/Ii;, 
where 
fi = p—wo(1-—24/7)), 
Ay = fifo—x?, 


A = (ajaza;)-3. 


The electrical coefficients ¢;, A; are 





a 2d : 
i cl J; ’ 
> _ i(1—e/es) 

A = — f 

: a7,” 


where 
f; = e—e(1—24/1;), 





SEES OPT RE eT PEEP 27TH 
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and 


I, = I,(0). 


(This notation will be adopted for the other elliptic functions defined above.) 


THE FIRST ORDER TERM 
It is found that 
ByJ, = (SiR2+ixasS2) (asx? — Ri Re)", 
BoJ2 = (S2Ri—ixas3S;) (a3’x?— Ri Re)“, 
B;J3A4; =O3(WiW2! — W2Wy') —ixT, Y(O2W1' —O,W2') —ixT,' Y’ 
(0\W.—02W,), 


2BR?(v) Ai = R3(O2Wy' —O,W2') —i(«/p)O3(a1 We’ +a2Wy’) + (x?/u)T,’ Y’ 
(a,;02+a20;), 


2B’ R*(v') Ax = R3(O,;W2—O2W,) +4(«/n)O3(a1W2+a2W,) — (x?/u)T,Y 
(a,;02+a:0)), 
where 


Si = Q2(u3— Mo) (nom3+ €ocA 112) +a3(u — Ho) (n3m2— eA Is) 
+ ixas3(n3m1+ eocAols) , 


Ss = Q;(u3— Mo) (nym3— eocAoI;) +a3(u—o) (n3mi+ eocA eI) 
- ixa3(n3m2— €ocA Is), 


Ri = 42(uo— ms) +43(uo—H#) —2u0A/Si, 

Re = a3(uo—H) +41 (uo — Ms) — 2A /Je, 

Rz = (uo—s) (ai+a2) —2uo0A /Js, 

O1 = (1—po/us)asnsms—ai[ (1 —Ho/m) nim —i(«/p) (myme+ oc Asli), 
O2 = (1—o/us)as3m3—Ao[(1—po/n)mom2+7(x/p) (nym, —e0cAsl2)], 
Os = (u—po)[aon2m,+aynym2+ e0cA 3(a1]1— Gel 2)] — ix (a2n2m2—aynym)), 


a,[(uo—u) T+ u0A /v) ei as{ (uo—mus)7>+u0A /v) ba: 3 Si 


me u(a,+r) us(as+) 
seg Ae a 
-s Qo+tv ai+y . 
Ay = —R3(WiWe! — W2W1')+ (x?/u) [Ty Y’ (a1W2+02W1) 
—T,Y(a,W2'+a2W;’)). 


The internal coefficients g,, are calculated from 
£1 * —nym,+ cA 3I2+ BJs, 
£21 Ziat ces, 


ll 
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and cyclic permutations thereof. The coefficients g;; are given by 
£ii = —nm;+ (a3/a;) (1 —po/us)n3m3+g(v)+g(v’), 


where 


ee ete 


The electric terms are 


zs anil aes 7 4 thoes 4g | 
kc vT y+ te (T,+A4/v) }. 


€—€)[Zo(aanalst+asmsl2) +mocA 1(a3/3—a2T2)] 
(e—€0) (@2+a3) +2€04 / Ji 


with B:, B; cyclic permutations of By. 


BJ, ian 


B = Z,(e.—e) 2. ndlertry* | 6r—r')[¢—e0) 00d /v] 


B’ = same as B with v and p 


, 


interchanged. 
Ji2 = —Z nol, —pocA sle+ B:J; 


and cyclic permutations thereof. The remaining g;; may be obtained from: 


Ji2—GJ2a1 = b3CC3, 

G23—G32 = C(uci—ixce), 

G31 —Jis = C(tKc; +c), 
. 2BR*(v)T, , 2B’R*(v’)T,’ 
fiec= ~Karen 


aitv a;+v’ 


THE SECOND ORDER TERM 


Computation of the second order term for a general ellipsoid would involve 
considerable work; however, a reasonable simplification occurs if the calcu- 
lation is done for a spheroid with a;=a2. In this case the quantities v,, v;’ 
defined above become 


y= ve HHS vw = —A), 
py’ = Vo! = — (1/5) (a;+4as), 
v3! = = (1/5) (3a,+2a3). 


Furthermore it may be assumed without loss of generality that the direction 
of incidence lies in the x1—x; plane, and thus m.=0. 
The following abbreviations are introduced: 


B; = Mola co —2py;"! —2A (aw?K,)—"], 


B,/ = same as §; with v;’, K,’ replacing », K,, 
Ni = €¢[usGox—ugs2+ixgail, 
No = €€[ugsi—wsGistixGs2], 
N3 = €¢{u(Gi2—Ge1) +2xGss], 
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Ar = (81—3/a1)?— (3x/a1)?, 
As = }ar'{3us[4+ (a1/as) (us/m)]—Bs'[2as+3a.ms/x]}, 
As = [83— (2u/a1) —u3/a3]?— (2«/n)?, 

As = 3uslstuols, 


Qi a Mit 34(k/p)No— deusc(1 — 41/03) G23 
+ }uo[m,—2nym3n3(1 —4a,/as)| 
—pol{AA i+ 4A Wi+ fe0cBi[Ji(1 —a;/a3)+A/as}} ’ 












Q2 = = 1No— 41(k/m)Ni+ASeusc(1 —1/03)91st+ fuome 
—po[AA ot 14 oli+ LeocB.A /ay), 






X; = $ €u3CG23/A3st tu0(min3?/a1+2m3nin3/a3) +uol — jar! (Ail 
+3¢0cByI33) +3 0c BiJ1/a3— (1/20) (11 +J]3)A1/ay}, 


X, = — 5 €u306913/ast+ 4 uomens"/ay — pol ja, ! (A 213+ €9¢B2J1) 
+ (1/20) (71 +J3)A2/ay}, 


Zi = 1myns?— 3 (A113+3¢€0cBil33) + (1/20) Arh, 
2. = Lmons? — }(Aol3+ e0¢BoJ:) + (1/20) A oli, 
Py = }ceai{ (e—€0)(3+41/a3)Gs2+ €(1 —us/u)G23—ie(K/u) Gar}, 









P. = 4c{ (€—€) (3a3+a1) G31 — €€1 (1 —3/u) G13 —t€01(K/u) G2}, 


Ps = 4c[4(1+a3/a1) (€— €0)Go1+2€(a3/a1) (G12 -—Gar) +3.N3/u] 
+ (3+2a3/a1)[— 4mgny?+3 (1A 3+ 2e0cBsJ3) — (1/10)13As], 


Ly = (a2—a;)[(ai+v)B+ (ai+v’)B’], 

Sy = (4+43/a1)(8:'Z;-—X )+a371(Qi—Bi'P;), «1 = 1,2; 
ry = — }ix[4ar-'+us(asu)~"], 

ro = —pay!— (11/4)psas-'+ 381 [4+ams(asu)"], 










Ty = ixecay *(G12+Go1) +uecay Gss—msecas Goo— eouoc x Lj,J;/a; 





—pomngn;/a1+ BseocL iJ, 





T2 = 2eca;—"[u (G12t+Ge1) = ixGss]+2eouocBs(Jiar~ —J3a3;“") 
— B3€9¢L1J1 — wo(mgny7a3-'+2myn3n,a1"'), 









Ag [ (83 — Quay! — pss!) T2 — 4ixT1/ay], 
— (1/5) [Pit (a3—a1)Z1]— 3 (Ai + ec Bis) + (1/10) Ail, 
— (1/5)[P2+ (a3—a1)Z2] — 3 (A oT + eo Bol) + (1/10) Aol, 
— (1/5) (Qi: + (a@s—a1)X1]— (1/20) wo 1 (121, +513) 






s 
ll 










_ duoeocas !ByI3, 
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= (1/5) [(Q2+ (a3—a1)X2] aa (1/20) 40A 2(12I,+5a3I3/a1) 


—dyoeccar' Beli, 


— (1/20) ¢[2(a1+as) (€— €0)G21 t+ €(@2—a1) (G12 —Gar) — 1 (x/m) €aiGss] 
+ (1/5) (as—ai)[}mgmy?— $A gl — €0¢B3J3+ (1/10)A als] 
— (3/20)A 31 —heccBslh, 


— (1/20) woms[ms°+ (a1/as) (1 —ms*)]+ (1/20) (a1/as) (us/u) Ns 


— (1/30) eomocB;(Ji+4a,J3 /a3) oe (1/50) woA 3[7 (1; —a;I3/a3) —3I3 
—22a,I;/a;]. 


The coefficients then are 


di33 = 


d122 


dit 
d233 
doit 


d222 
A3d333 


Audie3 = 


d311 
d322 


e1 


(riS2+72S1) (71? +727), 


—43(us/u)disst+4N1/u—}$eCGs2— tA2! (81 —3y0/a1)[ (81-3 /a1) 
Xm n,? —3ixmsn;*/ay], 


—dy22— (us/u)disst+ Ni/n, 
(ro52—7151) (11?-+ 727)“, 
—}(us/m)dosst+ tN 2/ut :ecG§sit As! (B1—3y0/a1) 
X [(61—3u/a1) mony? + 31nmn,?/a4], 
—do11— (us/u)des3+ No/u, 


3a37!{ tu3N3/ut tuoms(4n3a3/a1— 17) +uoA 3[ (1/10) 11 (7 —2a3/a1) 
= 47303/a1]+uoeccB3(J3—Jia3/a1)} —B;'P3, 


(83—2yay-! — was!) T14+-1xT2/a1, 
3(T— (us/u)ds3st+ N3/u], 
$(— T—(us/u)dsss+ N3/u], 


— (1/10) (@3—@1u3/m)dis3+2A07'{ filuoMi (1 —24/I1) —a1M3] 
+ix[uoM2(1 —2A/TI;) —aiM,]} ’ 


— (1/10) (a3—aips/u)d233+2A0! {filuoMe2(1—2A/lh) —a1.M,!) 
—tx[uoM,(1 —2A/I;) —a,M3}}, 


— (1/10) (a3 —ayu3/pm)d333— As! (2uol1Ms+a3I3Me), 
M+ }3e1+ (1/20) (a3—a@ims/u)diss, 
M2+ $e2+ (1/20) (a3 —@1u3/u)d233, 
As! (usMs—a3M¢), 
(1/5) (@s—a@1) (4+01/a3)—!{ [1+ teims(@su)-"]disst+as-'P1 
— (4+4;/a3)Zi}, 
(1/5) (a3— a1) (44+41/a3)~1{ [1+ fa.ms(asu)—"]des3+a37'P2 
— (4+4:/a3)Z3}, 





en ee eee 


a 
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3K;'E;’ (1/5) (@3—a@1) (8-+2a3/a1){ [3 (a3/a1) + $us/uldsss+P3}, 
Ey = E, = E3 = 0, 
EK = dy23— €9¢L,Ji. 

The electric components may be obtained from the above by making 
the substitutions given in the previous paper. The detailed results will not, 
however, be given here. 

APPARENT SUSCEPTIBILITY AS A FUNCTION OF PARTICLE SIZE 
Of current interest is the apparent susceptibility tensor which is defined 


as the ratio of the internal magnetization to the applied magnetic field. Thus, 
denoting the apparent susceptibility components by f—wo and &, one has 


b:— woh = (f— wo) m1 —ikm2, 
bs — woh, =1Km1+ (@—po)mM2, 


where 4; is the internal magnetic induction. Using the static approximations 
for b; and h; it is found that 


Quod : 2 
i— po = ae [((u—mo)fi—« J, 
(7) 1440 
_  Qpowd 
. Ta, Uituo—al. 


These relations are true when a;=4@2; if ai a» it is found that the apparent 
susceptibility no longer has the antisymmetric form of the ordinary sus- 
ceptibility matrix. 

From equation (7) it follows that, for a spheroid, 


Quod +«— 
eg ee SEO, aR BO 
(8) fitk—po l, ae 


This reduces to the values given by Spencer, Ault, and LeCraw (1956) in the 
particular cases of sphere, cylinder, and disk. 

To find how the apparent susceptibility behaves as a function of particle 
size, higher order approximations to the interior fields are necessary. It is 
clear, however, that since the higher order magnetic components are functions 
of position within the ferrite, some other approach than the one above must 
be used. Two methods of handling this problem have been devised. Artman 
(1956, 1957) defines the apparent permeability in terms of the exterior fields, 
using the dipole moment of the particle as the fundamental quantity. The 
other approach, used by Tompkins and Spencer (1957), makes use of the cavity 
perturbation formula to define the apparent permeability. In this paper both 
methods are employed. 

Consider two incident waves, 17 and 2% proceeding in the +x; direction, 
whose only non-vanishing field components are described by 


ise: mol, f=1 


2H: ms, = 1,1, = —1. 
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Under this form of excitation one finds that, from equation (6), 
AyI,= — (utK—po) (fit), 
and thus, writing Yo=ftKk—yo, and using equation (8), it follows that 
Yo= — 2A A. 


Taking the third term of the series into account, the apparent susceptibility 
at finite size yY may be defined by 


(9) y= — 2A (A 1—k*D,), 
so that for a spheroid, using the value of D,, it is found that 
(10) v = vot (1/5) wok? (1 +a;/a3){ (Wo/mo)?[401(€/€0) A—*(11 +I)? 
+ (€—€0) (@313—@11;)?(2a3J1A Ri) + 3A-! (6a,11(1+<a;/as) 
+ 2a313—a111?/A)|+ (Wo/uo)[(€/€0)a1 (Li +Ts)/A — (€— €0) (@303— ail) 
X (1—ai/as) (RiJ1)~!— 3 (601 — 2a,1;/A + (4/5)as+ (6/5)a12/as)] 
+ai(e— €0) [e071 +2A (RiJi)—]} ’ 


where 
R, = (€—€0) (@1+@3)Ji +2604. 


For a cylinder of radius a, length c, a « ¢ formula (10) becomes 


(11) Wb = oll — (2/25) k2c?]+ (1/20) wok2a?[ (1+ €/€0) (2+Yo/uo)? 
+ (8/5) (Yo/uo) siti 8e(e+ €9)7! (24+3y0/puo)]+O0(R?a'/c). 


For the disk of radius a, thickness c, a > c we obtain 


(12) = Poll — (3/25)k2a?]+ (1/20) mpuok?ac[3 (Wo/mo)?+ (1 — €0/€) (Yo/mo)] 
+O0(k?c?). 


The case of greatest current interest is that of the sphere, for which we find 


(13) = vot wok%a*[ (1/90) (5+€/e0) (3+ 20/10)? —Wo/mo—3/5 
—43(€—€) (2e+3e0)—'], 


a being the radius of the sphere. 
If the static magnetic field at which yo is resonant is denoted by Hp and the 
resonant field of y by H; then from equation (13) it is found that 


(14) H =~ Ho+4nrM (5+ €/e€0) (2/45) Ra’, 


where 47M is the magnetization. Formula (14) is derived under the assump- 
tions that the imaginary part of Wo is much larger than its real part near 
resonance, and that the resonance shift H1H— H) is much less than half the line 


width. 
The method of Tompkins and Spencer makes use of the cavity perturbation 


formula 


(15) 3 t = (utk—ypo) fn . A *dr 
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where f is the frequency, E is a constant, and % hare the circularly polarized 
incident and internal magnetic fields. The integration is taken over the volume 
of the ferrite. Equation (15) holds when the sample is placed at a zero of 
electric field. For a sphere of radius a, it is found that 






(6) BE a 88 eae) JF ert (durtdiatdaster) |} 
“7. 3 U+K—Mo a é1 10 111 12213361) |f- 


As a — 0, one has, 








Defining the apparent susceptibility at radius a by 





3E_ of 


(17) y Geet f 






and substituting for the coefficients in equation (16), one obtains 
(18) WY = Pot ka*pol (1/45) (5+ €/€0) (83 +2~0/u0) —3/5). 


This formula, though different in form from equation (13), leads to approxi- 
mately the same resonance shift formula as equation (14). 








DISCUSSION OF RESULTS 
It is of some interest to compare our results on the susceptibility size effect 


with those of other investigators. Artman (1956, 1957),* using an approxi- 
mation to the actual ferrite matrix, obtains a formula which in our notation is 


y = Yot (1 /90) (€/€o) wok?a? (3+ 2Yo/ uo)? +O0(R4), 











(19) 









from which he derives the resonance shift formula 


(20) H =~ Hyt+4aM (c/€5) (2/45) Ra’. 















Thus our formula (14) predicts a resonance shift greater than Artman’s by 
(5e9+e)/e, which for most ferrites is ~ 1.4. 
Artman has also proposed the empirical formula 


(21) H =~ Ho+4nM (e/e0) (4/45) ka?, 







which fits well with the available experimental data. It will be seen that our 
result, while not in agreement with experiment, is somewhat better than 
Artman’s theoretical formula. 

Ament and Passerini (1956) have also used the exterior fields to compute 
the size effect. They find a resonance shift formula with a coefficient 9/8 that 
of equation (21). However, it would appear that their result applies to 
intrinsic rather than apparent susceptibility. Indeed, when their analysis is 
applied to the apparent susceptibility, a resonance shift with a coefficient 
















*There are some errors in equations (23) and (24) of Artman’s 1957 paper. The mistakes 
do not appear, however, in the 1956 paper. 
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(eo+e¢)/e that of equation (20) is predicted. This point seems to have been 
overlooked by Artman in quoting the Ament—Passerini results. 
The formula derived by Tompkins and Spencer (1957) has leading terms 
which are 
(22) Y = Pot (1/45) k?a?(€/€0)Yo(3+2p0/u0). 
This should be compared with equation (18). From equation (22) can be 
derived a resonance shift formula which is almost identical with equation (20). 
In conclusion, it is expected that equations (13) and (18), being based on a 
solution rigorous to O(k?) for plane wave excitation, will be more accurate 
than the other formulas quoted above; and this is to some extent borne out 
by the better agreement of the resonance shift formula (14) with the experi- 
mental result of Artman (equation (21)). As far as choosing between the two 
approaches to the problem, as represented by equations (13) and (18), it 
would seem that the cavity perturbation formulation adheres more closely 
to the way in which the ‘‘actual’’ apparent susceptibility is measured, although 
there is little difference in the predicted resonance shifts to O(k?). 


APPENDIX 


Values of the Elliptic Integrals 
When a;=d», the various elliptic integrals reduce to elementary functions. 


Thus if a3 > ay 


4 4 

a 3 +(a3—a;) 

Pa tates ij 
(a3—a) 8 as (a3—a,)' 


When a; < a 
I = 2(a,:—a3)~4cos!(a3/a1)}. 


The remaining functions are then found from the relations 
2, 1¢ 2A, 
1 
2 a,l; = 7 


Jy = (I2—Is3)(a3—a2)~", ete., 
8lutJot+J3 = 2A /ay, etc., 
T.=4D Marty) '-A/y, 
“ A(t tg te) a 
K = (J2—Js)(a3—a2) 
For a sphere of radius a 
fe Qa. Wee sepa, 
J, =14= T, = T/ = (2/5)a>, 
K = K, = (2/7)a~. 
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EXPERIMENTAL EVIDENCE FOR AXIAL VECTOR INTERACTION 
IN THE DISINTEGRATION OF Pr‘! 


R. L. Granam, J. S. GEIGER, AND T. A. EAstwoop 


ABSTRACT 


A detailed study of the disintegration scheme of Pr has been carried out using 
lens spectrometers, scintillation spectrometers, and fast coincidence techniques 
The 2293 kev B-component (1.3%) studied in coincidence with y 691 has a pre- 
dominantly first forbidden unique shape and the 6-y directional correlation is 
strongly anisotropic. The results indicate a probable spin-sequence 0—(8) 
2+(7)0+, but the accuracy is not sufficient to exclude completely a 1— assign- 
ment for the Pr'* ground state. The 8 803 component (1.0%) observed in 
coincidence with y 2181 has an allowed shape and the 6-y directional correlation 
is isotropic. The shape of the total B-spectrum was studied in a double lens 
spectrometer having good anti-scattering properties using sources of Ce!*# oxide 
sublimed onto 200 pg/cm? Al leaf. Subtraction of a unique (B;;) shape 2293 kev 
1.3% component and an allowed shape 803 kev 1.0% component yields the ground 
state @-spectrum. Analyses were made using the electronic functions of 
Zerianova which assume uniform charge distribution with p = 1.2 A” 10-83cm 
as well as the point charge functions of Rose with p = 1.41 A? 10-8 cm. The 
observed spectrum shape can be explained on 0— to 0+ selection rules only by 
the axial vector interaction. Limits on the permissible tensor or pseudoscalar ad- 
mixtures with the axial vector are discussed. No acceptable fit is possible with 
tensor, pseudoscalar, or any TP admixture (interfering or not) using Rose and 
Osborn’s pseudoscalar formalism. 


I. INTRODUCTION 


Since the non-conservation of parity in B-decay was demonstrated, by Wu et 
al. (1957), the nature of the B-interaction has been the subject of renewed 
speculation. The He® recoil experiment by Rustad and Ruby (1955) resulted 
in the general acceptance of the tensor interaction as the principal contributor 
to GT decay. The more recent recoil experiment on A**, when considered in 
conjunction with the recoil experiments on Ne” (Herrmannsfeldt et al. 1957) 
and the neutron (Robson 1955 and 1958), and the measurement of the - 
asymmetry in aligned neutron decay (Burgy et al. 1958), require for their 
explanation a dominance of the axial vector interaction in GT decay. Re- 
examination of the He® measurements (Rustad 1958) revealed several pre- 
viously unsuspected sources of large systematic error which rencer the inter- 
pretation of this experiment ambiguous. The measurements of the Pr'* ground 
state spectrum presented in this paper provide additional evidence in favor of 
a dominating axial vector contribution in GT decay. 

The study of 6-spectrum shapes does not, in general, provide conclusive 
evidence about the nature of the @-interactions. However, for O— to 0+ 
B-spectra the interactions which can contribute, 7, A, and P, give markedly 
different spectrum shapes. A single nuclear matrix element is responsible for 
the possible pure tensor or pure pseudoscalar contribution to this type of 
transition and unique spectral shapes are predicted for these cases. 

‘Manuscript received May 8, 1958. 
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Pr'## (17 minutes) appears to be the most convenient nucleus for an ex- 
tensive study of a O— to 0+ #-spectrum shape. Its decay scheme is shown 
in Fig. 1. Pr has a long-lived (285-day), low energy (0.3 Mev) precursor 
Ce, The ground state 8-transition, 0— to 0+, is energetic (2.98 Mev) re- 
quiring a minimum of corrections for counter sensitivity and source thickness. 
The lower energy components of the total B-spectrum are weak and can be 
corrected for with adequate precision for shape analysis. 


144 144 
Pr Nd kev ASSIGNMENT 


2984 O- 


803 t10 
1.0% 


2293 15 
1.3% 
(0.26+.06%) 


2984110 (0.78t0.16%) 
97.7% 


0 O+ 


Fic. 1. Disintegration scheme for Pr'*4 showing radiation intensities derived in this work. 


Experimental evidence for a 0— assignment for Pr'* has been based on 
studies of the 2293 kev 6-component taken in coincidence with y 691. Ketelle 
and Brosi (see Laubitz 1956) and Laubitz (1956) both find agreement with a 
first forbidden unique shape. The 691 kev level in Nd'* has been assumed to be 
2+, as for almost all first excited states of even—even nuclei (Scharff-Goldhaber 
1953). The y 1490—y 691 angular correlation studies of Alburger and Kraushaar 
(1952) and Steffen (1954) are consistent with a 1—2—0 spin-sequence, but 
provide no unambiguous information (see Section II). Brazos and Roberts 
(cf. Steffen 1954) find evidence from polarization-correlation studies for a 
parity change in this sequence which leads to the 1— assignment of Fig. 1. 

Various values of the 6 branching ratios for the 803, 2293, and 2984 kev 
components have been deduced experimentally. The y-ray intensity measure- 
ments by Emmerich, Auth, and Kurbatov (1954) give 1.9%, 0.8%, and 97.4% 
respectively and those of Kreger and Cook (1954) give 1.26%, 1.44%, and 
97.3%. Estimates made from Kurie plot analysis (cf. Porter and Cook 1952; 
Alburger and Kraushaar 1952) can have only qualitative significance since 
the shape of the dominating ground state B-component is in question. 

Magnetic analyses of the Pr'* 8-spectrum have been made by four sets of 
experimenters. Porter and Cook (1952) and Emmerich, Auth, and Kurbatov 
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(1954) found their results consistent with a 96% ground state spectrum 
having an allowed shape. Alburger and Kraushaar (1952) and Laubitz (1956), 
however, found more low energy electrons, i.e., their Kurie plots were concave 
upward. Alburger and Kraushaar interpreted their spectrum in terms of an 
allowed shape 90% ground state component as well as two lower energy 
components each with ~5% intensity. Laubitz, on the other hand, concluded 
that the ground state spectrum had a non-allowed shape. 

The spectral shape predicted for a O— to 0+ transition assuming pure 
tensor (7) interaction disagrees with all of the observed shapes. Laubitz 
introduced an interfering pseudoscalar (P) contribution with a very strong 
coupling constant, Cp/Cr; = 107, in order to obtain theoretical agreement 
with his data. Zerianova (1956) found that a similar PT admixture was 
necessary to fit the essentially allowed shape data of Emmerich, Auth, and 
Kurbatov (1954). More recently Alaga et al. (1957a) have extended the 
pseudoscalar theory of Rose and Osborn to include nuclear force effects. 
They also found that a large pseudoscalar contribution was necessary in order 
to fit qualitatively the experimental data of Laubitz with a tensor pseudoscalar 
admixture. 

The experimental work described below was undertaken since the spectrum 
shapes derived from the various earlier experiments are not in agreement. The 
directional correlation experiments were carried out in an attempt to establish, 
in a less ambiguous manner, the assignment of 0— for the ground state of 
Pr’, The evidence described below strongly supports this assignment but 
the possibility of a 1— assignment cannot be completely excluded. It is 
possible, therefore, that the conclusions from the spectrum shape analysis 
may be invalid since it is quite easy (cf. Laubitz 1956) to account for the 
ground state spectrum shape on the basis of a 1 — 0, yes, type transition. 


Il. DIRECTIONAL CORRELATION MEASUREMENTS 


(a) Experimental Apparatus and Procedure 

The directional correlation measurements were carried out using the 
experimental arrangements shown schematically in the insets of Figs. 2 and 
3b. y-Rays were recorded in 2 in. X2 in. NaI(T1) crystals and the 6-rays in a 
3cm diam. X 1 cm thick trans-stilbene crystal covered by 2 layers of 2 mg/cm? 
Al foil. These scintillators were mounted on RCA-type 6342 photomultipliers. 
The photomultipliers were operated at high voltages (> 2000 v) and were 
connected to a fast-slow coincidence circuit (27> ~3X10-§ second) similar 
to that described by Bell, Graham, and Petch (1952). 

Gain stabilization was used on the counters in order to minimize systematic 
errors. The stabilization method used for the 691 kev y-counter was the 
following. A collimated beam of y-rays from a Hg™* source near the y 691 
counter produced a strong photopeak at 280 kev which was monitored con- 
tinuously by a dual differential pulse height analyzer connected to two 
counting-rate meters. The two pulse height windows were set one on each half 
of the intense 280 kev photopeak. The two ratemeter outputs, connected in 
opposition, provided a difference signal which when applied to a variable 
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Fic. 2. The directional correlation between y 691 and y 1490 taken with the experimental 


arrangement shown in the inset. The data taken in the 180° to 270° sector were averaged 
with those taken in the 180° to 90° sector. The combined data have been corrected for random 
coincidences and for the small variations in the individual counting rates with angular setting 
(< 1%). Standard deviations are shown by the error bars. The theoretical curves allow for the 
experimentally measured counter solid angles. The coefficients listed are the least squares 
fits to the experimental data after correction for counter solid angles. 
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(< 2%). The standard deviations are shown by the error bars. The theoretical curves in 3(a) 
(at left) allow for the measured counter solid angles and for the 8 energy window shown 
assuming pure )_|B;;|2(,-2) 8-radiation. The coefficients listed are the least squares fits cor- 
rected for counter solid angles. The theoretical values shown are for the 8 energy window used 
in this experiment assuming pure > ,|B;;|*(z~2) radiation. 
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gain first amplifier stage kept the over-all gain of the system constant. In this 
way the over-all gain was kept constant to ~1% over long periods of time. 
This system can accommodate photomultiplier gain drifts of as much as 50%. 
A similar arrangement using the 661 kev y-ray from a Cs'*” source was used to 
stabilize the y 1490 counter. 

A different type of gain stabilizer was used with the 8-counter. Instead of the 
differential method used for the y-counters a simpler integral method was 
employed. Here a single counting-rate meter is fed from a discriminator set 
high on the 6 pulse height spectrum. The counting-rate meter output voltage 
is compared with a preset voltage and the difference is applied as above to a 
variable gain amplifier stage to maintain a constant counting rate. This 
method is applicable only to long-lived sources where decay corrections are 
unimportant and when background fluctuations can be ignored. 

The directional correlation apparatus was run continuously in a 17-step 
automatic program. At the end of each 10-minute counting interval the angular 
setting was advanced through a 15° increment. The angular range was 
from 90° to 270° and provision was included for measuring random coincidence 
rates at the 105°, 180°, and 255° settings. The coincidence data were printed, 
along with a code indicating the angular and coincidence delay settings, 
‘by a solenoid-operated electric typewriter. The counter channel rates were 
monitored separately with counting-rate meters and pen-recorders. 

The sources were in each case carefully aligned and the net coincidence 
rates were corrected for positioning errors by normalizing to constant counting 
rates in both counters for each angular position. The corrected true coincidence 
counting rates for pairs of angles symmetric about 180° were averaged. 

The averaged data were analyzed by the least squares method to obtain the 
coefficients Ao, As, and A, of the Legendre polynomials appearing in the 
assumed expression for the correlation, 


N(0) = Ap t+A2P2(cos 0)+A4P,(cos 6). 


The resulting coefficients were then corrected for the finite solid angle sub- 
tended by the two counters in the manner suggested by Lawson and Frauen- 
felder (1953). 


(b) Preparation of the Source Material 

Cerium 144 and praseodymium 144 in radioactive equilibrium were used for 
all the experiments. A ‘‘cerium concentrate”’ extracted from fission products 
several years old was supplied by the Commercial Products Division of Atomic 
Energy of Canada Limited. The concentrate also contained some La, Pb, Zn, 
and other elements so that further purification was necessary. The trivalent 
cations were separated from the other elements by fluoride and hydroxide 
precipitations. They were then separated from each other by cation exchange 
using ammonium a-hydroxy isobutyric acid solution (Choppin and Silva 1956). 
Those fractions of the eluate which contained Ce! were adjusted to pH 1.0 
by the addition of conc. HCl and were passed through another cation exchange 
column. The absorbed cerium was eluted with 7 VN HCI (Diamond et al. 1954). 
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The eluate was evaporated to dryness and the cerium chloride was redissolved 
in 0.1 V HCl. Sources were prepared from this solution. 


(c) y 1490-y 691: Results and Discussion 

The experimental arrangement used to measure the y 1490—y 691 cascade 
correlation is indicated in the inset of Fig. 2. Both of the 2 in. X2 in. Nal(T1) 
counters subtended a solid angle of about 1.5% of 4 at the source. The 50 uc 
Ce'—Pr' source used was prepared by evaporating Ce™ chloride to dryness 
in a 1 mm I.D. glass container. The counter outputs were fed through pulse 
height selectors adjusted to record the respective photopeaks. 

The experimental data yield the counter-corrected coefficient ratios A./A» 
= —0.24+0.02 and 44/A,) = —0.02+0.04. The more probable spin-sequences 
which, with suitable multipolarities for the y-rays, give theoretical correlation 
coefficients agreeing with those observed are listed in Table I. Steffen (1954) 
has measured this directional correlation and has obtained the values 
A2/Ay = —0.210+0.015 and Ay/Ao < 0.01. He interprets his results in 


TABLE | 
y 1490-7 691 DIRECTIONAL CORRELATION 


y 1490 y 691 
Spin-sequence multipolarity* multipolarity A2/é 





0 OD —0.25 
Q —0.25 
D 


0 —0.25 
0 D+(+6)0 

0 D+(—0.22)?0 —0.24 
0 D+(+0.037)Q —0.24 


Experimental results (this work) —0.24+0.02 —0.02+0.04 





*Amplitude and phase according to Biedenharn and Rose, 1953, Revs. Modern Phys. 25, 729. 


terms of a 1—(y 1490)2+(y 691)0+ transition sequence in which y 1490 
is E1 radiation with a 0.1% admixture of M2 radiation. His results are con- 
sistent with our present measurements and the 1— — 2+ —0-+ spin-sequence 
is incorporated in the decay scheme of Fig. 1. 


(d) B 2293-y 691: Results and Discussion 

The 8-y directional correlations were measured with the counter arrange- 
ment shown in Fig. 3(6). The 2 in.X2 in. Nal(Tl) y-crystal and the 3 cm 
diam. X1 cm thick trans-stilbene 8-crystal subtended solid angles at the 
source of about 1.5% of 4x. A 10 uc open source of Ce!*—Pr'!, prepared by 
evaporating a solution of Ce™* chloride on 20 mg/cm? Al foil, was used. The 
source was mounted on the axis of rotation of the moving 6-counter and rotated 
with it. 

The observed 6 2293-y 691 correlation is shown in Fig. 3(a). The net 
coincidence count plotted for each angular setting is the sum of the counts 
for that position and for its equivalent position in the 180°-270° quadrant, 
e.g., 120° and 240°. The coefficients for this correlation, obtained by the least 
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squares method and corrected for finite counter solid angles, are A2/Ao 
= +0.404+0.012 and A,/Ay = —0.015+0.019. These values were obtained 
with a B-energy window of 1700 to 2300 kev. As a check on the experimental 
procedure separate analyses have been made of the data obtained in the two 
quadrants. The counter-corrected results are, for 90°-180°: A2/Ay = +0.408 
+0.018, A,/Ao = —0.019+0.029; and for 180°-270°:A:/Ay = +0.393 
+0.018, A,/Ao = —0.030+0.029. In these calculations no allowance has been 
made for the effect of back-scattering in the source material and backing or 
for air scattering of the 6-rays between the source and the counter. The 
angular response of the $-counter, measured with a narrowly collimated 
(< 2°) beam of 8-rays from the Ce!4—Pr'™ source and with identical bias set- 
tings, showed no evidence of air scattering. 

The coefficients, A »/Ao, one would expect under our experimental conditions 
have been evaluated for the 0— (8) 2+(y) 0+ and 1—(8) 2+(y) 0+ sequences 
for all first forbidden matrix element combinations listed by Biedenharn and 
Rose (1953). Allowance was made in each probable case for the 8-energy 
window used in the experiment (1700-2300 kev). 

For the 0—(6) 2+(y7) 0+ case, first forbidden selection rules permit only 
tensor interaction with pure > ;,|B4,|2, or axial vector interaction with pure 
‘> «;|Bi,;|?. The theoretical results for both of these are A2/Ay = +0.460. The 
theoretical coefficient A4/Ao is zero for all first forbidden 8-ray correlations. 
The discrepancy between the observed value of A2/Ao (0.404+0.012) and 
the value predicted for the 0—(8)2+(7)0+ sequence may result from an 
attenuation of the directional correlation as observed experimentally due to 
(a) back-scattering from the relatively thick source mounting, (6) undetected 
large angle air scattering of the 6 particles, (c) nuclear precession in the 691 
kev state (which is dependent on the solid state conditions of the source). 
Furthermore, the theoretical coefficients listed by Biedenharn and Rose 
(1953) are strictly valid only for low Z. It has been pointed out by Alaga et al. 
(19576) that relativistic nuclear corrections to the theory may be several 
per cent. It should also be noted that the present theory takes no account of 
the finite distribution of nuclear charge. 

For the spin-sequence 1—(8)2+(7)0+, the coefficients A2/A » which agree 
with experiment are realized for transitions involving mixed 8-matrix elements. 
For tensor interaction these are > B‘, with {Ba and also >B4, with fpexr 
and for axial vector Bj, with [éXr. Two of the cases are illustrated in Fig. 
4. In all cases it is possible to get agreement with the observed experimental 
correlation anisotropy, if the B,; contribution to the decay intensity is ~92% 
or ~40% and the interfering contribution is ~8% or ~60%. 

The 40% B,; possibility is excluded by the shape factor analysis discussed 
in Section III. The effect of attenuation is to reduce the B,; contribution for 
the other possibility below 92%. Our shape factor results do not exclude such 
an admixture. There are, however, no known cases where the B;; (L = 2) 
contributions compete favorably with Z = 1 first forbidden contributions. 
Our reasons for choosing 0— as the correct assignment for the ground state 
of Pr" are discussed in detail in Section VI. 
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(e) B 803-y 2181: Results and Discussion 

The 6 803-y 2181 directional correlation was measured in the same manner 
as that for 8 2293-y 691. The energy calibration of the 8-counter was de- 
termined with the 625 kev K conversion peak from an open Cs" source. For 
this correlation a B-energy window of 550-785 kev was used. The counter 
angular response for this energy window was measured using a strong Cs! 
source and the same collimator as used for the 8 2293 tests. The response 
curve was broader than that for the 8 2293 energy window, and revealed 
some large angle scattering. 

The results for this 8 803—y 2181 directional correlation are displayed in 
Fig. 3 (0). The coefficients deduced from these results, 42/49 = +0.02+0.02 
and A,/Ay = —0.05+0.04, are corrected for counter solid angles but not 
for the effects of large angle scattering. If this 8-transition is first forbidden, 
small but significant anisotropies are expected for several of the possibilities. 
For a spin-sequence 0—(8)1+(7)0+ the experimental results exclude pure 
scalar interaction, and pure vector unless the nuclear matrix element fal? 
dominates. Both tensor and axial vector predictions are consistent with the 
experiment. For the spin-sequence 1 — (8)1+(7)0+ the anisotropies predicted 
for all interactions are small and lie within the experimental limits. No 
anisotropy is expected for allowed 8-spectra, e.g. for O—(8)1—(y)0+ or for 
1—(8)1—(y)0+. 

This experiment provides no conclusive evidence for spin and parity assign- 
ments but is consistent with the decay scheme proposed in Fig. 1. 


III, COINCIDENCE BETA SPECTRA 


These experiments were carried out using the double lens, point focus 
spectrometer shown in the inset of Fig. 5. It was adjusted to a resolution of 
2.7% as shown by the width of the K 133 conversion line shown in Fig. 5 
(inset). The source was prepared by subliming, in vacuum, Ce™ oxide from a 
tungsten filament at ~2000°C onto an aluminum foil backing of ~800 
ug/cm*. The quality of the source is evident from the profile of the K 134 
conversion line which has an energy of ~90 kev. It does not have the low 
energy tail which is characteristic of a thick source. The @ scintillation counter 
and the 2 in.X2 in. Nal (Tl) ¥ scintillation counter were connected to a 
fast-slow (279 ~~ 210-8 second) coincidence circuit (Bell, Graham, and 
Petch 1952). The 8-counter sensitivity was deduced from integral pulse 
height bias curves taken at several spectrometer current settings. At the bias 
setting used the variation of sensitivity was small (< 3%) over the energy 
range under study. The y-counter was pulled back about 1 in. from the source 
and was covered with 3/16in. Pb and 1/16in. Al in order to shield it from 
8-rays and to attenuate the intense low energy y-rays from the Ce! decay. 

The coincidence data were automatically recorded using a solenoid-operated 
electric typewriter. A 10-step repeating cycle changed the current setting of 
the spectrometer in equal increments at 10-minute intervals. The random 
coincidence rates for these experiments were negligible and were occasionally 
checked manually. The y-counter circuit was gain stabilized using an Hg? 
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Fic. 5. Analysis of the 8-ray spectrum of Pr'4 obtained by coincidence methods using 
the double lens spectrometer shown in the inset. The Fermi plot of the component in coin- 
cidence with y 691 has a first forbidden unique shape while that in coincidence with y 2181 
is allowed. The error bars indicate the standard deviations of the counts. The subcomponents 
are plotted to arbitrary scales. The coincidence data are plotted without applying the direc- 
tional correlation corrections discussed in the text. 


y-source, in the manner discussed in Section I], in order to avoid errors due 
to drift and the imperfectly shielded spectrometer magnetic field. 

The 8-spectrometer was calibrated with a ~200 uc source of thorium 
active deposit using the F, J, L, 7, and X lines whose energies are accurately 
known (Siegbahn and Gerholm 1955). 

Fig. 5 shows allowed Kurie plots of the total spectrum and the two sub- 
components. The 8 2293 kev spectrum component was observed in coincidence 
with the 7 691 photopeak and the 8 803 kev spectrum component with a y 
pulse height window including both the y 2181 and the y 1490 photopeaks. 
Corrections have been made for source decay, background, and the variation 
of the B-counter sensitivity with energy. No allowance has been made for back- 
scattering. 
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The allowed Kurie plot of the 8 2293 component (uncorrected set of points 
in Fig. 5) shows the “a” shape characteristic of a first forbidden ‘‘unique”’ 
spectrum. The matrix element > ;;|B;,;|? (LZ = 2) therefore makes the dominant 
contribution to this transition. The 6 2293—y 691 directional correlation result 
of Section II suggests that there is an ~8 or 60% L = 1 8-contribution to this 
transition. In order to determine the Z = 1 6-contribution permitted by the 
experimental spectrum shape, a detailed analysis of these data was performed 
after they had been corrected for back-scattering. This correction was based 
on the discrepancy between the ground state spectrum shape obtained with 
this spectrometer and that found using the double lens, ring focus spectrometer 
described in Section V. The ratio between these ground state spectrum shapes, 
scaled in energy by the end-point energy ratio (2293 to 2984) was applied to 
the 6 2293 coincidence data as the back-scattering correction. The maximum 
correction here is < 6%. 

The corrected coincidence 8-spectrum data do not represent the true shape 
of this component because the anisotropy of the 8-y directional correlation 
varies with B-energy; the exact manner depends on the matrix element ad- 
mixture assumed. Spectrum shape factors (see Section V) were deduced from 
the data for two assumptions, (a) pure >>/B,,|? contribution, and (6) 86% 
> |B,,|? (L = 2) with 14% (ZL = 1) allowed shape contribution. The angular 
correlation effects were removed from the data using the correlation coefficients 
measured in Section II, the angular apertures appropriate to the coincidence 
arrangement, and the theoretical energy dependence of the correlation co- 
efficients for the admixture assumed. The resulting data were then treated 
with the appropriate theoretical shape factor correction and plotted versus 
energy. The shape factor for the correct admixture will be a horizontal straight 
line. The least squares fitted slopes to our data on the two assumptions are: 
for (a) —(3.4+4.6)% per Mev, and for (6) +(3.7+4.6)% per Mev. Since 
zero slope corresponds to an exact fit we can only conclude that the 6 2293 
component contains more than ~75% > B;;. 

The Kurie plot of the 8 803 coincidence data is consistent with an allowed 
spectrum shape. No corrections have been applied to the data shown in Fig. 5. 
The slight drop in counter efficiency at low energies and the rise at low energies 
due to back-scattering, in any case, tend to cancel each other. An allowed 
shape for this spectrum is consistent with the correlation result of Section II. 





IV. BRANCHING RATIOS 


In order to arrive at an accurate experimental result for the shape of the 
main ground state 8-spectrum from Pr! it is necessary to know the intensities 
of the two B-subcomponents as well as their shapes. As pointed out in Section 
I, earlier experimenters do not agree on these intensities and it seemed worth 
while to make a new determination. The method of measurement used, that 
most readily carried out with the apparatus at hand, was to determine the 
relative intensities of the 134 (Ce decay), 691, 1490, and 2293 kev y-rays 
with a calibrated scintillation counter and to obtain the intensity of the K 134 
conversion electron line relative to the Pr™ total 8-spectrum intensity using 





GRAHAM ET AL.: AXIAL VECTOR INTERACTION 1097 


a magnetic lens spectrometer. These results, together with the conversion 
coefficient for y 134, determine the absolute intensities of the y-rays occurring 
in the Pr'* decay and hence the 6 branching ratios. 

The experimental values for the y 134 K conversion coefficient, ax, are of 
low accuracy. These measurements, however, together with those of the A/L 
ratios indicate an M1 multipolarity for this y-ray. Pullman and Axel (1956) 
find ax = 0.7140.25 and K/L = 7.140.5. Emmerich et al. (1954) find 
K/L = 8, Porter and Cook (1952) find K/L = ~8.8, and our experimental 
value is K/L+M = 7.2+1.0. These results are consistent with pure 1/1 but 
do not exclude E2 admixture up to ~15%. The limits of error are also too 
large to rule out completely an admixture of ~90% El with ~10% M2. 
Cork et al. (1954) find K:Z1:L2:L3: = 53:10:0: ~1. The interpretation of 
this result is not clear, but theoretically one expects L2:L3 > 1 for all M1 
radiation. Sliv’s M, conversion coefficient ag = 0.49 (Sliv and Band 1956) 
was chosen in preference to that of Rose since Sliv’s values are in better 
agreement with experiment (cf. Ewan et a/. 1957). While their results suggest 
that Sliv’s values of 1/1 conversion coefficients are some 10% high, such a 
discrepancy has no serious effect on the interpretation of our spectral shape 
results. 

The relative y-ray intensities were obtained using a packaged 4 in. X4 in. 
diam. NalI(T1) scintillation counter mounted on a 5 in. Dumont phototube. 
The scintillation spectra obtained with the ~50 uc source of Ce'*—Pr™ are 
shown in Fig. 6 (a) and (0). The variation with energy of the total absorption 
peak efficiency was determined using sources of Na”, Na*, Zn ®, and Cs! under 
the same conditions. After allowing for absorption in the glass source con- 
tainer and Nal(TI) crystal container and applying escape peak corrections 
(cf. Bell 1955) we obtain the relative quantum intensities listed in Table 
II. The absolute intensities listed in Table I] are obtained according to 


TABLE II 
QUANTUM INTENSITIES 


~ 


Quantum intensity per Pr! 6-disintegration 
Relative intensities in % 


Photon This Alburgerand Kreger and This Emmerich Kreger and Pullman 
energy work*  Krauschaar? Cook* work’ et al. Cook? and Axel* 


134 6.75+1.0 8.6 10.6 +1.1 15.4 8.3 
691 1.0 1.0 1.0 1.57+0.28 1.25+0.2 1.79 
1490 0.17+0.02 0.4 0.20 0.2640.06 0.50 (c) 0.35 
2181 0.50+0.05 i 0.51 0.78+0.16 1.37 (c) 0.91 


“By scintillation spectrometer method. 
’By external conversion method. 
Values deduced using Alburger and Krauschaar'’s relative intensities. 

4Assuming the 134 kev level is fed in 22% of the disintegrations; cf. Emmerich et al. (1954). 
*From disintegration scheme considerations, 

/ Assuming y 134aK = 0.49. 


the relation: J,/Ig = (1y/TIy134)(1/ari34)(K134/[3). The value of AK134/ 
Ig(Pr'4) = 0.052, which we use here, was obtained using the lens spectrometer 
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Fic. 6. Scintillation spectra of the y-radiations from Ce! and Pr!“ in equilibrium taken 
with a 4 in. X4 in. diam. Nal(TI) packaged crystal mounted on a Dumont-type 6364 photo- 
multiplier. The data were recorded with an Atomic Instruments 20-channel pulse height 
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analyzer. The background has been subtracted. Typical standard deviations are indicated 
by the error bars. Calibration curves are shown for comparison in broken lines. Plot a shows the 
low energy region using the geometry shown in the inset. Plot b shows the high energy region. 
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shown in Fig. 5. Porter and Cook’s (1952) value is 0.058. For comparison 
some relative and absolute intensities of other workers are listed as well. 

The 8 branching ratios deduced from these absolute y intensity measure- 
ments for the Pr'4 decay scheme as shown in Fig. 1 are 6 803 (1.04+0.22)% 
and 8 2293 (1.3140.34)%, implying a (97.65+0.44)% intensity for the 
8 2984 ground state transition. 

V. GROUND STATE 86-SPECTRUM SHAPE 

In order to study the total 6-spectrum with as much freedom from scattering 

as possible, the double lens spectrometer shown in Fig. 7 was used. It was 


adjusted for a transmission of 0.3% and a resolution of 0.5%. The vacuum 
chamber and spiral and ring-focus baffle assembly have been carefully designed 
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DOUBLE LENS MAGNETIC SPECTROMETER 


Fic. 7. Cross-sectional diagram of the double lens spectrometer used for the detailed 
examination of the Pr! 8-spectrum shape. Modifications for these experiments which differ 
from the diagram are the use of (1) a smaller source mount ring and (2) a larger light pipe. 


to prevent the detection of scattered electrons. A detailed study carried 
out by Elliott et a/. (1954) with a strong P® source has shown that at a setting 
of 5 line widths beyond the spectrum end-point the scattered rate is « 10-° 
of the peak rate on the spectrum and at higher momentum settings could not 
be detected. 

The source used for this study was made by subliming Ce! oxide from a 
tungsten ribbon filament at ~2000° C onto 200 yg/cm? aluminum backing. 
The resulting deposit of 2 wc was just visible and had an estimated surface 
density of < 100 ug/cm*. The foil was supported on a 7/8 in. I.D., 0.060 in. 
thick brass ring which is smaller than the one illustrated in Fig. 7. As shown, 
the region behind the source is free of scattering materials but it is possible 
that some of the 8-rays scattered from the inside edges of the brass support 
ring could reach the counter. 

The detector is a 16 mm diam. by 2 mm thick disk of anthracene set in the 
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end of a 2-ft long lucite pipe which leads to a Dumont 6292 photomultiplier. 
The anthracene crystal is covered with a 200 yg/cm? aluminum foil light 
reflector. The photomultiplier was carefully selected for efficiency and low 
noise and was operated at room temperature. The discriminator bias was set 
above the noise and after pulse region and bias curves were run at many 
spectrometer current settings. The detector efficiency at the bias setting chosen 
for this spectrum study was deduced from these bias curves in a manner 
similar to that described by Porter et a/. (1956). 

In order to obtain the shape of the ground state spectrum the contributions 
due to the 803 kev and 2293 kev 8-components must be subtracted from the 
total spectrum. The intensities of the components used in this subtraction 
were those deduced from the y-studies in Section IV. It was assumed that the 
803 kev component had a theoretically allowed shape and that the 2293 
kev component had a pure first forbidden unique >>| B;;|? shape as suggested by 
the coincidence spectrum studies described in Section III. The spectrum 
shape was analyzed using / values obtained from the compilation of Dzelepov 
and Zerianova (1956). These F values incorporate the factor Zo, which is 
usually evaluated separately, and allow for uniform nuclear charge distribution 
and for screening by the orbital electrons. The F values used in this analysis 
were those appropriate to a nuclear radius p = 1.21078 A! cm and were 
obtained from Dzelepov and |Zerianova’s tabulated values, p = 1.5X10-" 
1" em, by scaling in the prescribed manner (for this case F corresponds to 
FyLo where Fy and Ly have the definitions given by Rose et al. (1953)). The 
ratio of F(po = 1.2) to FoLo(po = 1.41) varies less than 0.5% over the energy 
range of interest, i.e., 300 to 3000 kev. 

To compare experiment with theory we divide the observed momentum 
distribution (.V = counting rate/Bp) by the theoretically allowed distribution 
Fp*q?. This ratio is plotted as a function of 6 kinetic energy in Fig. 8 and is 
hereinafter referred to as a shape factor. The shape near the high energy 
end-point is very sensitive to the end-point assumed because of the q? factor. 
We have chosen the end-point value (2984 kev) to minimize the fluctuation 
above 2500 kev. The effect of different end-point assumptions is illustrated 
in Fig. 8 by broken lines. 

As a check on possible systematic errors in the procedures used for the 
Pr work the P® 14+ — 0+, 8-ray spectrum was examined in the same 
apparatus and was analyzed using the same procedure. Two sources were 
used: one was made by evaporating P®CI solution on 800 ug/cm? Al backing 
and the second was prepared by subliming P®CI from a tungsten filament 
onto ~200ug/cm? foil. The brass support rings were identical with that of 
the Pr'#* source. Consistent shape factor plots were obtained with these two 
sources. That for the sublimed source is shown in Fig. 9 together with the 
data of Porter et al. (1956). Our results indicate slightly more low energy 
electrons than found by Porter et a/. The reason for the discrepancy is not 
fully understood at present. Even if the discrepancy is entirely due to 
systematic errors in our experiment the interpretation of the ground state 
spectrum (see next section) in terms of axial vector interaction is not affected. 
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Fic. 8. Experimental shape factor for the ground state 8-spectrum of Pr! after sub- 
tracting the 1.3% 8 2293 and the 1.0% 8 803 components from the total spectrum. An allowed 
shape is a horizontal straight line. Errors which can arise from the wrong end-point assumption 
are illustrated by the broken lines at the right. Theoretical shape factors shown are for pure 
tensor, pure pseudoscalar, and two possibilities for pure axial vector interaction assuming a 
O— to 0+ spin-sequence. Typical standard deviations are shown by the error bars. 
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Fic.9. Experimental shape factor for P® obtained in the spectrometer shown in Fig. 7 using 
the same procedure as for Pr, The results of Porter et al. (1957), divided by Lo, are plotted 
as open circles for comparison. 

VI. DISCUSSION 

The spin assignment of Pr' has not been established with complete cer- 

tainty. While the assignment 0— has been favored by most authors to date, it 
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is clear from the coincidence studies above that a 1— assignment should also 
be considered. As has been discussed in Section III this assignment would 
require that the 6 2293 component have > 75% ¥)|B,,|? (L = 2). This tran- 
sition has a log fot = 7.9 and log ft = 8.0 which is typical of a first forbidden, 
pure L = 2 transition where |J;—J2| = 2 (cf. Wu 1955). No cases of L = 2 
transitions have been identified so far where the spin difference involved is 
only 1, i.e., |J;|—|J2| = +1, although such contributions are always possible 
except, of course, where one of the spins is zero. A pertinent example is the 
884 kev 8-component (1— to 2+) which occurs in the decay of Tm!” (Graham 
et al. 1952). This component has an unusually large ft value of 9.3 yet it 
has an allowed shape and is not consistent with any significant LZ = 2 contri- 
bution. This agrees with the theoretical expectation that the L = 2 contri- 
butions are inhibited by a factor of (aZ/2p)? ~ 10? relative to competing 
L = 1 contributions. For these reasons it is difficult to accept a 1— assignment 
for the ground state of Pr; in the discussion that follows we assume that the 
correct one is 0—. 

Only three interactions, tensor (7), axial vector (4), and pseudoscalar (P), 
may contribute to a first forbidden 0— to 0+ 8-transition. The shape factor 
formulae for the pure interaction Siz, Si4, Sip and two of the interference 
terms are taken from Rose and Osborn (1954). Lee-Whiting (1957) has 
deduced the interference formula S;74 for us. These formulae are: 
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where g = neutrino momentum, M = nucleon mass in units of the electron 
rest mass ~ 1836. The potential U = W—V.—q where V, = the potential 
of the field in which the 8-particle moves (units of moc?). We take V, = Ze/p 
where p = nuclear radius (units of h/moc = 3.84X107-'! cm). We assume, 
with Rose and Osborn (1954), that fpe.r which is the correct nuclear matrix 
element for tensor interaction ~ fe.r* 

If two interactions contribute to the transition the shape factor S for the 
spectrum is obtained by adding the pertinent formulae, e.g., for a 7 P mixture 
S = Sirpt+Sirp+Sip. The primed coupling constants denote the non-parity 
conserving contributions. The functions Lo, Mo, Vo, etc. have been defined 
by Greuling (1942). 

In our analyses the theoretical shape factors appropriate to the first forbidden 
0— — 0+ axial vector and tensor interactions, i.e., Siz, Si4, and Sip4, were 
evaluated using the accurate values of Mo, No, Po, Qo, and Ro calculated by 
Zerianova (1956) fora Z = 60 nucleus of radius p = 1.2X10-' A!cm. These 
include allowance for a uniform nuclear charge distribution. Her evaluation 
of these functions assumes that the ‘‘allowed”’ factor Lo is incorporated into 
the Fermi function F and Lp is therefore replaced by unity in evaluating the 
shape factor formulae. In addition the shape factors Siz, Sia, Sip, Sirp, and 
Siap have been evaluated using values of Lo, Mo, No, Po, Qo, and Ro from 
Rose’s tabulation where p = 1.41 A‘ 10-% cm. These shape factor values 
must be divided by Lo since this factor is already taken account of in F. Except 
for the regions of destructive interference, which are discussed below, the 
differences in the normalized shape factors computed from these two sets 
of functions is typically < 1%. 

The experimental points are in reasonable agreement with the allowed 
theoretical shape factor which appears as a horizontal line in Fig. 8. The 
theoretical shape factor for a O— to 0+ spectrum assuming pure tensor 
interaction has a strong positive slope, and is in obvious disagreement with 
the experimental results as evident in Fig. 8. The pure pseudoscalar shape 
factor has a similar slope and can also be rejected. 

Only a qualitative fit to the data can be obtained with a tensor—pseudoscalar 
combination using Rose’s and Osborn’s (1954) theory. This occurs when 
destructive interference is approached which requires a large pseudoscalar 
coupling constant, i.e., |Cp|/|C7|~10%. To illustrate this we show theoretical 

Dr. M. E. Rose has recently advised us that the factor V3 which appears in the final term 


of equation (6) resulted from an arithmetic error and should be unity. This change will alter 
slightly the pseudoscalar limits given in Table III for \ = 14 and A = 200, yes cases. 
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shape factors for 7 P admixtures in Fig. 10 for several ratios of Cp/Cr. All 
have been arbitrarily normalized to unity at 2100 kev. Our evaluation assumes 
time reversal invariance (C’s all real) and also 


Cp/Cp = +Cr/Cr, 


i.e., the interference term present, as expected on the two-component neutrino 
theory. In evaluating V., p was taken equal to the same nuclear radius 
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Fic. 10. Attempts to fit the Pr 0— to 0+ 8-spectrum shape factor with pseudoscalar 
tensor admixtures near the region of destructive interference. The formulation of Rose and 
Osborn (1954) was used and the tables of Rose et al. (1953) were used taking p = 1.41 A¥ 
10-3 cm. The theoretical shape factors are all arbitrarily normalized to unity at 2100 kev. 


assumed in the Rose ef al. (1953) tabulation, i.e., p = 4a A!#® = 1.41 A!’ 
10-" cm. The theoretical factors were divided by Lo to conform to the treat- 
ment used in Fig. 8 so that again an allowed shape is a horizontal straight line. 
It is apparent that no quantitative fit can be achieved. Since destructive 
interference is required, ~98%, one cannot calculate these factors from Rose’s 
4-figure tables to an accuracy of better than ~2%. 

Two nuclear matrix elements are present in the theoretical shape factor 
for pure axial vector interaction. The shape factor for | f-é.1r|? and for lSvsl? alone 
are both shown in Fig. 8 evaluated using Zerianova’s functions (1956). The 
shape of the interference term is midway between these two. The results 
seem to favor (fo.r|? but with the present accuracy either matrix element 
could be dominant. Admixtures are also acceptable both without interference 
and with interference provided that it is constructive, i.e., \ = i(f-vs)/ 
(f6.r) > 0. Pearson (1958) has shown that J is positive if one assumes a j—j 
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coupling single particle model with harmonic oscillator potentials of approxi- 
mately the same strength for parent and daughter states. Destructive inter- 
ference in the range \ = —3 to —40 gives rise to shape factor slopes which 
are inconsistent with the experimental results as shown in Fig. 11. 
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Fic. 11. Attempts to fit the Pr O— to 0+ £6-spectrum shape factor with pure axial 
vector for negative values of X = i{fys}/{f¢é.r}, i.e., destructive interference. The shapes for 
positive values of \ lie in the narrow region between \ = Oand A = © and are consistent with 
the experimental results. 


In Section V we noted that the shape factor for P® obtained in our equip- 
ment differed slightly from that of Porter et al. (1957). If we assume that the 
discrepancy is entirely due to a systematic error in our experiment it does not 
alter the conclusion regarding axial vector interaction except perhaps to 
favor dominance of the |f{é.r|? contribution. 

It should be pointed out that agreement with axial vector interaction is not 
sensitive to nuclear radius or finite distribution of nuclear charge. To illustrate 
this point we have also evaluated the theoretical shape factor appropriate 
to the conditions assumed in Fig. 8 using the Rose et a/. functions (1953). This 
is shown in Fig. 8 by the broken line Sialfo.r|? (Rose). As is evident, the 
experimental results are again consistent with this factor. 

It is of interest to see what admixture of tensor or pseudoscalar interaction 
with axial vector interaction is permitted by the experimental results. The 
slope of the axial vector shape factor varies violently in the region of de- 
structive interference of ifs with fo. r and it is difficult to make generalizations 
for this region. There are 7A or PA admixtures which are consistent with 
the experimental data for certain ranges of the ratio } = ifys four. The 
admixture required will be a sensitive function of \. The situation for con- 
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structive interference between ifys and fé.r is much clearer. Since the experi- 
mental data are consistent with any such axial vector shape factor one can 
then explore the amount of admixture of tensor or pseudoscalar interaction 
which is still consistent with experiment. We have carried out computations 
for a variety of values of \ and summarize our conclusions in Table III which 
lists the approximate limits of 7A and PA admixtures for three values of X. 


TABLE III 
Limits ON J—-A and P-A ADMIXTURES FOR THE 0— TO 0+ Pr! GROUND STATE SPECTRUM? 


Axial vector Tensor Pseudoscalar 
Interference 


A = {tfys/fé.r}> fyscontribution® term* Cr/Ca = Sir/Sia®@ = Cp/Ca  Sip/Sia4 


0 0 Yes <0.1 < .007 <15 <.009 
No <0.3 <.06 <S <.10 

14 50% Yes <0.2 <.007 <37 <.014 
No <0.8 <i <120 <.14 

200 99.5% Yes <2 <.012 <180 <.005 
No <7 <.15 <1000 <.20 





®Experimental data analyzed with the formalism and the functions of Rose et al. (1954) using p = 1.41 A”3 19718 
cm. Invariance under time reversal is assumed. 

’We here consider constructive interference only for reasons discussed in text. 

©f | frsl2Lo} /{ [fs]? Lot! [6.x]? (Mo— 3a Not hq? Lo)} at Eg ~1100 kev. 

4Numerical contributions evaluated at Eg ~1100 kev using equations 1, 2, and 3. 

*Siar7 for tensor and S14Pp for pseudoscalar admixtures. 


The limits were arrived at graphically. A typical example is shown in Fig. 12 
for the case \ = 0. It is clear from the table that no substantial contribution 
from either tensor or pseudoscalar interaction is permissible with axial vector 
with any positive value of X. 

Laubitz (1956) obtained an experimental shape factor which differs ap- 
preciably from that shown in Fig. 8. In contrast to our experimental shape 
factor his drops ~20% in the energy region 800 to 2800 kev. The drop above 
2200 kev would appear to have two causes. Laubitz did not apply a correction 
for the 6% resolution of his instrument which is important for his last three 
points. This may have influenced him in choosing too high an end-point energy 
value which in turn influences the values of neutrino momentum used in his 
analysis. If one makes the appropriate resolution correction and assumes an 
end-point energy ~20 kev lower the pronounced droop at the high energy 
end is largely eliminated. The rise at low energy may have been caused by 
back-scattering from the rear of the source mount. This is suggested by the 
ground state spectrum derived from our 3% lens spectrometer results shown 
in Fig. 5. This gives an experimental shape factor almost identical with that of 
Laubitz’ below 2000 kev. For these reasons we feel that his analysis of his 
experimental data does not provide good evidence for a strong pseudoscalar 
contribution in this 0O— to 0+ spectrum. 

Zerianova (1956) has evaluated more precise 7P shape factors in the 
region of destructive interference for Pr™ in attempting to fit the experi- 
mental results of Emmerich, Auth, and Kurbatov (1954). She has computed 
the pertinent electronic wave functions to 6-figure accuracy taking into 





CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 





T T T T “I 


N/p* q? Lo Fo—> 





500 1000 1500 2000 2500 3000 
B KINETIC ENERGY (kev) 


Fic. 12. Attempts to fit the Pr’ 0— to 0+ 8-spectrum shape factor with axial vector, 
= 0, and a small admixture of tensor interaction, both with the interference} terms, 


Cr/C4h = +Cr/Ca, 
and without the interference terms, 
Cr/C4 = —Cr/Ca. 


The theoretical curves were calculated using Zerianova's (1956) functions with p = 1.2 A‘ 
10-5 cm. The shape factors are arbitrarily normalized at 2200 kev. 


account screening and uniform distribution of nuclear charge, assuming a 
nuclear radius value p = 1.2 A!“ 10-"% cm. Her conclusion is similar, i.e. one 
cannot obtain quantitative agreement with an experimental shape factor that 
is essentially allowed using the pseudoscalar formulae of Rose and Osborn 
(1954). 

The value of log ft = 6.5 for the Pr'* ground state transition is typical of 
those of other first forbidden transitions. In that region of PT destructive 
interference which qualitatively fits the experimental data, e.g., for Cp/Cr 
= 107 (see Fig. 10), the presence of pseudoscalar interaction reduces the 
total decay probability to ~3% of that expected for the pure tensor inter- 
action. Hence the tensor part has a log ft ~ 5. Zerianova (1956) deduces that 
the cancellation is even more complete, i.e., S ~0.2% Siz, implying that the 
log ft for the tensor part is ~4. Log ft values as low as 5 are found only for 
highly favored first forbidden transitions which occur near closed shells (cf. Wu 
1955). Thus there is difficulty in accounting for the transition probability as 
well as for the spectrum shape with the Rose and Osborn pseudoscalar theory. 

There is lack of agreement regarding the validity of the Rose and Osborn 
theory for pseudoscalar interaction. A modification has been proposed by 
Alaga, Sips, and Tadic (1957(a)), who introduce a term into the Hamiltonian 
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for the B-interaction to take account of the influence of nuclear forces. In this 
way they derive new formulae for S;p and Sip. This gives an additional para- 
meter T' as well as Cp/Cr. So far they have only attempted to fit Laubitz’ 
experimental shape factor. Their fit with Cp/Cp ~10-15 requires a I of 
~30-40. These parameter values imply that the pseudoscalar contribution 
dominates the tensor contribution, Sp > Sz, whereas in the approximate 
fits obtained using Rose’s formulation the reverse is true, i.e., Sp < Sr. 

If one questions the method followed by Rose and Osborn in evaluating 
the nuclear moment SBys on which the pseudoscalar coupling acts one may 
return to the earlier treatment of Konopinski and Uhlenbeck (1941). This 
gives the following formulae for T and P shape factors (cf. Konopinski 1955). 


J Brs|"Lo, 





(|Cp|*+|Cp|*) 


(7) Sip 





(8) Sipp = 2Re(CrCp+CtCr )i(fBrys)([B6. 4) "{No+4qLo}, 








(\Cr|?+ |Cr|”) JBe : r|°{ Mot+3qNot+sq' Lo}. 


(9) Sir 





In order to fit our experimental data with these factors we require a domi- 
nating contribution of pseudoscalar interaction which has the allowed (Lo) 
shape. Because of the uncertainty regarding the relative magnitudes of the 
two nuclear matrix elements it is difficult to derive good limits on the Cp/Cr 
coupling constant ratio. If |{8ys|~|f@é.r| our data requires |Cp/Cr| > 10. 
In the light of Rose and Osborn’s arguments it seems more probable that 
\{Bys| « |f@6.r| which would require that |Cp/Cz| > 10 in order to account 
for our experimental data on the Konopinski Uhlenbeck formulation. In this 
event, however, one should also take into account the contributions from the 
smaller terms considered by Rose and Osborn. The present uncertainties 
regarding nuclear wave functions preclude a unique choice of the formulation 
for the pseudoscalar theory. 

This uncertainty in treatment of our data in terms of PT admixtures is in 
marked contrast to its simple explanation in terms of pure axial vector inter- 
action. 


VII. CONCLUSIONS 


The agreement between the experimental shape factor of a 0— to 0+ 
6-spectrum and the theoretical shape for axial vector interaction is not peculiar 
to Pr. The available experimental evidence favors a 0— assignment for 
the TI ground state. The TI?” 8-ray spectrum (0— to 0+) is unmasked by 
low energy components and has an allowed shape (Alburger and Friedlander 
1951) consistent with axial vector predictions. The 0— assignment for the 
ground state of 27-hour Ho’ has been established in the same experimental 
manner as described in this paper for Pr'* (Graham, Wolfson, and Clark 
1955). The 1854 kev 0O— to 0+ ground state 8-spectrum from Ho" is masked 
by an almost equally intense (48+4)% 1771 kev component which has a 
first forbidden ‘‘unique’’ B;; shape. The shape of the 1854 kev spectrum is not 
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consistent with pure tensor interaction, but is, within the rather large errors, 
consistent with an axial vector shape (Wolfson and Graham 1955). The 
agreement in spectrum shape of these experimental examples is more satis- 
factorily explained by pure axial vector interaction than by invoking in- 
terfering 7 P admixtures which are very sensitive to the pseudoscalar coupling 
constant and which, in any case, only give qualitative fits. 

This evidence from spectrum shapes for axial vector interaction in Gamow 
Teller decay is consistent with other recent experiments on the nature of the 
B-decay interaction. The electron recoil angular correlation results on the 
positron emitters A** and Ne reported by Herrmannsfeldt et al. (1957) 
are most simply accounted for by a VA combination of interactions. The 
circular polarization and resonant scattering results of Goldhaber et al. (1957) 
on the y-rays from Sm! following K capture in Eu'® indicate a negative 
helicity for the neutrinos emitted in this decay. On the assumption that this 
neutrino helicity is the same as that for positron emission they conclude that 
this Gamow Teller decay is axial vector. The evidence to date for negatron 
emitters, excluding the He® study, comes from studies of the neutron decay.* 
The electron recoil results of Robson (1955 and 1958) on the neutron decay are 
consistent with the combination ST or VA but not with SA or VT. The 
measurements of Burgy ef al. (1958) on the 8 and proton recoil asymmetries 
in the decay of polarized neutrons favor the vector —axial vector combination 
for this decay. 

The 8 2293~y 691 directional correlation results and the coincidence spectrum 
shape of the 6 2293 component strongly favor a O— spin assignment for 
the Pr'* ground state. The experimental shape of the 0— to 0+ Pr'* 6-spec- 
trum is in disagreement with that for the tensor interaction and with those 
for pseudoscalar and pseudoscalar tensor admixtures on the Rose and Osborn 
(1954) pseudoscalar theory. The axial vector interaction provides the only 
satisfactory explanation for this 0Q— to 0+ 8-spectrum shape. This conclusion 
is independent of assumptions regarding the presence or absence of the in- 
terference terms in the shape factor expressions. 
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* Note added in proof: ; 
Herrmannsfeldt et al. (preprint 1958) have recently measured the ion recoil spectrum following 
He® beta decay. Their results ‘‘clearly indicate that the axial vector interaction is dominant.”’ 
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THE FORWARD-SCATTERED RADIO SIGNAL FROM 
AN OVERDENSE METEOR TRAIL! 


P. A. ForsyTH 


ABSTRACT 


A recently presented expression for the forward-scattered signal from an over- 
dense meteor trail is tested by applying it to the signals scattered from a single 
meteor trail. The particular trail was selected carefully in order to permit the 
derivation of each of its parameters in several different ways. The internal con- 
sistency of the derived information is taken as a measure of the validity of the 
theoretical expression. In particular, the electron line density at each of two 
reflection points is derived by three different methods, of which two depend upon 
the new expression, and the other depends only upon well-established theory. 
The resulting agreement is within the experimental error and suggests that the 
new expression is an adequate working relation for use in experimental investi- 


gations. 


I. INTRODUCTION 


The reflection of radio waves by an ionized meteor trail is, admittedly, a 
complex process. However, in its study, a remarkable degree of success has 
been achieved by applying theoretical expressions derived by consideration of 
only two relatively simple models. The two mechanisms concerned should 
apply separately to trails of very low and very high ionization density. The 
first process, which applies to weakly ionized trails, involves the free pene- 
tration of the trail by the radio waves and independent scattering of the 
radiation by the individual electrons. The second process applies when the 
ionization density within the trail is sufficiently great to prevent the radio 
waves from passing through the trail; in such circumstances it is customary 
to treat the trail as if it were a conducting metal cylinder. The formulation of 
the theoretical expressions appropriate to observations made by means of 
radar equipment (the back-scatter case) has been summarized by Kaiser 
(1953) and by Eshleman (1955), and the history of the various treatments has 
been outlined in a recent paper (Hines and Forsyth 1957). 

In observational work it has been customary to assign each meteor trail to 
one of two classes according to the time-variation of amplitude in the reflected 
signal. If, over most of the lifetime of the trail, the low-density expression is 
appropriate then the trail is said to be ‘“‘underdense’”’ whereas if the high- 
density expression applies the trail is termed “‘overdense.’’ In fact, since the 
ionization of the trail is confined initially within a small radius, all but the 
most weakly ionized trails must behave as high-density trails for a short 
time after formation. On the other hand, because of the diffusive spread 
of the ionization, every trail eventually must become sufficiently diffuse to 
permit the free passage of the radio waves and hence must behave as a low- 
density trail near the end of its lifetime. 


‘Manuscript received March 12, 1958. 

*Contribution from the Radio Physics Laboratory of the Defence Research Board, Ottawa. 
The work was performed under project PCC No. D48-28-35-05. 
Can, J. Phys. Vol. 36 (1958) 
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In spite of the obvious limitations of the two separate models and of the 
rather arbitrary bridging which was used frequently to join the two sets of 
expressions in the transition region, the early back-scatter experiments 
achieved a surprising degree of success in the determination of various physical 
properties of the meteor trails and of the upper atmosphere (see, for example, 
the review by Lovell and Clegg 1952, and the more recent work by Lovell 
1954). Later, it was shown by Eshleman (1952) that a more exact treatment 
leads to very similar results but by this time a large body of evidence had 
been accumulated which, by its internal consistency, attested to the validity 
of the expressions. 

In recent years, studies of meteor trails have been carried out using radio 
systems in which the transmitter and receiver are separated by several 
hundreds of kilometers (see, for example, Eshleman and Manning 1954; 
Forsyth and Vogan 1955, 1956; McKinley and McNamara 1956; Hines and 
Vogan 1957). Under these conditions, the reflection process is described as 
forward-scattering in contrast to the back-scatter or radar case. In forward- 
scattering additional complications arise; however, a straightforward extension 
of the simple low-density formulae to include forward-scattering has been given 
by Eshleman (1952). Recently, a similar extension of the high-density formulae 
has been presented by Hines and Forsyth (1957), who also point out that, while 
the approach is oversimplified, it does not seem likely that a more exact 
treatment will lead to a simple result that is suitable for direct comparison 
with observation. 

It has been suggested that since neither of the models is likely to apply 
exactly to a given meteor trail, the resulting expressions should be expected 
to have only a statistical validity, describing adequately the average behavior 
of the signals scattered by a large number of trails. The first attempts to study 
statistically the duration of forward-scattered signals have revealed some 
small discrepancies (Forsyth and Vogan 1956; McKinley and McNamara 
1956) but it is not clear that these are due to a real failure of the theoretical 
expressions since the sampling procedures may have included some obser- 
vational selection. Another approach is to test the formulae by the detailed 
study of the signals scattered from individual trails. In this paper an attempt 
is made to derive as much information as possible from the recordings of 
signals scattered from a single meteor trail. The internal consistency of this 
information is taken as a measure of the validity of the formulae. 


Il. THE EXPERIMENTAL ARRANGEMENT AND SIGNAL SELECTION 

During March, 1956, recordings were made of the meteor signals received 
simultaneously at two stations in the Ottawa area, separated by about six 
kilometers. The transmitter was located at Greenwood, Nova Scotia, 860 km 
from Ottawa, and was operated at a frequency of 39.22 Mc/sec with an output 
power of 100 watts. The details of this circuit have been published previously 
by Vogan and Campbell (1957). The disposition of the receiving stations is 
shown in Fig. 1. Station 1 was located at the Radio Physics Laboratory and 
Station 2 was located at a remote site near the south shore of the Ottawa river. 
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The signals received at the remote station (2) were transmitted to the Radio 
Physics Laboratory (1) by means of a land line and were recorded on the same 
chart as were those received locally. The recordings were made with a chart 
speed of 10 mm/sec. The response of the recording system was such that the 
deflection of the recorded trace was very nearly proportional to the amplitude 
of the received signal. 


4 
STATION 2 y 
&p 


& GREENWOOD N.S. 
— 860 Km—= 





STATION | 





Fic. 1. Showing the disposition of the two receiving stations relative to the direction to the 
transmitter. 
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(STATION 1) 
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Fic. 2. The signals recorded at the two receiving stations at 1415 hours local time, 16 
March, 1956. 


In order to carry out the full analysis it was necessary to look for a pair of 
signals produced by a meteor that had fairly well-defined characteristics. In 
order to remove one of the ambiguities which enter the analysis it was necessary 
that the meteor trail lie in a plane inclined at an angle to the vertical of about 
45 degrees. This plane is more completely specified later. It was necessary 
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also that the signals received at both stations be of reasonably large amplitude 
and free from any significant fading. Finally, it was necessary that both 
signals exhibit clearly those characteristics which indicate that the parent 
meteor trail was intermediate in electron density between the typical ‘‘under- 
dense’ and “‘overdense”’ trails, that is, that the signals show a rounding-off 
of the peak of the amplitude and have, also, an exponential decay. In spite of 
these stringent criteria a suitable event was found in the records after only a 
few hours of recording. The original recording of the signal pair is reproduced 
in Fig. 2. The signals were recorded at approximately 1415 hours local time 
on 16 March, 1956. 


III. DISCUSSION OF THE THEORY 
The method to be used is based upon a fairly simple concept of the way 
in which the signal arises from an overdense meteor trail. This concept may 
be explained with reference to Fig. 3 where an idealized signal is drawn for the 
purpose of illustration. In this diagram the ordinate is the /ogarithm of the 


re 
a 
2 
= 
ee 
a 
= 
< 
© 
o 
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TIME 


Fic. 3. Idealized signal amplitude variation for a meteor trail which has an electron line 
density slightly above the ‘‘transition”’ value. 


signal amplitude and the abscissa is time. The initial sharp rise of the signal 
occurs as the meteor trail, in forming, bridges a principal Fresnel zone of the 
transmitter-receiver system. At formation, the electron density of the trail 
is too great to permit penetration by the radio waves, which are reflected at 
some critical radius. As this radius increases under the action of ambipolar 
diffusion the signal amplitude increases to a maximum value of Ao at time 7 
after the formation of the trail. Thereafter, as a result of further diffusion, 
the critical radius collapses to zero and the radio waves pass freely through the 
volume of the trail. From this point onward the individual electrons scatter 
independently and the signal amplitude decreases exponentially with time as 
the scattering volume increases and coherence of the scattered radiation is 
lost. By extending this decay rate backward in time (as indicated by the 
broken line in Fig. 3) to the initiation time of the signal it is possible to de- 
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termine a hypothetical amplitude A, which would apply if all the electrons 

in the principal Fresnel zone could scatter independently and coherently. 
The expressions appropriate to forward-scattering from an underdense 

meteor trail as given in the earlier paper (Hines and Forsyth 1957) are: 





{3} Ppr= | aS tre seed P7G7rGr < sin‘a a | ‘| 
) ® ~ | 162'RrRe(Rr+Rp)(1—cos 8 sin’¢)IL4m | 2 


‘ ‘ 
A" sec’ ¢ 


’ 


and for the overdense meteor trail :* 


P7GrGr Ty sin’ T poe” 
(3) Bg 8 ne z 
167° RrpRe(Rr+Rpe)(1—cos B sin’) IL4e 4m 


and 


1 
2 


P [oe gn° sec’ y 
4 eS ——— SS 
(4) . 4m 7 
where 


maximum received power, 
transmitted power, 
transmitting and receiving antenna gains (relative to an 
isotropic radiator), 
radio wavelength, 
distances from meteor trail to transmitter and receiver, 
angle between incident electric vector and the scattered ray, 
angle between axis of the meteor trail and the plane that 
contains the transmitter, reflection point, and receiver, 
one half the interior angle formed by the lines along Rp and 
Rr, 

Mo tr X 10-7 henry/meter, 

e = electronic charge, 

m = electronic mass, 

g = line density of free electrons in the meteor trail, 

7 = time constant of the exponential decay of signal amplitude, 

r’ = total duration of ‘“overdense’”’ signal, 

D = ambipolar diffusion coefficient, 

¢ = base of the natural logarithms. 


In the case of the overdense trail the maximum signal is received at a time 
7’/e after the initiation of the signal. 

The geometric factors that enter the expressions (1) and (2) now seem to be 

generally accepted. However, the reasoning that leads to similar factors 


*The expression (3) is equation (18) of the earlier paper. It should be noted that equation 
(17) of that paper should have a bracket at the end of the tirst line and another at the beginning 
of the second line 
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in (3) and (4) has not yet gained such support. Manning (1957) has suggested 
that, in particular, the factor, sec?y, that appears in (4) is quite inappropriate. 

For a signal from a ‘‘transition’’ meteor of the type represented by Fig. 3 
it is possible to measure rt and 7’ = Je. Eliminating D between (2) and (4), 


7’ mm : or 


) “a = 2.4110 


. 
which will provide a correct estimate of the line density only if the same 
geometric factor is used correctly in both (2) and (4). 

A similar situation exists with respect to the amplitudes predicted by (1) 
and (3). Again referring to Fig. 3, it is clear that for a given signal of this 
type, A, should correspond to the amplitude indicated by (1) and Ao to that 
indicated by (3). Combining these two, 


Ap _ aa nd Ae" 
Ay L4m q T 


4 1.33 
(6) q= ao(42) ’ 
“10 


where qo is the transition density obtained by equating (1) and (3), and is 
equal to 0.75 X 10" electrons/meter. This expression will yield the correct value 
of density only if the geometric factors in (1) are also applicable to (3). 
However, in this case, the ratio of the geometric factors in (1) and (3), if not 
equal to unity, would enter (6) raised only to the power 0.66. 

The two expressions (5) and (6) provide independent methods for obtaining 
the line density of the parent meteor trail and, if expressions (1) and (2) are 
accepted, then the accuracy of (3) and (4) can be tested. For this purpose, it 
is desirable to have a third method of analysis that is independent of (3) and 
(4). Such a method is available provided all the quantities necessary for a 
numerical solution of (1) can be determined. This approach is followed in the 
succeeding section where it will be seen that in the course of the analysis it is 
necessary to derive several other parameters of the trail in addition to the 
electron density. 


IV. ANALYSIS 


In Fig. 4 the signal pair under consideration is shown replotted on a loga- 
rithmic amplitude scale. The similarity between each of these signals and the 
idealized signal of Fig. 3 is apparent. In the case of signal 2 (recorded at 
Station 2) there is a slight deviation from the exponential decay law toward 
the end of the signal. This deviation probably is due to a slight distortion 
of the trail in the neighborhood of the corresponding reflection point, but 
it does not adversely affect the analysis. The significant features of each 
signal are easily identified and measured. 

The two signals may be associated with different points on the meteor 
trail. Before proceeding with the electron density analysis it is of interest 
to derive some preliminary information about the trail and the location of 
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the two reflection points. If the meteor trail is projected into a plane that is 
perpendicular to the transmitter-receiver axis, then the condition of specular 


SIGNAL AMPLITUDE ( 4 VOLTS) 


0 2 


I 
TIME (sec ) 


Fic. 4. The signals of Fig. 2 replotted against a logarithmic scale. 


reflection requires that the point in the projected trail that corresponds to 
the reflection point be the point of tangency with a circle drawn about the 
transmitter-—receiver axis. If it is assumed, as is reasonable in the present case, 
that the meteor trail was not formed near the transmitter or near the receiver 
but rather in the general vicinity of the mid-plane between the two (the 
region of validity of the cylindrical approximation of Hines (1955)), then the 
diagram of Fig. 5 (a) may be taken as a representation of the situation as 
projected into the plane that is perpendicular to the axis midway between the 
two transmitter-receiver axes. By reference to Fig. 5() it will be seen that 


Fic. 5. The geometry of the transmission paths illustrating (a) the projection of the 
paths into a vertical north-south plane, and (0) the derivation of relation (7). 
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y, the angle between the vertical and the meteor trail as projected into this 
plane, is given by 


(7) sin2y = 2B/A. 


There is an ambiguity in this expression, except when y = 45°, and it is for 
this reason that the angular restriction was included in the meteor selection 
criteria. For the present case, the value of B may be deduced from the decay 
rates measured for the two signals. The time taken for the amplitude of 
signal 1 to fall by a factor e is 

7, = 0.54 second, 
and for signal 2 

0.70 second, 


and the ratio 


3/7 = 0.77. 


The height difference that corresponds to this ratio of diffusion times may be 
estimated from the slope of the curve of logarithm of diffusion coefficient 
against height, given by Greenhow and Neufeld (1955). From their results, 
the difference in height of the two reflection points is found to be about 1.7 km. 
Estimating 4 = S/2 = 3 km it is found from (7) above that: 


y = 45° 


as required by the first selection criterion. It should be noted that the validity 
of this deduction depends upon the assumption that the value of 4 is approxi- 
mately S/2; that the meteor trail was near the mid-plane between the trans- 
mitter and the receivers. Since the maxima of the radiation patterns of the 
antennas intersect near this mid-plane, this assumption seems reasonable. The 
possibility that the meteor trail was near the transmitter would appear to be 
ruled out entirely, since, for any smaller value of 4, the right hand side of 
(7) would be much greater than unity. The fact that, for the value chosen, the 
right hand side of (7) is slightly greater than unity may indicate a slight 
displacement of the trail from the mid-plane toward the receivers, but the 
accuracy of measurement is not sufficient to support any such deduction. 

Turning now to the calculation of electron density, and using the subscripts 
1 and 2 to indicate the electron densities at the reflection points that correspond 
to the signals 1 and 2, from (5) 


gi: = 2.41X10" is 2.4110" G2) = 2.1410" /meter 


0.54 


71 


0.25 
0.70 
In practice there is some doubt concerning the initiation time of the signals, 
from which 7 should be measured. The finite rise time of the signal corresponds 
presumably to the time taken by the meteor to traverse the principal Fresnel 


qo = 2.4110" ( ) = 0.86 X 10'*/meter. 
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zone. In this analysis the initiation time of each signal is chosen at a point 
on the record near the top of the initial sharp rise. 

From equation (6) another pair of estimates of electron density may be 
made. By projecting the decay curves backward in time to the initiation time 
of the signal as indicated above, values of A, are found which lead to 


A 1.33 4 18 1.33 
qa = ar :) = 0.7510" i = 3.2210"*/meter 


and 


in = 1.3110" /meter. 


q2 = 0.75x10"( 


The final method of estimating the electron density makes use of equation 
(1), for which it is necessary to deduce the appropriate values for a, B, and ¢. 
A meteor with Y = 45° must occur approximately 100 km off to the side of the 
transmitter—receiver axis in order to satisfy the condition of specular reflection. 
Under the assumptions used earlier this leads to the approximate values 
a = 65° and ¢ = 70°. In an exact treatment, equation (1) should be corrected 
also for the polarization of the receiving antenna. In the present instance this 
correction is small. 

The estimate of 8 is obtained by consideration of the meteor velocity. The 
difference in height of the two reflection points has been estimated at 1.7 km. 
The difference in initiation times for the two signals was 0.16 second, which 
leads to a vertical component of velocity of about 11 km/sec. The component of 
velocity in the plane which is perpendicular to the transmitter—receiver axis 
was, therefore, about 15 km/sec. The results of radar studies of meteor 
velocities (see, for example, Lovell 1954) indicate that in the neighborhood of 
the earth meteors travel with heliocentric velocities close to the parabolic 
limit of 42.1 km/sec. The velocity of the earth in its orbit is 29.8 km/sec 
so that with respect to an observer on the earth the meteor velocity is a simple 
function of the angle between the radiant point of the meteor and the apex of 
the earth’s way, and it has a maximum value of about 72 km/sec when this 
angle is zero. 

In the present circumstances the time of observation (1415 hours local 
time), the time of year (16 March), the direction of the path (nearly east-west), 
the latitude of the path (45° N), and the value of y (45°) serve to fix the 
angular relationship between the half-plane that contains the possible direc- 
tions of the meteor trail (i.e., the half-plane in which 8 is measured) and the 
apex of the earth’s way. Since each direction in this half-plane makes a definite 
angle with the apex direction, it is possible to determine the meteor velocity 
as a function of 8. In the present case these computations are greatly simplified 
because the plane in which @ is measured is parallel to the equatorial plane of 
the earth. The velocities deduced in this way must also be corrected for the 
gravitational pull of the earth. Finally, it is necessary to determine which of 
the permitted combinations of velocity and angle 8 could produce the observed 
component of velocity when projected into the direction 8 = 90°. These com- 
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putations have been carried out using graphical methods and they result in 
the following solutions: 


B = 16.5°, V = 53 km/sec (meteor travelling toward transmitter) ; 


B 55 / 17.5 km/sec (meteor travelling toward receiver). 


In order to choose between the two solutions given in the last paragraph 
it is necessary to consider the height at which the trail was formed. The results 
of an experimental determination of the variation of height of maximum 
ionization with velocity have been presented by Kaiser (1953). From Fig. 16 
of Kaiser’s paper it may be deduced that the height of maximum ionization 
for a meteor velocity of 53 km/sec is about 99.5 km and that for 17.5 km/sec 
is about 86 km. According to Greenhow and Neufeld (1955) the diffusion 
coefficients for the two heights are respectively 8.9 and 1.25 m?/sec. Taking 
the previously deduced value of ¢ (70°), the observed decay rate for signal 1 
corresponds to a diffusion coefficient of 5.8 m*/sec or a height of about 96 km. 
The relative amplitudes of the two recorded signals suggest that the point of 
maximum ionization on the trail was above, rather than below the reflection 
points, which leads to a consistent picture only if the higher value of velocity 
is accepted. It is concluded therefore that the appropriate value of 8 is 16.5°. 

The remaining numerical quantities required for a solution of equation 


(1) are: 


A = 7.65 m, 


2 


Hol = 0.885X10-"'m, 
4m 


R, = Re = 4.5X%106* m, 

Pr 100 w, 

Gr = Gr = 12 (this is the numerical antenna gain in the appropriate 
direction reduced by the amount of the cable losses), 

Pr = A,?/50 (where A, is expressed in volts and Pz, is in watts), 

Ay 4.18 10-* v for signal 1 and 1.72 X10-* v for signal 2, 

hence 


QQ = 2.2X10'*/meter, 
g2 = 0.89 10"/meter. 


The results of the three methods of determining the electron density are 


summarized in Table I. 


TABLE I 
ELECTRON LINE DENSITY AT TWO POINTS ON THE TRAIL AS DETERMINED BY VARIOUS METHODS 


Determination method Equation (1) Equation (5) Equation (6) 


g:(electrons/meter X 107"4) 2.9 2.1 3.2 
g2(electrons/meter X 10~"*) 0.89 0.86 1.3 
Ratio 9:/q2 2.5 2.5 2.5 
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V. DISCUSSION OF RESULTS 

The degree to which the results presented in Table I can be said to support 
the theoretical expressions quoted in Section III is, of course, dependent upon 
the nature of the modifications to these expressions that are contemplated. 
In what follows it is assumed that the validity of equations (1) and (2) is not 
in doubt, or more generally it might be said that the theoretical limitations 
of these expressions are assumed to produce errors that, in the present instance, 
are much smaller than the errors of measurement. An examination of Figs. 2, 3, 
and 4 suggests that the accuracy of the measured quantities cannot be better 
than a few per cent. On the other hand, it is doubtful that there are any very 
large errors of measurement in view of the remarkable agreement among the 
ratios listed in the last line of Table I. In addition, there are the errors that 
result both from the extrapolation of measured values and from the un- 
certainty about the correct initiation time in deriving the values for A, and 7. 
Together, these errors probably lead to a total uncertainty of the order of 
50% in the derived values of q. 

Table I indicates that equation (5) produces substantially the same values 
of electron density as does equation (1). Equation (5) is derived from equations 
(2) and (4), each of which contain the geometric factor, sec?y. The appearance 
of this factor in equation (2) is not in question but its appearance in (4) has 
been questioned most recently by Manning, who has suggested that, if a factor 
that involves secg is to be included, then the exponent should be no greater 
than about 0.5. In the present instance this would lead to an increase by a 
factor of 5 in the values of g deduced from (5). It does not seem likely that the 
experimental uncertainties would permit so large a change. 

The values of q in the last column of Table I were derived from equation (6) 
which in turn was derived from equations (1) and (3). Once again, the 
geometric factors that appear in (1) are not in doubt. The ratio of the geometric 
factors that appear in (3) and (1) would appear in (6) raised to the power 0.66. 
If it is assumed that in the overdense case there is no forward enhancement 
of the signal and that the geometric factor should be unity, then the values of 
q should be about one-third of those given in the last column of Table I. Here 
again it is not likely that so large a modification is acceptable, but an angular 
dependence which is somewhat more moderate than that of (3) cannot be 
precluded on the basis of the present measurements. Any change in the geo- 
metric factors of (3) would lead to a new value of go, the “transition” density. 
The value used is the same as that which applies in the case of back-scatter 
and the present results suggest that no great error is introduced thereby. 

A final check on the consistency of the results is provided by the ratio of the 
electron densities found for the two reflection points. The theory of the 
variation of electron density along the trail as a function of height has been 
reviewed by Kaiser (1953). A suitable value.for the scale height, /7, at meteoric 
levels is 7 km. From the earlier analysis, the two reflection points are located 
below the maximum of ionization and are separated by 0.25H/. From Kaiser’s 
treatment, for two points with this separation a ratio of 2.5 in electron 
density occurs only when the two points are located about 0.5/7 and 0.75/17 
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below the height of maximum ionization (displacement of these points by 
0.1H in either direction leads to a variation of more than 0.5 in the electron 
density ratio). It has been deduced already that the reflection point for signal 
1 was located at a height of about 96 km. By this argument the height of 
maximum ionization is found to be about 99.5 km. As indicated earlier, the 
appropriate height of maximum ionization for the deduced velocity of 53 
km/sec is also about 99.5 km. The agreement is most satisfactory and serves to 
verify the argument which led to the angular values in the analysis. 


VI. CONCLUDING REMARKS 


The application of theoretical formulae to the signals scattered from a single 
meteor trail cannot positively verify those formulae; however, a quantitative 
description of the meteor trail has been built up through the use of the equations. 
This description shows a significant degree of internal consistency and also is 
in good agreement with the known properties of meteor trails. The manner 
in which each of the various expressions has been used to derive the properties 
of a particular trail is similar to the way in which these same expressions 
would be applied separately in statistical studies. In presenting the formula 
for the forward-scattered signal from an overdense trail (Hines and Forsyth 
1957) it was pointed out that the expression could hardly be exact, but might 
prove useful as a working relation in the prediction and interpretation of 
forward-scatter data. The success of the present analysis would appear to 
justify this hope. Certainly, no similar simple relation which would provide 
a better basis for the analysis of the present data appears to be available. 

In closing it may be appropriate to point out another limitation of the 
present analysis. It is now well known that after formation a meteor trail 
tends to be distorted sinuously by the action of large-scale turbulence in the 
upper atmosphere. After this distortion has proceeded for some time the 
scattered radio energy is no longer concentrated in well-defined directions 
and the amplitude of the signal, as received at a fixed point, fluctuates. This 
fading is characteristic of the signals received from distorted trails and such 
trails were excluded specifically from the analysis of this paper. For those statis- 
tical studies that do not exclude the fading signals, it is necessary to consider 
a further treatment to take account of the distorted trails. While fading 
occurs both in back-scattered and in forward-scattered signals, some forward- 
scattered signals last longer than do the corresponding back-scattered signals, 
so that distortion has a proportionately greater chance of becoming effective. 
The fact that distortion sometimes modifies the reflecting properties of the 
meteor trails does not invalidate the present analysis, but it does indicate 
that further analysis would be necessary in any comprehensive treatment of 


the statistical problem. 
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LETTERS TO THE EDITOR 


Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of tssue. No proof will be sent to the authors. 


The Kinematic Viscosity of Liquid Helium II 


Previous studies with a rotating cylinder viscometer in liquid helium II (Heikkila and 
Hallett 1955; Woods and Hallett 1958) have been concerned with the measurement of the 
steady deflections of the inner cylinder while the outer rotated at constant angular velocity; 
such studies yield values of the dynamic viscosity, 7, of the liquid. 

Also of interest are the transient motions of the inner cylinder which occur when the rotation 
of the outer cylinder is started or stopped, because the time, constant for the decay of such 
transients, depends on the kinematic viscosity, » = n/p, of the liquid. Since in liquid helium II, 
v might be determined either by p, the total density of the liquid (as suggested by Donnelly 
1958), or by pn, the density of the normal component, such studies are of considerable interest. 

The apparatus used differed from that described by Woods and Hallett (1958) in two 
respects. A new outer cylinder was made giving a gap width, d, of 2.5 mm, thereby giving 
transient times of convenient duration. Also the R.F. system detecting the motion of the 
inner cylinder was replaced by the standard lamp, scale, and mirror system using windows 
constructed in the dewar walls. 

Experimentally the deflection, ¢(¢), of the inner cylinder was measured as a function of 
time, ¢, from the starting and stopping of the outer cylinder. The appropriate transient solution 
of the Navier-Stokes equation gives the theoretical variation of @ with vt rather than t. Com- 
parison between the experimental and theoretical transients can be made most easily if, for 
selected values of ¢(t)/¢( © ) (¢( ©) being the steady state deflection), the corresponding experi- 
mental times, ¢, and theoretical values, vt, are plotted against each other. The resulting curve 
should be a straight line whose slope gives the effective value of »v. A typical curve obtained 
in this way at 2.16° K (1955 L scale) where pa = p, is shown in Fig. 1 and has the form of the 
expected straight line. From the slope of this line » = 1.6 X 10~* stokes in good agreement 
with n/p = 1.67 X 1074 stokes. Thus the theoretical analysis fits the experimental observations 
reasonably well. The Reynolds number of the motion, R = vd/v, where v is the linear velocity 
of the outer cylinder, was 340 in this case. 

At 0.925° K, however, the transient times, both for starting and for stopping, were much 
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less than 16.4 seconds, the free period of oscillation of the inner cylinder. Since at this tempera- 
ture 7 = 55+2 micropoise, it would be expected that the time required for ¢(t)/¢( © ) to reach 
0.5 would be about 25 seconds if vy = n/p, or about 0.1 second if vy = n/pn. Thus it can be 
seen that the effective density determining v is much less than p and is possibly pr. These 
experiments were carried out at Reynolds numbers < 230, the critical value for this apparatus. 

Between 0.925° K and 2.16° K the viscosity of helium II is much smaller and, in order to 
obtain measurable deflections, it was necessary to exceed the critical velocity and use Reynolds 
numbers of the order of 900. The results showed that the transient times for starting were 
much less than 16 seconds suggesting a very small value for the effective density; the transient 
times for stopping were, however, longer and easily measured, but (t)/o( ©) varied with 
time in a manner different from that predicted theoretically. It is possible that these stopping 
transients could be described by assuming that the effective density determining v increases 
from a value of the order of p,z to a value approaching p as the Reynolds number increases. 
Such a at high Reynolds numbers would then be in agreement with Donnelly’s result 
that vy = 7 
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